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INTRODUCTION 

The  anniversary  conference,  celebrating  the  100th  anniversary  of 
V.  I.  Lenin’s  birth,  devoted  to  the  topic  "Adaptive  Systems"  took 
place  at  the  Moscow  Aviation  Institute  on  March  26-27,  1970.  The 
conference  was  organized  by  the  Chair  of  "Aircraft  Automatic  Control 
Systems"  at  the  Institute. 

The  conference  attracted  a  large  number  of  specialists  in  the 
area  of  adaptive  systems.  It  was  attended  by  professors,  graduate 
students,  engineers  associated  with  the  Laboratory  Chair,  and  repre¬ 
sentatives  of  a  number  of  scientific  research  organizations. 

The  present  edition  contains  lectures  and  communications  that 
were  delivered  at  the  conference.  They  deal  with  various  aspects  of 
the  eybex’netlcs  of  aircraft. 


Zh.  S.  Ageyev,  B.  V.  Viktorov,  B.  V.  Kirsanov,  N.  A.  Shokolo 
took  part  in  the  preparation  and  the  editing  of  the  present  volume. 


FTD-HC-2 3- 1018-72 


V 


INTRODUCTORY  GREETING 
Academician  B.  N.  Petrov 


Comrades! 

Our  country  celebrates  a  momentous  anniversary.  One  hundred  years 
have  passed  from  the  birth  of  V.  I.  Lenin,  the  founder  of  our  party, 
the  creator  of  the  world's  first  socialist  nation. 

V.  I.  Lenin  made  an  enormous  contribution  to  the  development  of 
Marxism,  and  to  the  development  of  social  and  natural  sciences.  A 
clear  illustration  of  Lenin's  deep  foresight  is  Included  in  his 
famous  statement:  "The  electron  is  also  inexhaustible  like  the  atom, 
nature  is  infinite."* 

Lenin  b^ilt  the  foundation  for  the  scientific  management  of 
government  and  social  development.  He  devoted  a  great  deal  of  atten¬ 
tion  to  scientific  and  technological  progress,  and  to  the  scientific 
organization  of  work.  He  emphasized  that  the  proletariat  should  take 
advantage  of  the  best  elements  in  what  has  been  created  during  the 
centuries  of  the  development  of  man's  civilization:  "...  Everything 
that  has  been  conquered  by  science,  technology,  all  the  improvements, 
the  entire  knowledge  of  specialists,  all  this  should  serve  the  united 
proletariat."** 

*V.  I.  Lenin,  Complete  works,  Vol.  38,  p.  26. 

**V.  I.  Lenin,  Complete  works,  Vol.  18,  p.  275. 
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In  an  outline  of  the  plan  of  scientific  and  technological 
development,  written  by  Lenin  in  April  of  1918,  he  placed  a  whole 
series  of  the  most  important  tasks  before  the  members  of  the  Academy 
of  Sciences. 

Lenin's  GOELRO*plan  was  born  in  a  difficult  period  for  the  young 
Soviet  state.  It  was,  however,  not  only  the  first  state  economic  plan 
designed  for  a  number  of  years,  but  also  the  first  state  plan  of  scien¬ 
tific  research.  The  plan  was  full  of  optimism.  It  was  based  on  a 
deep  belief  in  the  primary  role  of  electrification  as  a  basis  of 
technological  progress  and  development  of  the  national  economy. 

Our  nation,  guided  by  the  Communist  Party,  has  been  following 
Lenin's  path  for  over  50  years  already.  The  first  socialist  state 
created  a  powerful  Industry,  leading  technology,  and  advanced  agri¬ 
culture.  Soviet  science  has  flourished  as  never  before,  and  in  a 
number  of  important  fields  it  has  taken  the  leading  role  in  the  world. 

The  Soviet  people  made  a  great  contribution  to  the  progress  of  modern 
science  and  technology.  Our  nation  showed  man  the  way  to  the  universe, 
and  the  way  to  a  peaceful  use  of  atomic  energy.  Our  country  is  a 
leading  power  in  aviation. 

The  development  of  aerospace  technology  in  recent  years  has 
introduced  a  whole  series  of  completely  new  problems  in  the  theory 
and  technology  of  aircraft  control.  The  traditional  methods  of  automatic 
flight  and  attitude  control  of  aircraft  are  no  longer  sufficient  to 
satisfy  the  increasing  demands  on  control  systems.  In  order  to  solve 
the  control  problems  arising  in  that  area,  it  is  necessary  to  use  ideas 
and  methods  that  one  can  call  "cybernetic"  in  the  t':dl  sense  of  the 
word.  The  development  of  microelectronics  and  computers,  the  creation 
of  control  devices  enables  us  to  realize  effective  control  algorithms 
and  to  create  automatic  control  systems  capable  of  controlling  both  an 
individual  aircraft  and  their  complicated  combinations. 

The  wide  range  of  flight  conditions,  the  great  range  of  effectiveness 

of  control  surfaces,  and  the  increasingly  complicated  tasks  facing  air- 

craft,  gave  rise  to  the  emergence  of  a  new  class  of  control  systems, 

"Translators’  Note:  This  designates  State  Commission  for  Electrification 
of  Russia. 
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namely  adaptive  systems).  The  latter  are  capable  of  providing  the 
required  control  quality  under  changing  flight  conditions,  with  the 
parameters  of  the  controlled  objects  varying  within  a  wide  range. 

Such  systems,  as  we  know,  include  self-adjusting  automatic  pilots, 
systems  of  variable  structure,  parametrically  invariant  systems,  and 
others.  The  theory  of  terminal  control  has  progressed  a  great  deal. 

The  foundations  of  the  theory  of  a  new  class  of  systems,  systems  of 
coordinate-parametric  control,  are  being  created.  Control  theory  makes 
effective  use  of  the  methods  of  game  theory  and  statistical  decision 
making,  operations  research,  theory  of  mass  servicing,  methods  of 
linear  and  nonlinear  programming,  and  many  other  theories  Included  in 
modern  theory  of  large  systems. 

The  theory  of  optimal  control  has  developed  to  a  great  extent. 

The  idea  of  optimization  is  the  fundamental  and  decisive  concept  in 
the  design  of  aircraft  and  power  plant  control  systems. 

This  conference  was  arranged  in  order  to  provide  a  comprehensive 
view  of  some  aspects  of  the  work  being  done  in  the  area  of  the  theory 
and  methods  of  design  of  adaptive  contro..  systems,  and  to  outline  the 
fundamental  objectives  in  further  development  of  this  important  class 
of  modern  aircraft  control  systems. 

It  is  my  pleasure  now  to  open  this  anniversary  conference  of  the 
Chair  of  Automatic  Control  Systems  associated  with  the  Moscow  S. 
Ordzhonikidze  Aviation  Institute,  devoted  to  the  theory  and  methods 
of  adaptive  control,  and  dedicated  to  the  100th  anniversary  of  V.  I. 
Lenin's  birth. 
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IMPLEMENTATION  OP  ADAPTIVE  CONTROL  ALGORITHMS 
BY  MEANS  OP  ON-BOARD  DIGITAL  COMPUTERS 

Alekseyev,  K.  B.j  Teryayev,  Ye.  D.j 
and  Ukolov,  I.  S. 

The  solution  of  present-day  control  problems  requires  the  use  of 
digital  computers  in  control  systems.  At  the  present  time  the  problem 
can  be  solved  If  digital  computers  are  used. 

The  use  of  computers  in  control  systems  assumes  particular 
importance  and  timeliness  in  connection  with  the  technical  execution 
of  adaptive  (self-adjusting)  control  systems.  In  this  area  one  is 
faced  with  the  specific  problems  of  studying  the  dynamics  of  automatic 
control  systems  that  include  computers,  the  problem  of  determining  the 
basic  characteristics  of  on-board  digital  computers,  and  their  inter¬ 
action  with  the  dynamic  characteristics  of  aircraft. 

In  the  present  report,  we  discuss  certain  problems  that  arise  with 
adaptive  control  algorithms  in  automatic  control  systems. 

An  adaptive  control  algorithm  will  be  defined  as  a  changing 
program  of  operations  that  takes  into  account  the  actual  state  of  the 
surrounding  medium,  the  parameters  of  the  input  and  output  coordinates 
of  the  object,  and  optimizes  in  a  certain  sense  the  value  of  a  given 
functional  of  the  quality  of  the  control  system. 
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Similarly,  an  adaptive  control  system  will  be  defined  as  a  system 
that  realizes  an  adaptive  algorithm  in  the  sense  described  above. 

On-board  digital  computers  offer  a  great  deal  of  promise  as  far 
as  their  application  in  aircraft  control  systems  is  concerned.  The 
effectiveness  of  using  on-board  computers  in  control  systems  depends 
on  whether  it  will  be  possible  to  use  radically  new  methods  of  solving 
the  control  problems,  methods  that  would  result  in  better  performance 
parameters  for  a  control  system  as  a  whole. 

The  use  of  on-board  digital  computers  is  based  on  the  following 
considerations: 

a.  Great  logic  and  computational  capabilities,  permitting  the 
solution  of  diverse  and  complicated  problems  with  practically  unlimited 
accuracy; 


b.  Stability  of  operation  within  a  wide  range  of  working 
conditions; 

c.  High  speed,  permitting  the  simultaneous  solution  of  several 
control  problems; 

d.  High  level  of  adaptation  and  self-control  of  the  control 
systems. 

On-board  digital  computers  offer  great  possibilities  to  the 
designers  of  control  systems  when  they  are  involved  in  designing  com¬ 
plicated  adaptive  control  algorithms.  The  reliability  of  control 
systems  is  always  improved  both  as  a  result  of  a  programmed  correction 
of  errors  in  the  operation  of  certain  reserved  elements  of  the  control 
systems,  and  by  controlling  the  apparatus  directly  during  its  operation. 

In  the  design  of  control  systems  involving  on-board  digital 
computers,  we  are  faced  with  the  following  technical  tasks: 
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1.  Analysis  of  the  structure  of  control  systems  for  the  purpose 
of  finding  those  versions  that  possess  highest  reliability  and  are 
least  costly; 

2.  Determination  of  rational  requirements  on  the  on-board 
computers; 


3*  Organization  of  efficient  parallel  programming  of  individual 
problems  solved  with  the  aid  of  an  on-board  computer; 

4.  Execution  of  the  on-board  computer  programs  in  real  time, 
while  satisfying  the  requirement  that  different  levels  of  reliability 
be  achieved  when  solving  problems  of  various  degrees  of  importance 
(l.e.,  various  degrees  of  Influence  on  flight  safety),  and  while 
satisfying  the  time  requirements  of  various  computer  users; 

5.  Taking  advantage  of  the  multichannel  structure  of  a  control 
system  and  an  on-board  computer  in  order  to  use  programming  to  Increase 
the  reliability  of  the  most  Important  oontrol  algorithms; 

6.  Debugging  computer  programs  using  general-purpose  digital 
computers  and  simulation  programs; 

7.  Debugging  computer  programs  using  an  analog-digital  (hybrid) 

unit. 


An  on-board  digital  computer  used  as  a  component  of  aircraft 
control  systems  must  perform  a  large  number  of  operations  needed  to 
solve  the  following  functional  problems: 

1.  Navigation  and  control  of  aircraft  in  various  stages  of 
flight; 

2.  Collection  and  processing  of  information  about  the 
surrounding  medium; 
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3.  Solution  of  auxiliary  problems;  control  of  the  operating 
conditions  of  individual  units,  selection  of  optimal  flight  conditions, 
etc. 


The  basic  feature  of  control  systems  and  control  based  on  the 
use  of  on-board  computers  is  the  possibility  of  flexible,  adaptive  and 
automatic  changes  to  the  control  program  and  the  possibility  of  making 
decisions  based  on  an  assessment  of  the  flight  situation. 

On-board  digital  computers  may  be  classified  according  to 
different  basic  characteristics  which  assign  them  to  a  given  class. 

One  possible  classification  is  based  on  the  location  and  purpose 
of  the  on-board  digital  computer  within  the  control  system: 

a.  Command  computers,  performing  strict  control  according  to 
a  preset  program,  and  independently  of  the  course  of  the  controlled 
process; 


b.  Compensating  computers,  used  in  loops  of  strict  load  control; 

c.  Information  computers,  where  the  information  about  the  state 
of  a  system  is  processed,  received  by  man,  and  finally  used  to  control 
the  system; 

d.  Control  computers,  used  in  a  closed  control  loop.  This 
form  of  application  is  basic  and  most  promising.  Their  location  in 
the  control  system  and  their  tasks  may  be  quite  diverse. 

The  use  of  on-board  digital  computers  in  control  systems  poses 
particularly  great  demands  on  such  parameters  of  a  computing  system 
as  weight,  volume,  reliability,  required  power,  range  of  operating 
temperatures,  radiation  resistance,  transmission  of  vibrations,  etc. 

In  addition,  the  on-board  computers  must  satisfy  a  number  of  specific 
demands.  A  computerized  control  system  should:  1)  operate  reliably 
without  preventive  maintenance  and  servicing  during  the  entire  period 
of  effective  service,  i.e.,  during  the  flight  of  an  aircraft.  This 
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means  that  an  on-board  computer,  consisting  of  thousands  of  components, 
must  have  an  acceptable  service  life  (In  the  case  of  an  on-board  com¬ 
puter  which  cannot  self-repair,  It  Is  equal  to  the  time  Interval  up  to 
the  first  failure),  2)  the  range  of  problems  solved  by  the  computer 
will  usually  remain  the  same  during  the  entire  service  life  of  the 
machine,  3)  Input  and  output  operations  must  be  carried  out  In  real 
time,  and  the  time  Intervals  must  be  strictly  specified.  An  on-board 
computer  must  be  capable  of  being  used  for  a  wide  range  of  problems, 
since  It  may  become  necessary  to  change  the  scale  of  the  telemetry 
data,  change  the  calculations  needed  to  prepare  the  Initial  data,  etc. 

When  synthesizing  a  computerized  control  system,  the  following 
parameters  are  given: 

1.  Program  running  time; 

2.  Speed,  defined  as  the  time  needed  to  perform  a  short 
operation; 

3.  System  of  commands; 

J».  Bit  size  of  numbers; 

5.  Storage  capacity. 

Program  running  time  Is  determined  by  the  maximum  permissible 
length  of  the  elementary  tije  Interval  for  the  most  dynamic  control 
algorithm  realized  using  the  on-board  oomputer. 

Programs  for  on-board  computers  have  a  complicated  structure  and 
consist  of  subprograms  that  perform  Individual  control  tasks,  and 
of  standard  subprograms  that  perform  typical  algorithms  used  In  control 
and  processing  of  numerical  Information;  a  subprogram  which  adjusts 
the  special  and  the  standard  subprograms  In  accordance  with  data 
(obtained  from  the  outside)  on  the  working  conditions,  the  time,  and 
external  control  signals,  and  an  Initializing  subprogram  that  reads 
the  initial  codes  into  the  operative  memory. 
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The  type  of  on-board  digital  computer  is  determined  by: 

1.  Complexity  of  algorithms; 

2.  Speed; 

3*  Accuracy  required  in  algorithms; 

4.  Allowable  weight  and  size; 

5.  Ease  of  adjustment,  so  that  other  problems  can  be  solved; 

6.  Freedom  from  interference; 

7*  Possibility  of  controlling  the  control  system. 

The  demands  placed  on  on-board  computers,  as  well  as  on  their 
structure  and  the  structure  of  the  input  and  output  units,  depend  on 
the  specific  problem  and  algorithm. 

Depending  on  the  type  and  purpose  of  an  aircraft,  an  on-board 
compi'ter  may  perform  the  following  tasks: 

1.  Data  processing; 

2.  Computation  of  controlled  quantities; 

3.  Generation  of  control  signals; 

4.  Control  functions. 

The  tasks  in  this  group,  in  particular  the  navigational  algorithms, 
require  a  higher  accuracy  of  calculations.  In  this  case,  the  computer 
word  length  should  be  no  less  than  20  -  24  bits,  and  the  speed  —  at 
least  5-8  thousand  operations  per  second. 

Cmtrol  of  the  on-board  equipment  by  means  of  a  computer  consists 
of  the  algorithms  of  automatic  control  and  those  for  predicting  defects 
*■’  .t  can  occur  while  the  basic  problems  are  being  solved.  In  order  to 
solve  the  control  problems,  one  needs  Increased  capacity  of  the  fixed 
and  operative  computer  memories. 

Certain  characteristics  of  foreign  on-board  third-generation 
computers  are  given  in  the  table  on  the  following  page. 

Normally,  an  on-board  computer  should  guarantee  the  realization 
of  a  large  number  of  algorithms,  and  thus  it  is  important  to  first 
consider  those  algorithms  that  constitute  the  largest  computational 
load  (n): 
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where  A  and  L  are  the  numbers  of  simple  arithmetical  and  logical 
operations. 

Those  algorithms  include  primarily  the  algorithms  used  In  data 
processing  and  the  computation  of  controlled  quantities. 

The  structure  of  the  on-board  computer  should  be  maximally  adapted 
to  the  composition  of  the  on-board  systems,  the  problems  being  solved, 
and  should  be  optimal  In  the  sense  defined  in  the  design  criterion. 

It  is  clear  that  it  is  necessary  to  have  a  certain  number  of  possible 
structural  versions. 

Unfortunately,  at  the  present  time  we  do  not  know  of  an  established 
procedure  for  solving  this  problem,  even  though  there  are  some 
approaches  that  permit  us  to  formalize  the  problem. 

The  increasing  demands  as  to  the  minimization  of  energy  losses 
incurred  in  orienting  flight  vehicles  on  trajectories  with  a  long 
flight  time  fan  the  order  of  one  year)  force  us  to  seek  new,  more 
economical  methods  of  controlling  the  motion  of  the  vehicle  about  its 
center  of  mass.  One  of  these  methods  involves  extensive  control  of 
the  attitude,  and  can  be  realized  by  using  a  computerized  control 
system.  The  method  consists  of  selecting  the  moments  applied  to  the 
vehicle  relative  to  the  axes  attached  to  it  in  such  a  way  that  the 
vehicle  is  set  in  motion  about  the  axis  of  equivalent  rotation  (the 
Euler  axis). 

As  compared  with  existing  methods;,  the  method  of  extensive  control 
enables  us,  other  conditions  being  equal,  to  solve  the  basic  problems 
of  orientation  with  much  smaller  energy  losses  (~  by  a  factor  of  2) 
and  much  more  rapidly  (•.  by  a  factor  of  2  -  3).  However,  extensive 
control  involves  a  relatively  large  number  of  computations,  which 
assumes  that  an  on-board  computer  is  used.  It  is  characteristic  that 
the  use  of  an  on-board  computer  in  systems  of  extensive  control  is  a 
necessary  condition  for  their  technical  implementation.  Moreover,  the 
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advisability  of  the  use  of  an  on-board  computer  Is  justified  by 
significant  energy  gains  and  higher  speed,  which  permits  us  to  lengthen 
the  active  life  of  a  vehicle,  and  Improve  its  tactical  and  engineering 
characteristics  without  additional  expenditures. 

Let  us  consider  the  principal  problems  involved  In  the  construction 
of  computerized  systems  of  extensive  control,  focussing  our  attention 
on  the  physical  oontent  of  the  problem  and  the  computer  functions 
resulting  from  It. 

The  formulation  of  the  control  problem  Is  as  follows.  Suppose 
we  are  given  the  Initial  (t  ■  0)  and  the  final  (t  ■  T)  angular  atti¬ 
tudes  of  a  vehicle,  determined  by  means  of  modified  Euler  angles  T,  v 
and  y,  subject  to  a  given  set  of  constraints. 

We  must  determine  a  program  of  the  time  dependence  on  the  bounded 
control  moments  Mi,  M«,  M*  applied  to  the  vehicle  relative  to  the  axes 
attached  to  It,  which  are  such  that  the  vehicle  Is  taken  from  Its 
initial  attitude  to  the  final  one  during: 

1.  a  minimum  time  T  , 

nt 

2.  a  given  time  T  with  minimum  energy,  losses. 

We  assume  that  the  vehicle  is  an  Ideal  rigid  body,  and  the  effect 
oi"  external  perturbations  on  the  angular  motion  of  the  vehicle  In  the 
course  of  control  Is  negligibly  small. 


Figure  1-2. 
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Without  giving  detailed  calculations,  we  shall  only  state  that 
the  problem  of  extensive  control  has  an  exact  analytic  solution. 

This  fact  is  of  great  importance  in  constructing  the  computer  functions 


Let  us  analyze  the  sequence  of  operations  performed  by  the  on¬ 
board  computer.  The  geometry  of  the  problem  is  shown  in  Figure  2. 

We  shall  give  the  computational  operations  necessary  to  determine  the 
resulting  angle  of  turn  ♦  ,  direction  of  the  Euler  axis  l »  [V|,  v2l  v* 
(in  the  system  Oxj,  xa,  x*  and  OXi,  Xa,  X»),  and  the  sign  of  the  angle 
of  turn,  sgn  *.  Using  the  initial  measured  values  of  <p,  d  and  y,  we 
determine  the  matrix  of  direction  cosines  a  *  ||ay  ■  || ;  t,  o'  *  1,  2,  3, 
which  we  then  use  to  calculate 


f  -  arc  cos  ~  («„  +  «a  +  aM  -  1) 
V[  - 

|«»  v,  j|  ay-i 


(2) 


The  last  operation  sgn  4  involves  the  largest  amount  of  labor. 


gyroscopic  moment. 


Figure  3. 
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Knowledge  of  the  direction  of  the  Euler  axis  gives  us  the 
direction  of  the  kinetic  moment  vector  and  the  axis  of  the  gyroscopic 
moment  (Figure  3). 


The  directions  of  these  vectors  can  be  found  from 

j 


1, 


If 


At- 


|  /.  o  o 
0  /,  0  ; 
I  0  0  /,  I 


then 


where 


1, 


VrH  ~  i+l)5  vl+ly*M 


(3) 


(4) 


(5) 


Here  t  are  the  unit  vectors  of  the  attached  coordinate  system.  We 

B  j 

note  that  the  unit  vector  of  the  Euler  axis  is  »  „ 

Zl  Vr 
o-i 

Knowledge  of  the  orientation  of  l and  lr  enables  us  to 
specify  the  direction  of  of  the  axis  of  the  resultant 

moment  vector  M,  and  write  implicit  expressions  (in  terms  of  the 
angle  0)  for  its  direction  cosines  (Figure  4). 


P«"“(a»  +  4#tg0)cos9 
a*  -  7J5 :  I  aH  |  -  aj 

K-lfcU 
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(6) 


Denoting  the  maximum  allowable  value  of  the  moment  applied  to 
the  vehicle  along  an  attached  axis  by  Mms  and  the  maximum  value  of  a 
direction  cosine  by  ft,  we  find  that  the  maximum  length  of  the  vector 
M  is  equal  to 

(7) 

Then  the  equations  of  motion  of  the  vehicle  can  be  written  using  the 
condition  that  the  dynamic  ^eection  of  the  vehicle  to  the  applied 
rontrol  moment  be  compensated:  The  equations  of  motion  are: 


Jocose 


Mnu 

/*  |««±Mg6l 


(8) 


au  tge 
77 1  a,±»,  1C»  I 


(9) 


Ti'ip  subsequent  computational  procedure  used  to  synthesize  the  system 
is  given  without  explanations. 


Synthesis  of  the  System 
r 

0)  /  —  |  dt  -*  min 
/.-isrcth/th^*  for  1) 


(10) 

(ID 
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or  /,-iarcthl/th^Il 

*  B  T  Rlml 

(12) 

r.- i[«,cth]/S^‘+.,c .h/tl^jfor  2) 

(13) 

For  0  <  t  <  ti  and  1) 

p* 

UK) 

■»  «,  a, 

mj  a»‘+b»lF,'s*ln, 

Mf+i  -  Mm,  afb!llU 

(15) 

and  2) 

«£l«eT 

1 

*? 

1 

* 

5= 

(16) 

i(  u  SUi+^i^tln*^ 

">rfl "  ’  *«*»«» 

(17) 

For  t4  <  t  <  T  and  2) 

P« 

(18) 

..  •  -/on 

M*+t "  A,«* — — 

(19) 

and  2)  Af,  — Al^k 

P« 

(20) 

A*,+/  -  ^  (SH<4-».wsh,U-B(/-/fn| 

** o,  ch*(x  ~  b  (/ —V,)J 

(21) 

i  “  1,2  x  «■  arctgthn^ 

(22) 
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(23) 


T  3 

Od-t 


(24) 

(25) 


For  ti^t<ti+tj 

M  •»  tliJ  nti 


(26) 


A*i- if*  thJ i- 1,2,3 


(27) 


Figure  5. 

The  existence  of  the  analytic  solution  of  the  problem  involving 
the  control  of  a  spatial  turn  maneuver  makes  it  possible  to  use  an 
on-board  computer  as  a  command  computing  unit  (Figure  5).  Such  a  use 
of  the  computer  excludes  any  feedback  (rigid  program).  However,  this 
does  not  mean  that  feedback  coupling  will  occur  in  the  controlling 
part  of  the  system  and  when  eliminating  small  angular  deviations. 
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»■  *  f  l«X¥3V>c%e» 


When  performing  a  turn  by  <J>  =  180°  the  effect  of  the  external 
disturbances  may  result  in  an  inadmissibly  large  error  in  the  final 
attitude  of  the  vehicle.  Therefore,  it  becomes  necessary  to  introduce 
a  correction,  which  can  be  effected  by  maintaining  a  programmed  angular 
velocity  of  the  vehicle  (Figure  6).  Then  the  computer  functions  are 
somewhat  widened,  but  the  purpose  of  the  computer  remains  the  same. 


Figure  8.  BCDP  function. 


As  for  its  structure,  an  on-board  computer  containing  arithmetic 
and  logic  units  may  be  subdivided  into  three  blocks  which  are  designed 
to  compute: 

1.  Kinematic  parameters; 

2.  Dynamic  parameters; 

3.  Programmed  values  of  controlled  moments. 

A  general  characteristic  of  one  of  these  blocks  is  shown  in 
Figure  7.  We  can  see  that  extensive  control  using  simple  computing 
devices  is  impossible.  The  functions  performed  by  the  other  blocks 
are  shown  in  Figures  8  and  9- 
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Figure  9.  Functions  of  BCCM. 

An  important  question  arises  in  the  control  of  the  attitude  with 
adaptation.  This  type  of  control  occurs  when  the  moments  of  inertia 
of  a  vehicle  either  change  dgnificantly  or  are  to  be  determined  in 
the  course  of  the  flight.  A  study  of  this  problem  shows  that  in  this 
case  it  is  also  advisable  to  retain  the  above-mentioned  function  of 
the  on-board  computer  that  ensures  optimal  control.  However,  there 
appears  a  new  computer  function  which  involves  the  problem  of  identi¬ 
fying  the  object.  It  is  characteristic  that  this  function  is  performed 
by  the  computer  In  a  closed  control  loop. 
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ELEMENTS  OP  ADAPTIVE  AUTOMATIC  CONTROL  SYSTEMS 
Braslavskiy,  D.  A.  and  Yakubovich,  A.  M. 

Adaptive  control  systems  are  constructed  using  both  general-  and 
special-purpose  elements.  The  general-purpose  elements,  used  to  con¬ 
struct  both  the  adaptive  and  linear  control  systems,  are  in  the  form 
of  linear  converters.  They  include  data  sensors  (angle-angular 
velocity,  altitude  sensors,  etc.)  and  amplifying  converters  (operation 
amplifiers,  power  amplifiers,  differentiators,  integrators,  etc.). 

The  special-purpose  elements  are  mainly  used  in  adaptive  control 
systems.  They  Include  the  following  basic  elements: 

-  test  pulse  generators, 

-  narrow-band  filters,  which  select  signals  in  a  given  narrow 
frequency  band, 

-  multipliers, 

-  lntegro-multipliers, 

-  functional  multipliers, 

-  quorum  elements. 

The  basic  requirements  placed  on  the  above  elements  are  given  in 
the  table. 

Pulse  generators  and  narrow-band  filters  can  be  constructed  using 
electronic  amplifiers  and  threshold  elements  in  combination  with  RC 
networks.  Pulse  generators  and  narrow-band  filters  having  a  higher 
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No. 


Element 


Basic  Requirements 


1 


2 


3 


4 


5 

6 


Pulse  Generator 


Narrow-band  filter 


Multipliers 


Integro-multiplier 


Frequency  range  2  -  30  Hz 
Frequency  stability  2  -  5* 

Amplitude  stability  5  -  10% 

Resonance  frequency  2  -  30  Hz 
Frequency  stability  2-5* 

Amplitude  stability  5  -  10* 

Null  level  accuracy  0.2  -  0.5* 
Multiplication  accuracy  3  -  5* 

3  dB  level  frequency  band  100  -  200  Hz 

Integration  time  constant  5  -  50  sec 
Null  level  accuracy  of  input  0.2  -  0.5* 
Multiplication  accuracy  15  -  20* 


Null  level  accuracy  0.2  -  0.5* 

Functional  multiplier  Accuracy  of  functional  conversion  5  -  10* 

3  dB  level  frequency  band  100  -  200  Hz 


Quorum  element 


Accuracy  of  conversion  1-2* 
Threshold  accuracy  5  -  10* 


stability  of  characteristics  (1  -  3*  as  to  frequency,  3  -  6*  as  to 
amplitude)  are  constructed  using  electromechanical  devices  with 
mechanical  moving  systems. 

The  multipliers  are  constructed  using  various  techniques  which 
involve  both  direct  and  indirect  transformation  of  signals.  It  should 
be  noted  that  the  multipliers  for  closed  adaptive  systems  can  be 
constructed  more  simply  than  others,  since  the  requirements  on  the 
accuracy  of  multiplication  can  be  much  lower.  Integro-multlpllers 
which  perform  the  operation  z(t)  ■  y(t)  •  Jx(t)dt,  where  y(t),  x(t) 
are  input  coordinates,  are  used  in  systems  used  with  a  standard  model. 
The  integro-multlpllers  may  be  constructed  according  to  various 
principles.  As  an  example  we  can  mention  an  electrochemical  integro- 
multiplier  which  is  called  a  mlnlstor.  Low  temperature  stability  is 
the  minister's  basic  drawback.  We  shall  consider  two  important  and 
promising  elements  in  more  detail:  the  functional  multiplier  and 
the  quorum  element. 
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Functional  Multipliers 


The  functional  multipliers  form  the  product  of  the  input  coordinates 
y»...yN  entering  through  N  channels  and  of  the  functions  f » (x) . . ,fN(x) 
of  the  coordinate  x,  so  that  the  output  coordinates  of  the  channels , 
zj . . .  z  ,  are  equal  to 

The  most  common  functions  used  in  automatic  control  systems  include 
algebraic  and  transcendental  functions  of  the  forms  x°,  e“ox.  In  x, 

1  -  x°,  (1  -  x)°,  where  a  is  an  integral  or  fractional  positive  number. 

One  of  the  most  promising  versions  of  the  functional  multiplier 
is  the  discrete  dynamic  version  [1,  2]  with  analog  input  and  output, 
and  pulse  signals  used  in  all  Intermediate  transformations. 

A  functional  multiplier  (Figures  1  and  2)  consists  of  a  control 
channel  and  N  channels,  x  is  the  input  to  the  control  channel,  and 
yA  is  the  input  to  the  ith  channel.  The  unit  consists  of  the  following 
subunits:  linear  dynamic  network,  comparator,  pulse  keys,  delay  lines, 
and  memory  capacitors.  Let  us  analyze  the  operation  of  the  unit. 

Suppose  that  in  the  control  tract  in  a  linear  dynamic  network  the 

-i 

process  occurs  involving  a  discharge  of  the  capacitor  C» 

charged  by  the  standard  voltage  Uq  >  x.  At  a  time  Tn  when  the 
equality  i8  satisfied,  the  comparator  sends  a  fixing  pulse. 

This  pulse  goes  through  the  delay  line  and  forms  the  tact  pulse.  The 
tact  pulse  at  a  time  t  closes  the  pulse  switch  through  which  C»  is 
then  charged  by  the  standard  voltage  Uq.  The  process  is  then  repeated. 
The  interval  between  the  tact  pulse  and  the  fixing  pulse  is 

7’„ - fln-i  (2) 

Simultaneously  with  the  charging  of  Ci  by  the  tact  pulses  through 
the  pulse  switches,  the  capacitors  Ci  in  the  linear  dynamic  networks 
LN2j  become  charged  (or  discharged).  The  amplitude  to  which  the  capa- 

f*  Vt 

citor  in  the  i  channel  is  charged  is  proportional  to  y^  The 
transient  processes  W/(0“  yrn, (t)  in  the  channels  are  periodic.  Their 
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start  coincides  with  the  start  of  ,the  process  i^.e*  in  the  control 
channel.  At  a  time  Tn  the  fixing  pulses  connect  the  output  capacitors 
CR  to  LN2i...LN2n>  through  the  switches  PS»,  ...»  LSN>  storing  the 
signal  from  the  output  of  LN2i  in  the  ith  channel  equal  to 

■*i  -  ®i (7'«)  «i (T„)  -  yru, (-  Hni)  ( 3 ) 

where 

(4) 


Figure  1. 

1.  control  channel}  2.  channel  1}  3.  channel  N; 

4.  fixing  pulse;  5.  comparator}  6.  LDN*; 

7.  pulse  switch}  8.  tact  pulse}  9.  delay  line. 

Using  various  transient  processes  u.(t)  in  the  channels,  we 
obtain  various  f1(x).  For  example,  if  t~  ,  then  f^(x)  »  x, 

and  »  yi(x),  for  _ ,4.  /,(*>- 1 -i. 

#LDN  designates  linear  dynamic  network 

24 


X 


0  ~t 


Figure  3. 

*CC  designates 
current  clipper. 


Quorum  Element 


The  quorum  element  [3]  is  an  analog  logic  unit  with  N  input 
•yN  (where  N  is  odd),  and  one  output  y  equals 

y~l/|-fi/»+  •••  +IJm 
7  in  ' 


s  y» 
(5) 


for |y-y,|< A,  |y-yy|>Ay;  j^m+  ,  where  /«>4-A/  and 

are  threshold  values.  In  the  special  case  when  0,  > 


y«meil(yi ...  yN), 

i.e.,  is  a  median  of  the  signals  yi...yN.  According  to  [4],  Equation 
(5)  represents  an  algorithm  which  is  close  to  being  optimal  for  adap¬ 
tive  systems.  A  diagram  of  the  quorum  element  for  N  electric  signals 
is  given  in  Figure  3. 
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The  basic  unit  of  the  quorum  element  is  in  the  form  of  a  current 
clipper.  The  volt-ampere  characteristic  of  the  current  clipper  is 
shown  in  Figure  4.  Let  us  analyze  the  operation  of  the  device  for 
N  *  3,  assuming  for  simplicity  that  the  load  resistance  is  *  «. 

Three  modes  of  operation  are  possible:  averaging,  blocking  the  circuit 
of  one  of  the  signals,  and  generation  of  a  median. 


If  the  input  voltages  ui,  U2,  uj  are  such  that 

| «/  —  |  <  «0.  i-  1.2,3, 

then  the  balancing  currents  in  the  circuits,  Ij,  I2,  Ij,  are  smaller 
than  the  limiting  current  IQ.  In  this  case  the  current  clippers  are 
linear  resistances,  and 


Let  us  consider  the  mode  in  which  circuit  1  is  blocked.  Let 


«i  >  "2  >  «a  ,  and 


|«, -«|>  «o 


(6) 


l«2-«K«o;  \Ui-U\<  Mo. 


(7) 


When  the  inequality  (6)  is  satisfied,  the  current  in  the  circuit 
1  attains  the  maximum  value  In.  This  current  splits  into  identical 
currents  in  circuits  2  and  3,  equal  to  f  .  In  addition  to  -f  ,  the 
channels  2  and  3  carry  an  additional  current  due  to  the  fact  that  U2 
is  different  from  us.  If  the  condition  (7)  is  satisfied,  then  the 
resultant  currents  in  the  circuits  2  and  3  will  not  reach  the  maximum 
value  Iq,  The  dynamic  resistance  of  the  nonlinear  current  clipper  CC» 
is  practically  infinite,  and  the  operating  currents  in  the  elements 
CC2,  CC3  remain  in  the  region  of  small  linear  resistances.  The  cir¬ 
cuit  1  is  blocked.  The  device  averages  the  signals  U2  and  u*,  i.e., 

a~u,+u‘ 

The  mode  in  which  a  median  is  generated  will  occur  f or  lwi  “  ui\  >  % 

«j -  w,|>w0 .  Suppose  that  iti>'h>u3  .  Then  the  balancing 
current  IQ  is  shunted  through  the  circuits  1  and  3»  and  /,«=  —  1,  .  In 
the  second  circuit  there  is  no  current.  Since  CCi  and  CC»  are 
saturated,  the  dynamic  resistance  in  those  circuits  will  be  close  to 
infinity.  Since  the  resistances  in  the  circuit  2  are  small,  u  ■  U2. 
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It  will  be  noted  that  the  quorum  element  has  been  constructed  according 
to  a  generalized  structural  diagram  given  in  [33.  [33  discusses 

circuits  of  quorum  elements  for  the  case  when  input  signals  are  in 
the  form  of  displacements  and  pressures. 

In  conclusion,  we  note  that  a  quorum  element  is  a  new  basic 
element  of  adaptive  automatic  control  systems  that  enables  one  to 
solve  several  independent  problems,  including: 

-  realization  of  the  above-mentioned  majorizing  function, 

-  construction  of  adaptive  systems  of  variable  structure, 

-  construction  of  low-frequency  filters  with  small  phase 
distortions  (quorum  filters). 

Fundamental  Problems  and  Prospects  for  the  Further 
Development  of  Elements  of  Automatic  Control  Systems 

The  principal  trends  in  the  work  in  progress  are: 

1.  The  use  of  new  physical  principles  enabling  one  to  simplify 
elements  and  to  increase  their  reliability.  In  particular,  the  use 
of  optical-electronic  and  electrochemical  converters  in  the  construc¬ 
tion  of  multipliers  and  integro-multipliers  3eems  to  offer  promise. 

2.  Unit  miniaturization  of  the  elements  of  automatic  control 
systems. 

The  miniaturization  has  two  aspects: 

a.  A  reduction  of  the  weight  of  the  on-board  equipment  in 
order  to  improve  the  performance  parameters  of  flight  vehicles; 

b.  An  Increase  of  the  capacity  in  order  to  apply  the  method 
of  functional  redundancy  to  a  greater  extent. 

The  miniaturization  of  purely  electronic  elements  is  based  on 
the  use  of  integral  monocrystal  and  hybrid-film  microelectronic 
circuits. 


The  miniaturization  of  the  electromechanical  elements  is  in  some 
cases  based  on  the  use  of  new  technological  processes  without  changing 
the  construction  of  the  element  itself.  In  other  cases,  one  has  to 
develop  new  designs  of  the  elements,  and  sometimes  even  switch  to  new 
physical  principles. 

3.  Transition  from  analog  control  systems  to  discrete  ones  with 
the  maximum  use  of  universal  or  specialized  computers. 


As  digital  computers  become  more  sophisticated,  and  as  their 
speed  and  reliability  increase  while  the  size  becomes  smaller,  it 
seems  that  analog  systems  will  be  replaced  by  discrete  ones.  With 
a  transition  to  digital  computers,  the  nomenclature  of  the  converting 
elements  used  in  automatic  control  systems  will  be  greatly  reduced, 
but  in  all  cases  we  shall  have  to  retain  the  data  sensors  and  the 
I/O  units. 
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CERTAIN  PROBLEMS  OP  THE  ADAPTIVE  BANK 
CONTROL  OP  A  PLIGHT  VEHICLE 

Ageyev,  Zh.  C.;  Viktorov,  B.  V.  and  Ukolov,  I.  S. 

The  bank  control  of  the  effective  aerodynamic  quality  of  a  flight 
vehicle  as  we  know  presupposes  a  flight  with  a  constant  value  of  the 
equilibrium  angle  of  attack  which  is  achieved  by  proper  centering. 

The  control  of  the  effective  aerodynamic  lift-to-drag  ratio  is 
achieved  by  turning  the  vehicle  about  the  air  velocity  vector  by  a 
given  bank  angle.  In  an  emergency  situation  it  is  advisable  to  change 
the  bank  angle  of  the  vehicle  with  an  angular  velocity  on  the  order 
of  15  1/sec, 

However,  a  number  of  studies  (for  example,  [1]),  have  shown  that 
it  is  advisable  to  control  flights  by  simultaneously  changing  the 
angle  of  attack  and  the  bank  angle.  In  such  a  system,  a  change  in 
the  angle  of  attack  substantially  modifies  the  effectiveness  of  the 
bank  regulator.  In  addition,  there  are  several  well-known  dynamic 
problems  of  controlled  lateral  motion,  for  example,  due  to  the  fact 
that  the  air  velocity  vector  and  the  vector  of  the  controlling  moment 
of  bank  nozzles  do  not  coincide.  These  circumstances  may  necessitate 
the  introduction  of  additional  (including  also  adaptive)  structural 
couplings  in  the  control  system. 

This  article  is  devoted  to  certain  aspects  of  this  problem. 
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Figure  1  is  a  diagram  of  an  aircraft.  The  aerodynamics  of  flow 
over  an  axisymmetrlc  body  Includes  the  peculiarity  that  in  any  motion 
of  the  aircraft  relative  to  its  center  of  mass,  due  to  the  symmetry 
of  the  field  of  aerodynamic  forces,  an  instantaneous  plane  of  aero¬ 
dynamic  symmetry  appears.  This  plane  includes  the  linear  velocity 
vector  V,  the  geometric  axis  of  symmetry,  and  the  total  aerodynamic 
force  vector  R.  It  should  be  noted  that  the  control  of  the  attitude 
of  the  aerodynamic  symmetry  plane  is  the  basic  objective  of  controlling 
the  aerodynamic  lift-drag  ratio  of  aircraft,  by  varying  the  bank  angle. 

The  flight  trajectory  in  the  atmosphere  of  vehicles  of  this  type 
is  characterized  by  very  small  angular  velocities  of  the  linear 
velocity  vector  V.  Therefore,  in  the  study  of  the  angular  motion  of 
the  vehicle  it  is  convenient  to  use  the  semi -velocity  system  XYZ  as 
the  reference  system  for  measurement. 

Let  us  determine  the  attitude  of  a  vehicle  (or  in  other  words,  the 
attitude  of  the  attached  system  XiYjZi)  relative  to  the  following  semi¬ 
velocity  system  of  the  Euler  angles  <,  a  and  $  (Figure  2).  For  this 
combination  of  angles,  we  cleavay  have  the  following  kinematic 
relations 


*  “IhmK  sin  V  -UV.  cos  <?); 

(1) 

V-^.-ctgeK  sin  <?-  u»#1  cos  f); 

(2) 

a— <u,(  cbs  ?+“>#,  sin  a. 

(3) 
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In  writing  the  dynamic  equations  of  angular  motion,  we  make  the 
following  assumptions: 

the  inertia  ellipsoid  of  the  vehicle  is  sufficiently  close 
to  a  sphere,  i.e.,  it  is  permissible  to  consider  the  centrifugal 
moments  of  inertia  and  the  differences  of  the  axial  moments  of  inertia 
to  be  of  first  order  of  smallness  in  the  axial  ones; 

-  the  characteristics  of  the  servo  engines  do  not  have  any 
delay  and  are  linear;  their  firing  does  not  disturb  the  dynamics  of 
flow  over  the  vehicle; 

the  direction  of  the  linear  air  velocity  vector  and  the 
magnitude  of  the  ram  pressure  do  not  change  in  the  process  of  a 
command  turn; 

-  the  vehicle  is  symmetric  in  weight  and  aerodynamics; 

the  hypothesis  of  stationarity  is  valid; 

specific  damping  of  the  vehicle  is  negligibly  small. 
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In  accordance  with  these  assumptions  and  the  location  of  the 
servo  units  (see  Figure  2)  the  approximate  dynamic  equations  for  the 
vehicle  become 


- c„ (a)  •  F T •  sin 9  +  Mmx  +  Al( 

VV,  “  c»  («)  [°d  (*)  -  At  J  qSe  sin  ?  +  Mynpi/  ( 5 ) 

•'t",,  *  I c„ («)  (a)  -At] cos 9  - Ct (?)  Ft)  qSe  +  Af„  ( 6 ) 

where 

I  ,1,1  are  the  moments  of  inertia  of  the  fuselage  in  the 
x  1  y  1  z  1 

XiYiZisystem; 

Cn  (a)  is  a  coefficient  of  the  normal  component  of  the 
aerodynamic  force; 

C  (a)  is  a  coefficient  of  the  tangential  component  of 
the  aerodynamic  force; 

XT,  Yt  are  the  dimensionless  coordinates  of  the  center 
of  mass; 

Cd  is  the  dimensionless  coordinate  of  the  center  of 
pressure; 

q  is  the  ram  pressure; 
l  is  the  length  of  the  vehicle; 

S  is  the  area  of  the  midsection; 

Mq  ,  Mx,Mz  are  the  components  of  the  control  moments  in  the 
x  y  XjYiZj  system  axis. 

Equations  (1-6)  describe  the  total  angular  motion  of  the 
vehicle.  An  analysis  of  this  system  of  equations  shows  that  a 
statically  stable  attitude  of  the  vehicle  is  characterized  by  a  * 
-ctbai  and  $  *  0  for  any  x.  This  means  that  the  bank  motion  of  the 
vehicle  for  a  *  ”abal  and  $  “  0  is  similar  to  the  rotation  of  a  rigid 
body  in  a  vacuum 

*  „  .^vp  .  ( 7 ) 

sin  a6/yt 
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Prom  the  condition  of  ideal  turn  (performing  obvious  transfor¬ 
mations  of  (I  -  6)  for  a  =  an<^  4>  =  0)  we  can  determine  the 

magnitude  of  the  angle  of  placement  of  jet  engines  needed  for  bank 


control 


(8; 


The  structure  of  the  bank  control  unit  in  the  simplest  case  should 
contain  the  self-equalizing  and  angular  bank  velocity  signals,  i.e., 

Mynp  =  III  (y.„  —  *)  —  kA.  ( 9 ) 

However,  a  realization  of  the  signal  x  is  difficult  in  practice. 
Instead  of  x  one  normally  uses  a  signal  from  the  PS  oriented  properly 
relative  to  the  fuselage.  In  this  case,  the  instrument  signal  is 

wnj>«a  9  cos  ayct  +  i  cos  (aJi|;T  —  aa).  O- 0 ) 


The  dynamic  properties  of  the  channels  a  and  <t>  in  accordance  with 
the  ideas  of  bank  turn  require  the  introduction  of  active  damping 


and 


M,  ■=  — 

(ID 

(12) 

It  should  be  noted  that  for  <|>  ■  0  from  (3)  it  follows  that 

u>  «  a. 


(13) 


The  information  about  the  rate  of  change  of  the  angle  4>  may  to  a 
first  approximation  be  obtained  by  differentiating  the  ratio  of 
lateral  load3 


tg?- 


(1*0 


or  the  ratio  of  the  angular  readings  of  a  pendulum  suspended  at  the 
center  of  gravity  at  the  axis  of  rotation,  the  collinear  geometric 
axis  of  symmetry  of  the  vehicle. 


It  should  be  noted  that  the  effectiveness  of  damping  in  the 
channel  (the  location  of  the  servo  units  is  shown  in  Figure  2)  does 
not  depend  on  the  value  of  the  angle  of  attack. 

A  study  of  the  controlled  angular  motion  of  aircraft  [system  of 
equations  (1  -  6)]  shows  that  the  latter  can  be  decomposed  into  "longi¬ 
tudinal,"  described  by  Equations  (3)  and  (6),  and  "lateral,"  described 
by  Equations  (1-2)  and  (^  -  5). 
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Let  us  consider  the  lateral  motion  of  aircraft.  For  this  purpose, 
let  us  represent  the  motion  in  the  form  of  a  system  of  differential 
equations  solved  for  the  higher  derivatives  x  and  £,* 

*  -  4  }  +  B  [k{  (x*  -  x)  -  hj,  [?  cos  ayct  +  i  cos  (aycr  —  aj)j},  ^ 1 5 ) 


f  +  C  {kx  (x,  -  x)  -  kt  |f  cosayct  + 


+  x  COS  (ayCT  —  a6)]}  —  j-y. 

X% 


(16) 


where 


A  »  CnWVdW-'X’tUSt  . 

y„(-*tna6  * 

giAk; 

/^•sinaj 

q  _  ayct  _  slnctyef-ctgaa  . 

'll 


<0-  =>  Cn  +  £&  c)g  tt»  (°0  — *t]  gSe 

'tx  '*» 

When  condition  (8)  is  satisfied,  the  *  motion  does  not  depend  on  the 
bank  motion  since 

C  —  ft°s  am  _  ttgOfl^ingytr  )  ^  q 

\  7^  7 


Therefore  when  afcal  is  constant,  it  is  necessary  in  practice  to 
make  sure  that  the  vector  Mq  follows  so  that  Equation  (8)  can 
be  satisfied.  Then,  as  we  know,  the  effectiveness  of  the  bank 
regulator  remains  constant. 


Otherwise  the  motion  of  the  channels  becomes  coupled  in  both 
directions.  The  character  of  the  motion  becomes  very  complicated, 
and  the  easily  visualized  physical  considerations  used  in  the  construc¬ 
tion  of  the  structure  of  regulators  and  a  selection  of  their  parameters 
become  invalid. 


In  addition,  the  peculiarities  of  the  partial  motions  of  the 
system,  including  the  well-known  effect  of  the  bank  regulator  on  the 
dynamics  of  the  $  channel  for  Mx  ■  0,  can  be  clearly  explained  by  the 


•System  (1  -  6)  is  solved  for  x  and  $  under  the  assumption  that 
a  *  -“bal*  a  =  $  is  small  (i.e.,  sin  *  “  *  and  cos*  a  1),  and 

Mv,  M_,  are  determined  from  (9),  (11),  and  (12). 

X  s  c 


34 


the^y  of  differential  equations  with  small  parameters  which  multiply 
the  derivatives  [4].  We  shall  limit  ourselves  to  the  explanation  of 
this  effect,  setting  for  simplicity  ?=z  . 


Considering  the  relatively  high  frequency  of  the  eigen  motion  with 
respect  to  ♦,  it  is  permissible  to  assume  that  u>  is  "large."  Therefore, 
;i=-=-  is  below  always  considered  a  small  parameter.  Let  us  represent 
Equation  (16)  in  the  form  of  the  system:  $  =  z, 

— ,7  ?+i*C  {—liit—ki  |Acos(ay„-a0)+2COsaycT}},  (17 ) 

where  . 

System  (17)  contains  a  small  parameter  y  which  multiplies  the 
derivative  of  z.  In  addition,  the  system  also  contains  a  small  para¬ 
meter  on  the  right-hand  side  of  the  second  equation,  where  it  enters 
in  an  irregular  fashion:  at  y  ®  0  the  right-hand  side  is ’discontinuous. 
Following  the  procedure  used  in  [3], the  substitution  of  the  variables 
)-¥=*«  is  performed,  and  then  System  (17)  becomes  !*«■=*. 

I‘r  t 

l»z  =s-«+|iC  j  (A  cos  (aycr-afi)+zcos  a>CTJ}.  ( 18 ) 

In  turn  it  is  advisable  to  write  Equation  (15)  in  normal  form.  Thus 

x  «*  A  .  , 

(“A!ix— ^  [A cos (aycy  ag)-pzcosaycT]} ,  (iq) 

Let  us  consider  the  dynamics  of  the  coupled  system  (18  -  19) >  con¬ 
sisting  of  two  "fast"  (18)  and  two  "slow)  (19)  equations.  In  this 
we  shall  confine  our  attention  to  just  the  "fast"  partial  motion, 
generated  by  System  (18),  considering  the  fact  that  it  is  coupled 
with  (19). 


Let  us  introduce  a  new  independent  variable  ("fast"  time)  x  «  t/y 


du 


and  make  the  substitution  in  (18)  and  (19).  We  get  ax'**2' 


di 


J7=  {— A,x-*jldcos(ayot— aC3,)-fz.cosayJ> 


(18') 


and 


<JA 


*2  |Acos(aycT-30M)+z.cosa>c,]} 


d  t 


(19') 
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The  report  [3]  gives  and  Justifies  a  procedure  for  an  approximate 
analysis  of  the  "fast"  phase  of  the  solution  of  systems  such  as  (18'  - 
19').  With  certain  improvements  the  "slow"  system  (19')  is  replaced 
by  another  one  which  can  be  obtained  from  (19')  by  performing  the 
following  iterations: 

1.  "Slow"  variables  x  and  A  on  the  right-hand  side  of  (19') 
are  set  equal  to  their  initial  values  (here  —  zeroes),  i.e.,  (19') 
is  replaced  by  j~=0, 

A\S-U- fljJiAj- COS  Oyer  •?.  (20) 

Therefore,  taking  into  account  the  null  initial  conditions  and 

the  first  equation  in  (18') ,  *i(*)“°  and  A,  (t)  ® Arf u (*) &\ -B\t.k2 cos aVl.T .  u (-)  . 

0 

2.  "Slow"  variables  z  and  A  on  the  right-hand  side  of  (19') 
are  set  equal  to  MT)  and  A,(t)  ,  i.e.,  System  (19')  is  replaced  by 

57 « (?)  di-Bp-kt  cos  ayct  u  (t), 


‘-fc~AnJu-Bv.  •  k2  (cos  (ayu-a0)  •  [V  J  u  (?)  d?-  ,  v 

5  (21) 

—Buftj  cos  aycT  u  (*)) ) — .  cos  ayCT  •  z  . 

In  essence  we  use  here  Picard's  method  of  successive  approximations, 
but  only  in  application  to  the  "slow"  system. 

Let  us  stop  after,  for  example,  the  second  iteration  and  neglect 
the  integral  components  in  (21),  since  they  are  proportional  to  w*. 

Then,  System  (1?')  is  replaced  by 

~Bk2  cos  oyvT  •  i** •  //. 

^^A^u+Bikl  cos  («ytT-«0u)  cosaycT.|i*«—  (22 ) 

— flA2cosaycT.(w 
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The  system  of  equations  composed  of  (18')  and  (22)  is  structurally 
equivalent  to  the  scheme  shown  in  Figure  3:  it  describes  the  dynamics 
of  an  oscillating  circuit  enclosed  by  a  series  of  feedback  couplings. 

With  null  initial  conditions,  the  integral  feedbacks  reduce  to 
a  single  feedback  with  a  resultant  amplification  coefficient.  It  is 
important  to  note  that  the  amplification  coefficients  for  all  feedback 
couplings  are  small.  In  this  case  the  effect  of  each  feedback 
coupling  on  the  dynamics  of  the  oscillating  loop  can  be  taken  into 
account  independently  of  one  another  (this  follows  from  the  theory 
of  small  parameters). 

The  effect  of  the  feedback  couplings  I  and  II  (Figure  3)  is 
obvious:  depending  on  the  sign  of  C  they  respectively  damp  the 
oscillating  link  and  increase  its  frequency  (or  conversely).  Let  us 
consider  the  effect  of  the  resultant  integral  feedback  coupling. 

The  characteristic  polynomial  (upon  disconnecting  the  loops  I  and  II) 
in  terms  of  the  time  t  becomes 

(p'-+l)p+Cf,  (23) 

where  Cf  is  the  resultant  amplification  coefficient  of  the  integral 
couplings  (small  quantity). 

For  f  ■  0  one  of  the  oscillating  roots  of  the  polynomial  (23)  P 

•  • 

pj  «  s.  For  f  t  0  the  same  root  is  pi  »  s  +  Api.  Below,  f  and  Apj 
are  assumed  to  be  small.  Substitution  of  the  new  value  of  pi  into 
(23)  and  a  retention  of  the  terms  of  the  first  order  of  smallness 
yield  Api  «  fC/2.  Thus,  the  oscillating  root  of  the  system  moves  to 
the  left  or  to  the  right  depending  on  the  sign  of  C,  but  its  imaginary 
part  remains  the  same  *).  Therefore,  from  the  point  of  view  of  the 
oscillatory  partial  motion,  the  integrating  feedback  coupling  with 
the  amplification  coefficient  "Cf"  is  equivalent  co  a  differentiating 
feedback  coupling  with  the  amplification  coefficient  "-Cf." 


*  ) 

'  With  the  aid  of  the  Rouche  theorem  [5],  one  can  prove  the 
admissibility  of  the  calculations  as  the  first  approximation  in  i. 
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The  resulting  effect  (damping  or  anti-damping)  of  the  circuit  is 
determined  by  the  mutual  relationship  of  the  "weight"  coefficients  of 
the  feedbacks  II  or  III. 

In  particular,  for  the  apparatus  shown  in  Figure  2,  the  resulting 
effect  for  ctpl  >  afaal  will  be,  as  we  know,  in  the  form  of  damping  (and 
conversely) . 


Figure  3* 

It  is  clearly  advisable  (in  order  to  retain  the  working  capacity 
of  the  system  in  case  of  a  malfunction  of  the  damper  in  the  'V  channel) 
to  make  sure  that  in  a  steady-state  motion  opl  -  ab&1  is  positive  and 
equal  to  several  degrees  (up  to  10°).  The  latter  should  be  taken  into 
account  in  the  design  of  the  tracking  system  used  to  place  Mpl  in  the 
same  direction  as  V  for  the  adaptive  channel  of  a  system  for  controlling 
the  lateral  motion. 
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CONTROL  OP  THE  PLIGHT  OP  SPACECRAFT  DURING  ATMOSPHERIC 
RE-ENTRY  BY  MEANS  OP  ON-BOARD  COMPUTERS 

Ukolov,  I.  S.j  Mltroshin,  E.  I.;  Krayzman,  V.  Ye.; 
and  Oleynlchenko,  L.  G. 

Formulation  of  the  Problem 


During  atmospheric  re-entry  at  parabolic  velocity  and  a  flight 
range  on  the  order  of  7,000  -  10,000  km,  the  flight  of  a  re-entering 
vehicle  must  be  organized  in  such  a  way  that  its  trajectory  outside 
the  atmosphere  [1]  is  as  shown  In  Figure  1.  The  altitude  of  100  km 
is  taken  as  the  upper  boundary  of  the  atmosphere. 

The  most  Important  segment,  from 
the  point  of  view  of  maximum  accele¬ 
rations  and  its  effect  on  the 
longitudinal  range,  is  that  of  the 
first  immersion.  Therefore,  below 
the  control  system  operation  will  be 
considered  only  within  that  segment, 
and  the  accuracy  of  the  system  will  be  determined  with  respect  to  the 
total  range  of  the  first  immersion  and  the  elliptical  path  outside 
the  atmosphere. 
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The  system  of  the  equations  of  motion  for  a  re-entry  vehicle  in 
the  longitudinal  plane ,  written  in  the  velocity  coordinate  system, 
has  the  following  form 

Jr-  -nX‘g-gsM 

^-VsinO 

I  a) 

~jf  cos  6 

p-foe-"".  e-e»(TrfwJ>  “x-ipv'itf 

The  basic  task  of  the  control  system  considered  in  the  present  paper 
consists  of  bringing  the  re-entry  vehicle  to  the  point  with  a  given 
value  of  the  longitudinal  range  Lu,  irrespective  of  the  large  scatter 
of  the  initial  conditions  and  the  existing  disturbances.  The  motion 
is  controlled  by  changing  the  effective  aerodynamic  lift-to-drag  ratio 
of  the  vehicle 

(2) 


where  y  is  the  bank  angle. 


One  of  the  possible  methods  of  determining  the  free  fall  portion 
of  the  trajectory  consists  of  successive  refinement  of  the  lift-to- 
drag  ratio  of  the  vehicle  by  predicting  the  longitudinal  flight  range 
using  an  on-board  computer  to  solve  System  (1).  Thus,  for  example, 
if  we  take  the  simplest  flight  program  in  which  the  lift-to-drag  ratio 
is  constant  on  the  entire  trajectory,  then  the  iteration  process  used 
to  calculate  the  required  lift-to-drag  ratio  of  the  vehicle  in  the 
(i  +  l)th  step  is  done  using  the  following  relation 
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Another  method  for  computing  the  program  of  lift-to-drag  ratio 
variation  is  proposed  in  [1]. 
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When  the  motion  is  predicted  by  integrating  System  (1),  it  is 
necessary  to  somehow  determine  the  actual  value  of  the  atmospheric 
density  p(H)  and  the  correction 


r  (t)  «  P .actual  , 
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In  addition,  due  to  the  strong  dependence  of  the  exit  (from  the 
atmosphere)  parameters  on  the  variation  of  the  density  p  (H,  L),  it 

a 

is  necessary  to  extrapolate  5(t)  and  somehow  predict  its  behavior  over 
the  entire  remaining  portion  of  the  first  immersion. 


The  limited  possibilities  of  modern  on-board  computers  require 
a  very  thorough  study  of  the  methods  of  solving  the  boundary  problem, 
and  a  thorough  analysis  of  the  equations  used  to  predict  the  flight 
with  a  view  toward  their  simplification. 

In  order  to  lessen  the  demands  on  the  storage  capacity  and  speed 
of  the  on-board  computer,  it  is  advisable  to  construct  a  control 
system  that  combines  the  principles  of  flight  prediction  with  the  feed¬ 
back  principle.  In  this  case  during  the  time  interval  needed  by  the 
computer  to  calculate  a  new  programmed  motion,  the  object  is  controlled 
not  by  an  open  system,  but  instead  by  a  feedback  principle  with  the 
advantage  that  the  motion  becomes  stabilized  even  in  the  presence  of 
various  perturbations.  This  enables  us  to  eliminate  the  unpleasant 
situation  when  the  parameters  of  motion  at  the  end  of  the  interval 
needed  to  calculate  the  control  command  differ  significantly  from  the 
calculated  values  due  to  perturbations  and  using  an  open  system.  This 
results  in  an  instability  of  the  control  process  and  in  a  sharp  change 
of  the  controlling  action  in  the  consecutive  control  intervals. 

Selection  of  Equations  Describing  the 
Process  of  Spacecraft  Re-entry 

In  order  to  select  the  mathematical  model  which  will  then  be  used 
by  the  on-board  computer  to  predict  the  re-entry  process,  it  is 
desirable  to  simplify  the  equations  describing  the  dynamics  of  space¬ 
craft  descent,  or  if  possible,  to  use  approximate  analytic  solutions 
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of  the  equations  of  motion.  The  wide  range  of  variation  of  the  initial 
conditions  and  the  required  range,  the  large  values  of  the  perturbations, 
and  the  strong  dependence  of  the  elliptical  range  on  the  parameters  of 
motion  at  the  point  of  departure  from  the  atmosphere  complicate  the 
application  of  approximate  analytic  solutions. 


Let  us  replace  the  system  of  differential  equations  (1)  with  a 
simpler  system  (5),  using  the  flight  velocity  V  as  the  independent 
variable,  replacing  the  equations  for  the  variation  of  altitude  with 
the  equation  for  the  variation  of  the  acceleration,  and  assuming 
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const  [2]. 
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The  transition  from  System  (1)  to  System  (5)  has  been  made  using  the 
following  relation 


nx 


"y8(nx  +  sinfi)  — 
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Let  us  consider  the  effect  of  various  variations  of  the  density 
p  ■  Pn  •  C  on  Systems  (1)  and  (5). 


For  5  ■  1  we  obtain  the  standard  type  of  variation  of  the  actual 
density  pa  »  p  .  Both  Systems  (1)  and  (5)  are  equivalent  In  this  case. 
For  C  ■  0.5  ■  const  or  C  ■  1.5  ■  const  Equations  (6)  imply  that  per¬ 
turbations  of  this  type  in  no  way  affect  the  accuracy  of  System  (5). 

In  general  for  any  $  ■  const  It  is  only  necessary  to  state  correctly 
the  Initial  conditions  for  System  (5),  since  when  the  initial  condi¬ 
tions  on  nx  are  given  exactly,  $  is  automatically  taken  into  account. 
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For  perturbations  of  the  type  £  *  e 
obtain 
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using  Equation  (6)  we 


(7) 


For  k  *  ±  0.5  •  10”*  [^]  the  relative  variation  of  the  density  at  the 
altitude  H  «  100  km  is 

5  *  and  for  k  -  ±  1  5  =» 


i.e,  the  density  changes  more  than  by  a  factor  of  two  in  either 
direction  from  the  nominal  value.  When  the  basic  system  (1)  is  used 
with  such  perturbations,  the  right-hand  sides  of  System  (1)  change 
by  factors  1.5  -  2.7,  respectively.  When  System  (5)  is  used  with 
analogous  perturbations,  only  one  of  the  terms  in  the  second  equation 
changes  by  3,5-7%. 


When  the  density  perturbation  changes  according  to  the  law 
assumed  in  [1] 


jl  +  e^11  +  O  sin 


(8) 


the  coincidence  of  the  solutions  of  Systems  (1)  and  (5)  without  the 
use  of  information  about  the  character  of  the  perturbations  is  some¬ 
what  disturbed.  Figure  2  shows  the  plot  of  one  of  the  coordinates 
(altitude)  corresponding  to  the  exact 
solution  of  System  (1),  assuming  that 
we  know  the  character  of  the  density 
perturbation  (solid  lines),  and  the 
corresponding  solution  of  System  (5) 
with  an  additional  third  equation 
from  System  (1)  (dotted  lines), 
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obtained  without  using  any  information  about  the  character  of 


perturbations.  In  Equation  (8)  we  set  c  ■  ±1;  f 


n 
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Thus  the  use  of  System  (5)  instead  of  System  (1)  permits  us  in 
certain  cases  to  avoid  direct  measurement  of  the  flight  altitude  and 
the  atmospheric  density,  and  use  only  the  information  about  the 
acceleration,  which  considerably  simplifies  the  necessary  calculations. 
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It  should  be  noted  that  recently  other  simplified  forms  were  proposed 
for  the  description  of  the  center  of  mass  of  the  spacecraft  [L]. 

The  Solution  of  the  Boundary  Value  Problem 

The  problem  of  bringing  a  spacecraft  to  a  given  point  is  a 
boundary  value  problem  whose  solution  involves  great  computational 
difficulties.  In  addition ,  in  a  number  of  cases  various  additional 
requirements  are  Imposed  on  the  seleoted  control  functions  which 
complicates  the  solution  of  such  problems  even  more. 

As  we  know,  when  solving  optimisation  problems  in  accordance  with 
the  classical  methods  of  variational  calculus  or  the  method  based  on 
the  Pontryagln  maximum  principle*  it  becomes  necessary  to  solve  a 
system  of  differential  equations  whose  order  is  twice  that  of  the 
system  describing  the  physical  process*  and  Includes  two  boundary 
conditions. 

When  solving  such  problems  numerically*  one  needs  to  supplement 
the  Incomplete  boundary  conditions  at  the  initial  point*  then  to 
numerically  integrate  the  differential  equations  up  to  the  final 
point*  to  see  to  what  extent  the  boundary  conditions  at  the  final 
point  thus  obtained  correspond  to  the  given  conditions*  and  finally 
to  Improve  the  inaccurate  Initial  conditions.  This  process  must  be 
repeated  until  all  final  conditions  are  satisfied. 

These  methods  are  of  little  use  If  one  wants  to  obtain  numerical 
solutions  in  real  time  using  an  on-board  computer. 

Since  the  problem  must  be  solved  within  a  short  period  of  time* 
it  is  desirable  to  use  the  most  rapidly  converging  procedures.  For 
this  purpose*  we  use  a  method  which  Is  analogous  to  the  procedure 
of  the  X-matrlx  control  [5,  6].  This  permits  us  to  Improve  the  con¬ 
vergence  of  the  iterative  process  and  thus  reduce  the  number  of 
iterations*  since  when  integrating  the  system  in  the  forward  direction 
one  uses  the  information  about  the  current  deviations  of  the  phase 
coordinates  of  the  actual  motion  of  the  object  from  the  intermediate 
reference  trajectory. 
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If  System  (5)  is  written  in  a  general  form 


(9) 


where  x  is  the  state  vector  of  the  system,  then  the  equations  of 
motion  in  variations  with  respect  to  the  reference  trajectory  have  the 
form 


>B(V)bx  +  e(V)-U<. 


Adding  the  conjugate  system  of  equations  to  System  (10) 

<=<  (11) 

one  can  obtain  the  fundamental  Bliss  formula  [3,  4],  and  take  into 
account  the  change  in  the  final  point 

.  * 

(>/•  -  S  V-c)(U<)dV  +  Q.'*x)v  (12 ) 

h 

(the  prime  symbolizes  the  operation  of  transposition).  We  impose  on 
the  scalar  product  (X'.Ajc)  the  condition 

W-LX)Vltw-Ua.  (13) 

Taking  this  condition  into  account,  the  initial  conditions  for  System 
(11)  are  determined  in  the  following  way 


where  x*  is  che  phase  coordinate  defining  the  endpoint  of  the  control 

dL,\i  dx* 

interval,  ~ai  510(1  37  are  partial  derivatives  with  respect  to  the 
corresponding  coordinates,  and  jp  are  total  derivatives  determined 
by  the  corresponding  equations  of  System  (5). 


In  connection  with  the  fact  that  the  lift-to-drag  ratio  at  our 
disposal  is  small,  it  appears  advisable  to  choose  the  control  functions 
such  that  the  required  change  in  the  range  ALj 3  will  be  achieved  by 
minimizing  the  change  in  the  control  program  (in  the  mean-square  sense). 

Thus  we  arrive  at  a  variational  problem  involving  a  conditional 
extremum.  We  are  required  to  determine  the  law  of  variation  of  the 
control  action  (AK)  such  that  the  following  expression  will  be  minimized 

Vk 

\{M<y-dV  (15) 
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subject  to  the  constraints  (12)  and  (13).  This  type  of  problem  involv¬ 
ing  a  conditional  extremum  subject  to  constraints  in  an  integral  form 
is  called  isoperimetric . 


An  isoperimetric  problem  can  be  reduced  to  the  Lagrange  problem, 
i.e.,  to  the  minimization  of 

\(Afl)W 

(15) 

subject  to  the  condition  ^(A IU)J,-(X'-<?)A/C  •  In  order  to  minimize  the 
expression  in  (15),  we  construct  the  intermediate  function 

H  -  (A/C)’  +  p  (V)Wy  ~  M  ( 16 ) 

and  solve  the  corresponding  Euler  equation 

-0  (17) 
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Prom  Equations  (17),  we  obtain 

2(A/0-,*(lW'*«)-0  (18) 

o-a^M-o 

The  second  equation  in  (18)  implies  that  p  does  not  depend  on  V,  i.e., 
p  ■  const,  and  from  the  first  we  obtain 

A*— 5^  M -*(*'•*)  (19) 

where  pi  is  a  constant  to  be  determined. 


Using  (19),  (12),  and  (13),  we  obtain 

*>  (20) 
dLu-y&'cYdV+W'lx)?, 

whence 

dLn-(V-6x)y, 

* — r> - 

j  (V-e)W  '«*/ 

The  obtained  expression  (19)  for  the  variation  of  the  control  (AK) 
has  a  definite  physical  meaning.  The  variation  of  the  functional  (20), 
when  (AK)  is  selected  in  accordance  with  relation  (19) ,  always  has  the 
same  sign  as  the  sign  of  the  variation  of  control  AK. 

Equation  (19)  gives  the  direction  of  the  gradient,  i.e.,  the 
direction  of  the  maximum  effect  of  (AK)  on  ALn. 
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V/e  introduce  a  function  A(V)  by  means  of  the  equation 

(IV 


(22) 


for  which  the  boundary  condition  is  L  =0. 

vk 

Using  (19) ,  (21),  and  (22),  we  finally  obtain 

AK  (V)  -  -  A, (10  [AI»  +  (X'-Ajc)„J 


(23) 


where 


Equation  (23)  determines  the  control  law  for  a  closed  system,  which 
takes  into  account  the  effect  of  the  input  ALl3  and  the  deviations  of 
the  actual  motion  from  the  reference  trajectory. 


The  difference  between  this  method  of  solving  the  boundary- value 
problem  and  the  methods  discussed  in  [1],  is  that  in  solving  the 
boundary-value  problem  we  use  the  information  about  the  deviation  of 
the  actual  motion  from  the  nominal  motion,  and  in  forming  the  control 
AK  we  only  use  a  single  integration  of  the  fundamental  nonlinear 
system  (5).  Instead  of  a  second  integration  of  the  nonlinear  system 
describing  the  dynamics  of  a  descending  vehicle,  we  integrate  the 
associated  linear  system  (11)  and  Equation  (22)  in  order  to  determine 
new  control  coefficients. 

A  block  diagram  of  the  control  system  realizing  this  control 
algorithm  is  shown  in  Figure  3. 


Figure  3:  1  -  object;  2  -  data  sensors; 

3  -  information  conversion  block;  4  - 
prediction  block;  5  -  automatic  pilot. 
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The  information  about  the  parameters  of  motion,  obtained  by  means 
of  the  data  sensors,  is  fed  through  the  data  conversion  block  to  the 
parameters  of  motion  prediction  block  of  the  spacecraft  and  to  the 
automatic  pilot  which  determine  the  command  value  of  the  effective 
aerodynamic  lift-to-drag  ratio. 

Conclusion 

The  system  discussed  here  which  involves  an  on-board  computer 
connected  into  the  control  loop  can  be  applied  to  various  objects 
subject  to  large  perturbing  forces  and  a  large  area  of  scatter  of  the 
initial  conditions,  where  the  presently  known  systems  using  a  fixed 
nominal  trajectory  do  not  meet  our  requirements.  In  addition,  such 
control  systems  may  be  uped  in  such  Important  cases  as  when,  during 
the  operation  of  the  system,  it  is  necessary  to  change  the  trajectory 
of  motion  in  a  direction  unknown  beforehand  —  for  example,  when 
guiding  a  spacecraft  toward  a  moving  point.  Similar  control  systems 
can  also  be  applied  in  spacecraft  returning  from  other  planets.  A 
characteristic  feature  of  interplanetary  flights  is  the  fact  that 
spacecraft  re-enter  the  atmosphere  at  a  very  large  speed  which  is 
much  larger  than  the  parabolic  speed.  This  places  high  requirements 
on  the  accuracy  of  the  necessary  spacecraft  maneuvers  in  the  atmosphere. 

At  superparabolic  speeds,  the  operation  of  the  control  system 
subject  to  large  density  perturbations  is  greatly  improved  if  the 
control  system  include^  ait  on-board  digital  computer  that  can  calculate 
in  real  time  the  maneuver  required  when  perturbations  are  present. 
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PARAMETRICALLY  INVARIANT  AUTOMATIC  CONTROL  SYSTEMS 
WITH  A  LINEAR  PHYSICALLY  REALIZABLE  REGULATOR 

Sokolov,  N.  I.;  Makovlev,  V.  I.; 
and  Lipatov,  A.  V. 

When  synthesizing  adaptive  control  systems,  we  are  usually  faced 
with  the  problem  of  maintaining  a  system's  performance  indicators 
constant  or  changing  them  in  a  desired  direction  with  a  change  in  the 
parameters  of  the  object. 

There  are  papers  available  in  which  the  solution  of  the  problem 
is  attempted  by  using  systems  of  variable  structure,  self-adjusting 
systems  of  various  types,  and  systems  with  an  infinitely  large  ampli¬ 
fication  coefficient  (for  example,  1,  2,  3»)* 

In  this  paper  we  show  that  it  is  possible  to  construct  linear 
physically  realizable  regulators  with  a  rigid  structure  and  constant 
parameters  in  order  to  control  objects  whose  parameters  change  within 
a  wide  range.  It  is  shown  that  such  regulators  are  capable  of  keeping 
the  performance  indicators  of  a  dynamic  system  constant  within  any 
finite  range  of  the  parameters  of  the  controlled  object. 

It  is  assumed  that  the  object  is  controlled  by  means  of  some 
linear  actuator.  Suppose  that  the  object  is  described  by  a  differential 

4.  h 

equation  of  m  order 
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x(m)  +  S  «,(/)  *(,)  -  «o(OKo6(/)  r.  ( 1 ) 

lm  0 

and  that  the  actuator  is  described  by  a  differential  equation  of  the 
(n  -  m)th  order 

+  Tv-M.  ... 

/«o  ( 2 ) 

where  x  is  the  controlled  coordinate; 

$  is  the  actuator  output; 
u  is  the  control; 

a^t)  are  the  variable  parameters  of  the  object,  varying  within 
finite  limits; 

K  „(t)  is  the  variable  transfer  coefficient  of  the  object 
ob 

bi  are  constant  parameters  of  the  actuator. 

The  problem  is  to  control  the  object  in  such  a  way  that  the 
performance  indicators  of  the  system  will  remain  constant  while  the 
parameters  of  the  object  are  changing. 

The  control  must  be  effected  using  the  input  coordinate  of  the 
system  and  the  signals  generated  by  physically  realizable  correctors 
that  transform  the  output  coordinate  of  the  system  and  a  number  of 
measurable  intermediate  coordinates  of  the  system. 

In  the  case  of  complete  observability,  i.e.,  when  it  is  possible 
to  use  n  -  1  derivatives  of  the  output  coordinate  of  the  system,  one 
can  obtain  the  required  dynamic  properties  within  any  range  of  the 
object  parameters  [1,  3»  ^3. 

Under  the  actual  conditions,  it  is  not  possible  to  obtain  and 
use  "pure"  derivatives,  particularly  those  of  high  order.  Therefore, 
it  is  of  interest  to  consider  the  possibilities  inherent  in  linear 
systems  with  physically  realizable  correctors.  Physically  realizable 
correctors  will  be  defined  as  correctors  for  which  the  order  of  the 
numerator  of  the  transfer  function  does  not  exceed  the  order  of  the 
denominator. 
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Figure  1. 


Figure  1  is  a  diagram  of  the  control  system  for  which  the 
controlled  object  and  the  actuator  are  described  by  Equations  (1)  and 
(2),  and  the  physically  realizable  correctors  C0,  Ci,  C2,  Cs  are 
described  by  the  differential  equations 
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The  corrector  C»,  described  by  Equation  (3),  is  introduced  to 
eliminate  the  derivatives  on  the  right-hand  side  of  the  differential 
equation  for  the  closed  system  [5].  The  motion  of  the  closed  system 
is  then  described  by  the  equation 

P»+AP)  &  “  Co/CKoSOo^k^oi  ( 7 ) 

where  Pn+r,(p)  is  a  differential  operator  of  the  form 
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Depending  on  the  degree  j  of  the  operator  p,  the  form  in  which 
the  selected  parameters  ct  +  K^i,  Aj,  dt  enter  into  the  coefficients  Aj 
of  the  characteristic  polynomial  (8)  changes.  Three  forms  of  the 
coefficients  Aj  of  the  characteristic  polynomial  are  possible 


a)  m  +  r+  1  </<«  +  r 

mini/,  a)  min  If,  m] 

Aj~  2j  i  (Cl-i  4*  /C[*»-j) 

1  f«mu  (U.  )-t\  /-out  [0, 1-n+m]  1  ’ 

b)  r+  1  </<m  +  r 


min  I/,  a] 


min  If,  ai) 


litl/,  a]  min  (I, 

/-mi*  10,  l-r]  /— mi*  to.  /-a+m| 


min  I/,  m) ' 

s  W/-I 

/-mi«(U,  l-r | 


c)  0</<r 

minimal  mln(/.  m) 

At  «*  i  2  +  KPj-i)  A- 

'  /-mi>|0.  l-r]  I-idu [0,  /-a+m] 

min  |{,  m] 


+  2  4*  KKa(AtP^ 

/-ou*fo,  j-r\ 


I- 


(9) 


(10) 


(11) 


Let  us  find  the  relationship  between  the  smallest  necessary  order 
of  the  correctors  and  the  orders  of  the  equations  describing  the 
object  and  the  actuator.  In  order  to  uniquely  construct  the  coeffici¬ 
ents  Aj  of  the  characteristic  polynomial  (8),  it  is  necessary  that 
each  coefficient  of  the  characteristic  polynomial,  starting  with  the 
coefficient  An+  ,  contain  successively  a  new  parameter  out  of  the 
set  of  the  selected  parameters  Cj,  Aj,  dj.  The  number  of  the  arbi¬ 
trarily  selected  parameters  c^  is  equal  to  r.  The  number  of  the 
arbitrarily  selected  parameters  Aj  is  equal  to  (r+1),  and  the  number 
of  the  parameters  d^  is  equal  to  (r+1).  Since  the  parameter  Aj  enters 

the  coefficient  .A  ...  first,  the  value  of  r  must  be  chosen  such  that 
m+r 

the  parameters  will  form  the  coefficients  of  the  characteristic 
polynomial  from  An+  to  Am+r+i  inclusive,  i.e., 

or  r-(a  +  r— l)-(i»  +  r+ 1)+1  (12) 


r(  *■«  —  /«—  1. 


52 


Since  the  parameter  d^  enters  the  coefficient  Ar  first,  the  r 
parameters  c^  together  with  the  (r+1)  parameters  must  form  the 
coefficients  of  the  characteristic  polynomial  from  An+r__j  to  Ar+j 
inclusive,  i.e., 

r  +  (r+l)«(»  +  r-l)-(r+l)+l 
or 

,  n—2  (13) 

2  • 


We  obtain  two  possible  values  of  r:  ri  and  r2.  The  larger  is 
taken  as  the  required  value  of  r.  If  the  number  is  a  fraction,  which 
occurs  when  n  is  odd,  then  the  value  of  r  is  rounded  off  to  the 
nearest  higher  integer 


>  max 


xfft  —  m  —  1, 


(1U) 


We  shall  give  an  example  of  how  to  select  the  order  of  the 
correctors.  Suppose  that  the  order  of  the  equation  of  the  object  is 
m  ■  2,  and  that  the  actuator  has  the  order  n  -  m  *  3»  (n  *  5). 


Then  from  (12) 
and  from  (13) 
and  r  —  max  (r„  rt\  —  2. 


f|  «a  3  —  1  k  2, 


1.5 


If  the  order  r  of  the  corrector  is  selected  properly,  then  it 
becomes  possible  to  obtain  any  deBired  performance  for  any  possible 
mode  of  operation  of  the  system  described  by  Equation  (7).  However, 
due  to  the  variation  of  the  parameters  of  the  object  —  in  particular 
the  component  KQbao  in  the  case  of  other  operating  modes  —  it  is 
possible  for  the  dynamic  properties  of  the  system  to  deteriorate  all 
the  way  to  a  loss  of  stability. 


In  order  to  obtain  a  system  which  Is  operational  in  the  entire 
range  of  the  admissible  operating  conditions,  we  suggest  the  following 
approach  to  the  formulation  of  the  variable  characteristic  polynomial 
(8).  Each  of  the  coefficients  Aj  of  this  polynomial  contains  both 
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constant  terms  (due  to  the  fact  that  a^  ■  1),  and  variable  terms. 

By  properly  selecting  the  value  of  K,  we  can  achieve  a  situation  in 
which  the  constant  terms  will  be  larger  by  an  order  of  magnitude  than 
the  sum  of  the  variable  terms  in  each  coefficient  Aj  for  j  >  r.  Since 
the  terms  containing  the  component  KQbao  usually  vary  the  most,  it  is 
advisable  to  write  the  characteristic  polynomial  (8)  in  an  approximate 
form  reflecting  the  variability  of  this  component. 

FnJrr  (p)  -  + . . .  +  /W,+l  + 

+  W+47)p'  +  ...  +  (^+^)  (15) 


where 


Tk 


Koba» 

oba#)min 


and  changes  from  1  to  ^Koba®^niax 

Yraax  “  <Koba#>min 

A',  A'' are  constants. 

Prom  the  structure  of  the  characteristic  polynomial  (15),  we 
can  see  that  a  change  in  y  may  result  in  a  displacement  of  its  roots 
to  the  right  half-plane  of  the  root  plane  even  assuming  that  the 
system  has  a  sufficient  stability  margin  for  any  y.  The  effect  of 
y  on  the  roots  of  the  characteristic  polynomial  can  be  controlled  to 
suit  our  purposes  by  writing  the  polynomial  (15)  in  the  form  of  a 
product  of  two  polynomials:  one  of  order  n  -  1  and  the  other  of 
order  r  +  1: 


Fn+r (p) -  F„.x (p) Fr+X (p)  ~(/>'rl  +  B„.ip'ri  +  •  •  •  +  B0) X  (i5) 

X(^+'  +  Cf/r+  ...  +CU) 


If  the  natural  frequency  of  the  polynomial  F 


than  the  natural  frequency  of  the  polynomial  F 
one  order,  i.e., 

»Tc„. 


r+i 


n_  (p)  is  greater 
(p)  by  approximately 

(17) 


then  the  polynomials  F  (p)  and  pr+l(p)  m®y  be  approximately 
expressed  in  terms  of  the  coefficients  of  the  polynomial  (15)  in 
the  following  way  [6]: 
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Fn-i  ip)**  P*~l  +  An+f.ip"'2  +  •  •  •  +  A,+2p  +  A,+ ,; 


(18) 


Ff+i  (P)  ~  [Ar+lp'+'  +  (A;  +  A)i)p'  +  •  •  •  +  (A'0  +  7^1-  (!9 ) 

If  the  condition  (17)  is  satisfied  for  all  y  (and  this  can  always  be 
achieved  by  properly  choosing  cA  +  Kyl,  A^  dA  for  any  finite  range 
of  the  object  parameters  and  any  finite  speed),  then  the  performance 
indicators  of  the  system  are  completely  determined  by  the  polynomial 

(19) ,  and  a^  can  be  selected  such  that  the  inequality  A ' <*  A"  is 
satisfied,  and  then  the  polynomial  (19)  becomes 

Fr+i  (P) »  ^  l'WiP'+1  +AW  +  •  •  •  +  Al  7l  (20) 

One  can  always  achieve  a  situation  in  which  for  y  *  1  one  of 
the  roots  of  the  polynomial  (20)  is  so  far  from  the  imaginary  axis 
that  it  will  have  little  influence  on  the  dynamic  properties  of  the 
system.  According  to  [6],  the  value  of  this  root  o  is  approximately 
and  as  y  increases  the  root  will  affect  the  dynamic  proper- 

Ap. ft 

ties  of  the  system  even  less.  However,  the  remaining  roots  of  the 
polynomial  (20),  determining  the  performance  indicators  of  the  system, 
can  be  obtained  from  the  polynomial  Pr(p)  with  constant  coefficients 

P,  (P)  *  w  +  -+4I  (21) 

If  the  absolute  value  of  the  increasing  root  a  of  the  polynomial 

(20)  approaches  the  average  absolute  value  of  the  roots  of  the  poly¬ 
nomial  Pjj.^P)*  then  the  group  of  roots  with  the  largest  modulus  will 
be  determined  by  a  polynomial  of  degree  n  with  a  free  term  dependent 
on  y.  As  y  increases,  the  system  may  become  unstable  (for  n  >  2)  as 
a  result  of  an  increase  of  the  free  term  in  the  polynomial. 


Thus,  on  one  hand,  the  variable  root  a  must  be  greater  than  the 
roots  of  the  polynomial  Fr(p)  which  mainly  determines  the  performance 
indicators  of  the  system,  and  on  the  other  hand,  it  must  not  approach 
the  roots  of  the  polynomial  F^^p),  !•«.,  thus  for  all  y  we  must 
have  the  inequalities 


(22) 
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The  inequalities  in  (22)  are  a  sufficient  condition  (but  not 
necessary!)  for  the  stability  and  constancy  of  the  performance 
indicators  of  the  system  for  any  finite  range  of  the  parameters  of 
the  controlled  object. 


In  practice,  a  similar  result  can  be  obtained  with  less 
restricting  inequalities.  It  is  not  necessary  to  require  that  for 
any  y  the  right-hand  inequality  in  (22)  be  satisfied.  We  may  assume 
that  for  large  values  of  y  the  root  a  reach  the  region  of  the  roots 
with  the  largest  modulus,  corresponding  to  the  polynomial  (18).  As 
Y  Increases  further  a  pair  of  complex  conjugate  roots  appears  of  the 
polynomial  of  degree  n  with  a  free  term  dependent  on  y,  and  the  pair 
moves  toward  the  imaginary  axis  of  the  root  plane.  The  absolute 
value  of  this  pair  of  roots  remains  quite  large,  and  as  a  result 
these  roots  have  hardly  any  effect  on  the  performance  indicators  of 
the  system,  even  if  they  move  all  the  way  to  the  Imaginary  axis  in 
;he  root  plane.  The  character  of  motion  of  the  roots  as  y  changes 
is  shown  in  Figures  2  and  3. 
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Figure  2. 


The  polynomial  pr+1(P)  (20)  can  be  constructed  in  such  a  way 
that  —  as  y  increases  —  not  Just  one  but  two,  three,  or  more  roots 
become  larger.  The  condition  for  the  separation  of  the  absolute 
values  of  the  roots  (IT)  of  the  polynomials  (pn_j(P)  and  Fr+j(P)  then 
changes.  For  example,  in  the  case  when  two  roots  of  the  polynomial 
Fr+1(p)  increase  with  an  increase  in  y*  we  write 


(23) 
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which  is  a  less  restricting  condition  than  Condition  (22). 

Finally,  if  the  polynomial  Fr+l(P)  is  constructed  in  such  a 
way  that  as  y  Increases  sc  do  all  of  its  (r  +  1)  roots,  then  the 
performance  indicators  of  the  system  do  not  remain  constant,  but 
instead  vary  within  a  certain  range. 

Thus  linear  systems  with  a  finite  amplification  coefficient  with 
physically  realizable  correctors  can  in  principle  satisfy  the 
requirements  on  the  dynamic  properties  of  a  system  for  a  finite  range 
of  the  parameters  of  the  controlled  object.  Moreover,  the  structure 
and  the  parameters  of  the  regulator  remain  constant. 
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APPLICATION  OP  "APPROXIMATE"  STABILITY  CRITERIA 
TO  THE  SYNTHESIS  OP  ADAPTIVE  SYSTEMS 

Sokolov ,  N.  I.  and  Lipatov,  A.  V. 

When  synthesizing  adaptive  systems,  we  usually  deal  with  a  set 
of  systems  with  "frozen"  coefficients.  Usually  it  is  required  that 
the  entire  set  of  systems  as  a  whole  or  some  of  its  subsets  satisfy 
certain  performance  criteria.  In  the  analysis  and  synthesis  of 
systems,  one  has  to  know  the  relationships  between  the  performance 
on  one  hand,  and  the  variable  parameters  of  the  object  or  parameters 
of  the  regulator  (constant  or  variable),  on  the  other.  In  particular. 
It  may  become  necessary  to  establish  such  relations  between  the 
varying  parameters  of  the  object  and  the  varied  (or  constant)  para¬ 
meters  of  the  regulator  such  that  the  performance  remains  constant 
or  changes  according  to  a  certain  law. 

In  this  paper,  we  attempt  to  supply  the  means  for  establishing 
this  type  of  relationship  for  one  of  the  partial  performance 
criteria  of  a  system  —  namely,  stability.  The  application  of 
existing  stability  criteria  to  establish  such  relationships  is 
possible,  but  in  practice  this  is  impossible  for  systems  of  high 
order . 

There  is  a  need  for  approximate  but  simpler  stability  conditions 
which  could  eliminate  this  drawback.  The  approximate  stability 
conditions  will  be  defined  as  either  only  necessary  conditions  for 


58 


stability  or  only  sufficient  conditions.  Geometrically  this 
corresponds  to  replacement  of  the  stability  region  in  a  k-dimensional 
space  of  "k"  parameters  that  are  of  interest  to  us,  and  for  which  the 
equation  of  its  boundary  is  too  complicated,  with  another  region 
inscribed  into  the  exact  stability  region,  or  described  around  it, 
for  which  the  equation  of  the  boundary  has  a  simpler  form. 

One  of  such  "approximate"  simple  stability  conditions  involves 
the  necessary  condition  requiring  that  all  the  coefficients  of  the 
characteristic  polynomial  of  a  system  be  positive.  This  condition 
also  specifies  a  certain  region  in  parameter  space,  but  it  differs 
too  much  from  the  exact  stibility  region  for  systems  of  order  higher 
than  two.  This  makes  it  impossible  to  use  this  condition  by  itself. 

Let  us  formulate  two  theorems  about  the  necessary  conditions 
of  stability. 

Theorem  1.  In  order  that  the  polynomial 

F(*)- +  +  . ..  +  an-,  «,>0  (!) 

have  ail  roots  lying  in  the  left  half-plane,  it  is  necessary  that  the 
following  inequalities  be  satisfied 

>  1;  <  — 1,2 . n- 2.  (2) 


Theorem  2.  In  order  that  the  polynomial  (1)  have  all  roots 
lying  in  the  left  half-plane,  it  is  necessary  that  the  following 
inequalities  be  satisfied 

i-2,3 . n-2,  (3) 

where 


Ct 


(— 1)  «+!  +  Z 


h. 


l-l)‘+3 


) 


(._.l+i=Jl£=L)(<+  S=!jl=±y 


(4) 
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The  values  of  Ci  for  each  n  and  i  can  be  calculated  and  listed 
in  a  table. 


Both  conditions  involve  a  certain  parameter  Xi>  defined  as 

.  _  a‘°‘"  (  k  ) 


In  terms  of  this  parameter  conditions  (2)  and  (3)  can  be 
rewritten  as: 


)./>!;  {-1,2,...,  «-2 


(6) 


and 


M’l+i  ^  C/j  { —  2, 3* ,  n  —  2. 


(7) 


Each  characteristic  polynomial  of  order  "n"  is  characterized  oy 
n  -  2  parameters  X.  The  number  of  the  parameters,  X±,  is  smaller  by 
one  than  the  number  of  the  coefficients  determining  the  characteris¬ 
tic  polynomial  in  a  normalized  form.  Nevertheless,  it  can  be  shown 
that  the  stability  of  the  system  describing  this  polynomial  depends 
only  on  the  parameters  X^ 


Suppose  that  we  have  a  polynomial  of  degree  n  with  a  given  system 
of  parameters  X^,  Since  for  a  complete  description  of  the  polynomial 
in  this  case  we  need  two  or  more  parameters,  they  will  be  assigned 
arbitrary  values.  Let  us  write  the  polynomial  in  the  form 


F,(z)  -  *»  +  +  Bz  -*  +  + 


4.j4£L*.-»  + _ 

T  X;'X,X, 


./»— 8 


+ ... 


(8) 
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Here  A  and  B  are  arbitrary  parameters  which  define  a  set  of 
polynomials  having  the  same  system  of  parameters  The  only  con¬ 
ditions  on  A  and  B  are  that  they  must  be  greater  than  zero.  We  know 
that,  if  the  initial  polynomial  has  all  roots  in  the  left  half-plane, 
then  they  remain  there  upon  the  substitution  z  *  zlk  (by  the  scaling 
theorem),  and  when  all  even  or  odd  coefficients  are  multiplied  (or 
divided)  by  the  same  positive  number,  the  roots  also  remain  in  the 
left  half-plane  (by  Mikhaylov's  stability  test).  Therefore,  without 
changing  the  character  of  stability  of  the  polynomial  in  (8),  we 
make  the  following  transformations.  We  make  the  substitution  z  - 
Then  we  divide  the  even  coefficients  of  the  polynomial  thus  obtained 
by  (/B)n,  and  the  odd  ones  by  A(/B)n“1.  We  obtain  the  polynomial 
(dropping  the  primes ) : 


F„  (2)  =•  2"  +  2«->  +  gn-1  +  jl  2"-3  +  jJ-_  2«- 1  + 

+  -r — J"-5  4 - i—  2"“'5  + 

>;;u,  ww, 


(9) 


Thus,  it  is  clear  that  any  polynomial  in  the  set  of  possible 
polynomials  having  the  same  system  of  parameters  X.^  can  be  reduced 
to  a  polynomial  whose  coefficients  depend  on  X^  alone  by  means  of 
transformations  that  do  not  affect  its  stability.  Since  the  stability 
of  a  system  described  by  this  polynomial  depends  only  on  the  para¬ 
meters  X^,  the  stability  of  a  system  describing  the  original  poly¬ 
nomial  also  depends  on  the  parameters  X^ 

The  system  of  parameters  Xi  completely  describes  the  stability 
of  the  system  but  hardly  characterizes  its  quality.  For  example,  the 
two  characteristic  polynomials: 

/?5l(2)=  2,  +  52,  +  102)+102J  +  52+l 

Fc.(z)  =»  2s +502'+ 10z5+  100^+52  + 10  1 

are  described  by  the  same  system  of  parameters  X^  and  are  factored 
as  follows: 

Fm  (*)  =*  (*  +  !)3 

F «  (*)  -  (*'  +  ‘19.84)  (2-  +  0,0472  +  0,105)  (r  +  0,112  +  1,89)  ( 11 ) 
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The  character  of  roots  from  the  point  of  view  of  oscillations 
is  different,  but  from  the  point  of  view  of  stability,  these  two 
polynomials  are  equivalent.  The  systems  described  by  these  poly¬ 
nomials  will  also  have  the  same  stability  margin  in  any  of  the  para¬ 
meters  entering  the  coefficients  of  the  characteristic  polynomial  in 
identical  fashion.  We  shall  assume  that  any  parameter  enters  the 
coefficients  of  polynomials  in  identical  fashion  if  it  produces  iden¬ 
tical  fractions  in  identical  coefficients.  Thus,  the  free  term  in 
both  equations  can  be  increased  until  instability  occurs  the  same 
number  of  times.  This  example  shows  that  it  is  incorrect  to  judge 
the  stability  margin  of  a  system  on  the  basis  of  the  oscillatory 
character  of  the  roots  of  its  characteristic  polynomial. 

Since  the  system  of  parameters  A^^  determines  completely  the 
stability  of  a  system,  it  is  natural  to  construct  the  regions  of 
stability  of  the  system  in  the  (n  -  2)-dimersional  space  of  the 
parameters.  Then  the  necessary  conditions  for  stability  (2)  and 
i 3)  will  also  yield  a  certain  region  which  comprises  the  exact  region 
of  stability  and  may  serve  as  an  "approximate"  region  of  stability. 

Prom  the  point  of  view  of  the  regulator  synthesis,  it  would  be 
more  useful  to  have  ;ome  simple  sufficient  conditions  for  stability 
that  would  guarantee  the  stability  of  a  system.  It  turns  out  that 
such  conditions  may  be  obtained  also  by  using  the  parameters  A^ 

We  can  prove  the  following  theorem  about  the  sufficient  condi¬ 
tions  for  the  stability  of  a  system  having  the  characteristic 
polynomial  (1). 

Theorem  3.  If  all  roots  of  all  possible  polynomials  of  fifth 
degree,  formed  out  of  six  consecutive  coefficients  of  the  polynomial 
(1),  lie  in  the  left  h^lf-plane,  then  all  the  roots  of  the  polynomial 
(1)  also  lie  in  the  left  half-plane. 
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Thus,  one  can  obtain  a  stable  system  of  nfch  order  by  achieving 
the  stability  of  n  -  4  systems  of  fifth  degree,  which  is  much  simpler 
to  do.  The  stability  conditions  for  a  system  of  fifth  degree, 
expressed  in  terms  of  the  parameters  X  ,  have  the  form 


x,  >  l 

X1X2X3  +  2X1X3  +  XjX^Xj  —  XjXgX] — X1X3X3 — xfxl  —  1>0 


(12) 


As  we  can  see  from  Equation  '12),  these  conditions  are  also  quite 
complicated  which  will  hinder  their  use  in  practice  when  investigating 

1.U 

the  stability  of  systems  of  n  degree.  But  this  expression  can  be 

considerably  simplified.  We  can  use  the  ''approximate''  sufficient 

conditions  instead  of  the  exact,  necessary  and  sufficient,  conditions 

for  the  stability  of  a  system  of  fifth  degree.  By  investigating  the 

second  inequality  in  (12)  graphically  or  analytically,  one  can  3how 

that  it  is  valid  if  all  A.  >  2.144  ...  (i  -  1,  2,  3)  or  if  all 

a  +  <  0.89  ...  (i  ■  1,  2).  Using  this  and  Theorem  3,  we  can 

Ai  Ai+i 

formulate  the  following  two  theorems  about  the  sufficient  conditions 
for  the  stability  of  linear  stationary  systems. 

Theorem  4.  If  the  polynomial  in  (1)  satisfies 

X, > 2,144  . ,  i-  1,  2,  ....  »-2,  (13) 

then  all  its  roots  lie  in  the  left  half-plane. 

Theorem  5.  If  the  polynomial  in  (1)  satisfies 

1  +  J-  <  0,89  .. i=l,2 . n  — 3,  (14) 

Xj  Xj+i 


then  all  its  roots  lie  in  the  left  half-plane. 

Conditions  (13)  and  (14)  determine  a  certain  region  in  the  space 
of  the  parameters  A^  which  is  inscribed  into  the  exact  stability 
region,  determined  U3ing  the  exact  necessary  and  sufficient  conditions 
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for  stability.  This  region  may  also  be  used  as  an  approximate  region 
of  stability. 

Let  us  consider  several  examples  showing  how  "approximate" 
conditions  for  stability  are  used. 

Example  1.  Suppose  we  are  given  the  characteristic  equation  of 
a  system 


z*  -j- 122=  +  47z<  +  1082s  +  1222s  +  Kit  +  Kx  -  0  (15 ) 


We  are  required  to  find  the  "approximate"  regions  of  stability 
in  the  space  of  the  parameters  Kj  and  K2. 

In  this  case,  we  have  selected  a  simple  polynomial,  so  that  for 
comparison  we  can  obtain  the  "exact"  boundary  of  stability  using  the 
method  of  D-decomposition. 

The  necessary  conditions  for  stability  of  (2)  and  (3)  and  the 
requirement  that  the  coefficients  of  the  polynomial  be  positive  impose 
on  the  parameters  Ki  and  K2  the  following  conditions: 

0  <  /r,  <  105 

ATS  <242  (16) 

K,  >  0,89/e, 


If  even  one  of  these  conditions  is  not  satisfied,  the  system 
becomes  unstable.  On  the  other  hand,  if  all  the  conditions  are 
satisfied,  this  still  does  not  guarantee  stability. 

The  sufficient  conditions  for  stability  (13)  and  the  requirement 
that  the  coefficients  be  positive  impose  on  the  parameters  Kj  and  K2 
the  following  conditions: 


0  <Ki  <0,52  Kt 
Kj  <  100 


(17) 


6H 


If  these  two  conditions  are  satisfied ,  the  stability  of  a  system 
is  guaranteed. 

And,  finally,  the  sufficient  conditions  in  (lb)  yield 

0<Kl<K*-0,mKl  (18) 

Kt  <  109 

If  both  of  these  conditions  are  satisfied,  independently  of  the 
validity  of  the  inequalities  in  (17),  the  system  is  also  guaranteed 
to  be  stable. 

The  inequalities  (26)  -  (18)  determine  certain  regions  in  the 
space  of  the  parameters  Ki  and  K2.  In  the  same  space  we  can  also 
construct  the  exact  boundary  of  the  region  of  stability  by  using  the 
method  of  D-decomposition.  Three  regions,  constructed  using  the  suf¬ 
ficient,  necessary,  and  exact  conditions  for  stability,  respectively, 
are  shown  in  Figure  1.  The  sufficient  conditions  for  stability  deter¬ 
mine  a  region,  lying  completely  inside  the  exact  region  of  stability, 
and  the  boundary  of  the  region  obtained  using  the  necessary  conditions 
for  stability,  encompasses  the  exact  region  of  stability. 

Example  2.  In  the  preceding  example,  the  approximate  tests  for 
stability  were  used  to  solve  a  problem  that  can  also  be  solved  using 
the  exact  methods  (D-decomposition).  Let  us  make  the  problem  a  little 
more  complicated.  Suppose  we  are  given  the  following  characteristic 
control  function 


2«  +  12z5  +  47Z4  +  108z3  ■+•  122eJ  +  Kk#  +  /TAT,  -  0,  ( 19 ) 

where  K  is  a  variable  parameter  of  the  object,  and  K»  and  K2  are  the 
unknown  parameters  of  the  regulator.  The  method  of  D-composition  is 
difficult  to  use  here,  since  first  of  all  the  problem  is  three- 
dimensional  which  makes  it  hard  to  visualize  the  solution.  Secondly, 
it  does  not  reduce  *0  a  system  of  three  linear  equations  with  three 
unknowns.  However,  the  use  of  the  "approximate"  conditions  for 
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stability  gives  simple  inequalities 
that  do  not  require  geometrical 
interpretation.  Thus,  the  U3e  of 
the  necessary  conditions  for 
stability  (2)  and  (3)  yields: 


0<K,<-jjr 


Kt< 


242 


(20) 


/C,>  0,89/Ci 


Figure  1.  Exact  and  "approximate" 
regions  of  stability. 

1  —  region  of  stability  obtained 
using  the  sufficient  conditions 
for  stability  (17)  and  (18);  2  — 
exact  boundary  of  stability, 
obtained  using  the  method  of 
D-decompositionj  3  —  boundary  of 
stability  obtained  using  the 
necessary  conditions  for  stabi¬ 
lity  (16). 


The  sufficient  stability  con¬ 
ditions  (13)  yield 


0  </Cj  <  0,52/C, 


(21) 


Finally,  the  sufficient  con 
dltions  (l1!)  yield 


0</C,</C,  — 0.004 Klk- 

if  >  160 
A|  '  . 


(22) 


For  any  given  range  of  change  in  the  parameter  K,  we  may  find  the 
values  of  the  parameters  Ki  and  K2  necessary  for  stability. 

Thus,  it  may  be  seen  that  the  "approximate"  stability  criteria 
may  successfully  replace  the  "precise"  stability  criteria  for  high  order 
systems.  The  use  of  the  "approximate"  stability  criteria  yields  simple 
relationships  between  the  system  parameters,  which  makes  it  possible  to 
synthesize  the  regulator  simultaneously  for  a  certain  set  of  values  for 
individual  parameters  of  the  system.  The  use  of  "approximate"  stabi¬ 
lity  criteria  greatly  increases  the  productivity  of  the  engineer- 
designer. 
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CERTAIN  PROBLEMS  OP  THE  CONSTRUCTION  OF  PARAMETRICALLY 
INVARIANT  AUTOMATIC  CONTROL  SYSTEMS 

N.  I.  Sokolov 

The  control  of  nonstationary  objects  with  their  parameters 
changing  within  wide  limits  is  in  principle  possible  with  the  aid  of 
adaptive  automatic  control  systems  and  systems  equivalent  to  them. 

Here  adaptive  regulators  will  be  defineci  as  regulators  in  which 
certain  (adjustable)  parameters  are  functions  of  the  corresponding 
variable  parameters  of  a  nonstationary  controlled  object.  Consequently, 
the  operation  of  such  a  regulator  requires,  in  principle,  that  infor¬ 
mation  be  given  about  the  variable  parameters  of  the  controlled  object. 
Technically  the  variation  of  the  regulator  parameters  is  achieved  with 
the  aid  of  parametric  couplings.  A  parametric  coupling  is  a  necessary 
(but  not  sufficient)  trait  of  an  adaptive  system.  The  basic  problem 
in  the  construction  of  adaptive  automatic  control  systems  is  to  organ¬ 
ize  the  information  about  the  variation  of  the  parameters  of  the  non¬ 
stationary  controlled  object. 

In  adaptive  control  systems  with  open  loops  for  the  adjustment 
of  the  regulator  parameters,  the  information  about  the  change  of  the 
parameters  of  the  nonstationary  object  is  obtained  directly  by  meas¬ 
uring  the  external  conditions  that  result  in  the  nonstationarity  of 
the  object,  or  in  certain  cases  by  using  a  specially  programmed  unit. 
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In  both  cases  the  accuracy  with  which  the  problem  is  solved  is 
small. 

Adaptive  control  systems  with  closed  loops  for  the  adjustment  of 
the  regulator  parameters  offer  much  greater  possibilities.  The  infor¬ 
mation  about  the  changes  in  the  parameters  of  a  nonstationary  object 
is  obtained  from  the  signals  circulating  in  a  closed  system.  However, 
a  signal,  circulating  in  a  system,  has  no  natural  characteristics 
that  would  lead  in  a  technically  simple  fashion  to  information  about 
changes  in  the  parameters  of  an  object. 

On  the  other  hand,  we  cannot  exclude  the  theoretical  possibility 
of  building  such  characteristics  into  a  signal  by  selecting  structures 
and  parameters  of  the  basic  contour  of  a  control  system. 

The  best  known  methods  of  extracting  information  about  changes 
in  the  parameters  of  an  object  are  the  correlation  methods  with  the 
method  of  the  energy  balance  as  one  of  their  forms.  These  methods 
have  become  widespread  in  solving  engineering  problems.  However,  the 
energy  balance  method  is  applicable  only  in  the  case  when  the  para¬ 
meters  of  an  object  change  very  little  during  the  transition  process. 
The  technique  in  which  the  effect  of  changing  parameters  (parameter) 
is  translated  into  a  high-frequency  signal  by  properly  selecting  the 
structure  and  parameters  of  the  basic  circuit  can  theoretically  be 
used  to  extract  information  about  changes  in  a  certain  parameter  of 
an  object.  This  is  the  so-called  indirect  method  of  determining 
parameters  of  an  object. 

The  equivalent  adaptive  regulators  will  be  defined  as  those 
regulators  in  which  the  information  about  the  variable  parameters  of 
an  object  is  not  used  and  is  not  needed  to  provide  a  given  type  of 
control. 

The  usual  dynamic  and  static  properties  of  the  equivalent  adaptive 
regulators  are  achieved  by  introducing  proper  coordinate  links.  It 
is  also  theoretically  possible  to  use  parametric  links  in  such  regu¬ 
lators.  However,  the  parametric  linkage,  in  contrast  to  adaptive 
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systems,  changes  the  parameters  of  a  regulator  depending  on  the 
instantaneous  value  of  certain  phase  coordinates  of  the  control  system. 

The  equivalent  adaptive  control  systems  may  be  subdivided  in 
two  classes: 

1)  parametrically  invariant  systems  are  characterized  by  the 
fact  that,  when  the  parameters  of  a  control  object  change,  the  output 
coordinates  of  the  performance  indicators  of  the  system  remain  the 
same  or  change  in  a  desired  direction  within  given  limits; 

2)  ultra-coarse  systems,  characterized  by  the  fact  that  —  as 
the  parameters  of  an  object  vary  —  the  stability  of  a  control  system 
is  preserved,  and  the  performance  indicators  are  not  subject  to  control. 

These  classes  in  turn  may  be  subdivided  into  groups  that  differ 
from  one  another  in  design,  and  naturally,  in  their  capacity. 

Figure  1  shows  the  classification  table. 

The  definitions  given  here  and  the  classification  table  do  not 
pretend  to  be  complete,  have  purely  utilitarian  objectives,  and  are 
used  in  the  further  discussion  of  the  material. 

Let  us  briefly  discuss  the  basic  differences  and  capabilities  of 
the  structures  in  question. 

The  parametrically  invariant  control  systems  can  be  subdivided 
into  the  following  groups: 

a)  parametrically  invariant  compensating  systems  are  linear 
within  a  limited  range.  In  order  to  assure  stability  and  independence 
of  the  performance  indicators  when  the  parameters  of  an  object  undergo 
a  change,  one  introduces  inversely  parallel  correcting  devices  which 
comprise  the  object  and  certain  elements  of  the  regulator,  and  are 
in  the  form  of  physically  realizable  filters.  In  order  to  compensate 
for  the  lag  Introduced  by  physically  realizable  filters,  a  low  fre¬ 
quency  filter — whose  poles  are  equal  to  the  poles  of  the  inversely 
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parallel  correcting  devices 
—  is  introduced  in  the  direct 
circuit  between  the  measuring 
element  and  the  adder. 


The  structural  diagram  of 
the  parametrically  invariant 
compensating  systems  in  their 
general  form  is  given  in 
Figure  2,  where  Wi(p)  is  the 
measuring  element;  VMp)  is 
the  amplifier;  Wj(p)  is  the  drive,  Wo(t,  p)  is  the  control  object; 
and  Wki(p)  are  the  correcting  devices. 

The  poles  of  the  transfer  function  Wk#(p)  are  identical  with  the 
poles  of  the  functions  W^(p),  Wkj(p),  Wkj(P)>  Wkj(p)  (°r  at 

least  identical  with  the  smallest  ones  among  them). 

A  differential  equation  for  parametrically  invariant  compensating 
systems  can  in  general  be  written  as 


Figure  2. 


n-rtn  m+m-1 

AJ-1iT~  + 

JmQ 


.2 

J-o 


2  Cj(i)t&  +  2jKpK(i)b,Zjil>*‘ 

-KPK(l)x.;, 


d'x 


(1) 


where  n  is  the  order  of  the  differential  equation  for  an  uncorrected 
system,  m  is  the  order  of  the  differential  equation  describing  the 
processes  in  the  correcting  devices,  C j ( t )  are  the  coefficients  which 
contain  the  variable  parameters  of  the  object  and  their  derivatives, 

K(t )  is  the  product  of  the  amplification  factor  and  the  natural  fre¬ 
quency  of  the  nonstationary  object.  The  introduction  of  a  compensating 
filter  minimizes  the  coefficients  multiplying  the  derivatives  in  the 
equation.  By  the  same  token,  assuming  there  is  interference  and  non- 
linearities  of  the  saturation  type,  it  keeps  the  control  system 
operational  when  the  parameters  of  the  object  vary  within  a  wide  range. 
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b)  Systems  of  variable  structure  with  internal  feedbacks 
linking  the  control  object  and  certain  elements  of  the  regulator, 
technically  based  on  the  iji-cells  Involving  sufficiently  small  lags  in 
the  correcting  devices.  Also  theoretically  assure  parametric  invari¬ 
ance  when  the  parameters  of  an  object  vary  within  a  certain  range. 

A  diagram  of  a  system  of  variable  structure  is  shown  in  Figure  3, 
where  F  is  the  unit  in  which  the  switching  functions  are  formed,  bj 
is  the  coefficient  of  action  with  respect  to  the  derivatives. 

Taking  into  account  the  real  lag  involved  in  the  formation  of 
the  derivative  signals,  the  differential  equation  for  a  system  of 
variable  structure,  shown  in  Figure  3»  may  be  written  as 

«+*»  .  n+m-\  m 

2  Aj-$tr+  2  cy(0“jir  +  S  KpK(t)bjX 

M  M 

JmO 

b. 

where  b.  *  {.J1,  Since  here  we  deal  with  systems  of  variable 
J  Ja 

structure  which  are  characterized  by  sliding  regimes,  when  a  sliding 
regime  does  occur,  the  above  equation  may  be  replaced  with  a  linear 
one,  where  the  coefficients  are  equal  to  the  averages  bt «■» , 


* 
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The  linear  equation  thus  obtained  will  describe  the  motion  in  a 
hyperplane.  Past  movements  relative  to  the  hyperplane  are  not  repre¬ 
sented  by  this  equation.  However,  since  the  frequency  of  a  sliding 
regime  in  systems  of  variable  structure  of  this  type  is  much  higher 
than  the  band  of  the  working  frequencies  of  the  system,  these  fast 
oscillations  will  have  practically  no  effect  on  the  output  coordinate. 
Consequently,  the  systems  of  variable  structure  with  physically  reali¬ 
zable  correcting  devices  may  in  this  case  be  investigated  using  linear 
equations. 

The  occurrence  of  an  uncompensated  lag,  introduced  by  physically 
realizable  correcting  df /ices,  interferes  with  the  minimization  of 
the  coefficients  of  the  differential  equation  for  the  system.  There¬ 
fore,  when  there  is  interference  and  a  system  involves  additional  non- 
linearities,  its  capacity  is  less  than  that  of  a  parametrically  in¬ 
variant  compensating  system.  In  addition,  the  nonlinear  control  part 
(the  control  part  might  be  said  to  have  a  source  of  internal  inter¬ 
ference)  also  results  in  a  higher  sensitivity  to  Interference  than 
in  the  case  of  parametrically  invariant  compensating  systems. 

A  system  of  variable  structure  may  serve  sb  an  example  of  a 
parametrically  invariant  system  widely  utilizing  parametric  links  in 
the  phase  coordinates  due  to  which  the  system  acquires  certain  non¬ 
linear  properties.  In  general,  the  presence  of  parametric  links  does 
not  necessarily  result  in  nonlinear  properties,  which  can  be  directly 
seen  in  Figures  4  and  5.  Figure  5  shows  an  equivalent  model  of  the 
control  portion  of  a  system  of  variable  structure  [2],  A  system  of 
variable  structure  is  characterized  by  the  parametric  links  in  the 
coordinates  x,  x,  x...,  resulting  in  a  change  of  the  "bank"  of  relays. 

Figure  4  shows  the  control  portion  which  is  fairly  close  to  the 
control  portion  of  a  system  of  variable  structure,  with  the  exception 
that  here  the  modulus  is  taken  of  the  total  signal  ” q .  **  The  parametric 
linkage  does  not  disturb  the  linear  properties  of  the  system.  Thus, 
in  the  case  of  Figure  4,  we  have  a  technical  realization  of  a  linear 
system  composed  of  nonlinear  elements. 


c)  Systems  involving  a  special  nonlinear  filter.  These  systems 
are  conceptually  similar  to  systems  of  variable  structure.  They 
differ  from  the  latter  in  that  a  signal  entering  a  relay  element  is 
formed  by  one  operator,  whereas  a  signal  entering  modulus-type  non¬ 
linear  elements  is  formed  by  another  operator.  In  the  case  of  a 
single  frequency  input  signal,  one  can  form  any  amplitude  and  phase 
frequency  characteristics  (in  the  first  harmonic).  However,  in  the 
case  of  a  mixed  signal,  in  view  of  the  invalidity  of  the  superposition 
principle,  these  characteristics  do  not  hold.  If  in  such  a  system 
high  frequency  sliding  regimes  occur,  then  processes  In  the  system, 
just  as  in  a  system  of  variable  structure,  are  approximately  described 
by  linear  equations.  The  capabilities  of  such  systems  in  the  case 

of  a  nonstationary  automatic  control  are  similar  to  those  of  systems 
of  variable  structure. 

d)  Auto-oscillatory  control  systems  (with  a  passive  self¬ 
adjustment).  These  systems  perform  harmonic  linearization. 

For  a  certain  structure  of  the  regulator  with  respect  to  slow 
actions,  the  system  will  operate  as  a  linear  and  stationary  system, 
independently  of  the  variable  parameters  of  the  control  object.  The 
amplitude  of  the  auto-oscillations  will  vary  depending  on  the  changes 
in  the  parameters  of  the  object. 

The  amplitude  of  the  auto-oscillations  must  for  technical  reasons 
lie  within  a  certain  range,  which  considerably  limits  the  capabilities 
of  cuch  systems. 

The  systems  in  question  are  also  more  sensitive  to  noise  than 
linear  systems. 

e)  Nonlinear  systems  with  high-frequency  external  signals. 

Such  systems  perform  linearization  at  the  expense  of  external  high- 
frequency  signals.  The  system  has  quasllinear  properties  with  respect 
to  slowly  varying  signals.  The  capabilities  of  such  systems  approxi¬ 
mate  those  of  the  auto-oscillatory  systems. 
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We  have  in  mind  the  following  systems: 

a)  parametrically  invariant  compensating  systems; 

b)  systems  of  variable  structure; 

c)  systems  with  special  nonlinear  filters. 

We  can  formally  construct  a  combined  system  of  the  form  linear- 
variable  structure  auto-oscillatcry.  However,  this  type  of  system 
is  of  no  technical  interest,  since  all  problems  are  solved  more  com¬ 
pletely  by  a  combined  system:  parametrically  invariant  compensating- 
auto-oscillatory  (optimality  at  large  deviations  and  stability,  as 
well  as  given  values  of  performance  indicators  at  small  deviations). 

To  illustrate  the  parametrically  invariant  systems  and  ultra- 
coarse  control  systems,  we  shall  use  an  example  from  [3]  (Kokhenburger's 
system).  The  equation  of  the  control  object  will  be 

^“grX+So)'  ',here  m<*<o<io  W 


The  regulator  elements  were  assumed  to  be  inertia-free  Wi(p)  »  Kj 
Wi(p)  ■  Kt  W»(p);for  simplicity  we  assume  that  Kx  *  1,  Ks  «  1. 

The  control  system  to  be  designed  must  satisfy  the  following 
technical  requirements: 

a)  static  error  no  greater  than  5%i 

b)  maximum  regulation  time  no  greater  than  t  *  0.28  sec  (the 
process  is  considered  to  be  over  when  the  dynamic  error  is  less  than 
5*); 

c)  over-regulation  is  no  greater  than  15% • 
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A  Parr metrically  Invariant  Compensating  System 

1)  Excluding  the  static  error,  we  include  an  integral  link  in 
the  adder  (See  Figure  3.)* 

Then 


In  a  general  form  Wu (p) - gg-g-;  ^(^=*5—; 


The  transfer  function  for  che  system  will  be 

_ K, 5000X,  1 1)  _____ 

A  w  "  U»+5)»(P+i!00i  fl*i  Ip)  £  +7v'7;OOOK,  (/)  + 

1 

4  k,A,  (p)  p(p+ 5)*  (p+2O'J)+tf,50K)OK‘<  (/)  A*,  (p)  p 

In  case  of  correcting  devices  of  the  hype 

0,383p»  ■}•  29 p' ±  3850 p  +  7ft  000  _ 

W  "*  (p+150)  (p  +  8C)(p+-«b)  7 


(includes  the  object,  power  drive,  amplifier), 


^«(/0 


0,135p»  +  'i,84p+252 
(p+150)(p-r£m  (P+40)T 


(5) 


(6) 


(7) 


(includes  amplifier), 

WWHW'mW  7rav//+  i50)(r +so)(p+40tp  (8) 

(The  last  correcting  dv:''  ice  was  put  between  the  measuring  element 
and  the  adder)  and  K*  *  2140,  ...1  technical  requirements  are  satisfied. 

Let  us  consider  the  operational  ability  ot  the  system  i u  the 
presence  of  interference  sr.d  saturation-type  nonlineurities  in  the 
drive . 

Suppose  that  the  drive  becomes  saturated  by  an  input  signal  xdr  * 
16  V.  Consider  xdl,  with  noise  of  the  form  0.001  sin  2000  tin» 
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a)  The  noise  xin  ■  0.001  sin  2000  t^fl  acts  on  the  input  of 
the  system.  The  object  has  a  sufficieht  inertia  and  does  not  pass  a 
signal  of  frequency  w  «  2000.  Therefore,  at  that  frequency  the  system 
might  be  said  to  be  open  with  respect  to  the  principal  feedback,  and 
we  can  write  approximately 


*np 


80  000 


*Bx.n.  (P+  1»mP+80)(p+40)  P ' 


2  M0 

(Q.135p,+0.84p  +  252)2M0' 


1+  (p+T50)(p+8O)(p+4O) 
480  000-2140 


(p»+5b0p‘+23000p  +  i020000>  p 


(9) 


|xnpj«0t001.g^^b-0,0000000615  V 


(10) 


The  noise  at  the  input  to  the  drive  will  be  greatly  weakened. 


b)  The  noise  xQUt  n  ■  0.001  sin  2000  t  enters  the  input  of  the 
correcting  device  Wc  d  1  (p). 


2M0 


Xnv  0.383P* + 29p*  +3850p  +  78  000 

*  (P+150)(P+80)(p+46) *  ,  (0.135p*+0.84p+252)  2140 

(p-fi56)(p+80)(p+40) 


(0,383p,+29p,+3850p+76  000)2MQ 
p> + SOOp’-f  &  OOOp  + 1 020000 


(11) 


|*np|»0,001 .0,368.2140-0,785  V 


(12) 


The  jystera  will  be  operational  with  this  type  of  noise. 
2)  Correcting  devices  may  be  chosen  to  be  simpler. 


For  example,  the  above  dynamic  requirements  are  satisfied  for 


Wu(p)' 


0.054p+S.38 


x^-T-n.u o.  lyt  0,0204  XV/  /  v  45 

P  +  45  1  (/*)“ 7+5  >  ^*0  (^)“ 5+45) 


45 


(13) 
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The  operational  capacity  of  the  system  in  the  presence  of  noise  is 


a) 


*HP 


45-36  SOd  (p-f-5) 

*  P(P+45)(p+800) 


|*J*  0,001  -0,38-0,00038  V 


(14) 


b) 


•fnp  _ 


w.ww-w-^w  V 


(15) 


In  the  second  version,  the  scheme  is  also  operational,  but  its 
resistance  to  noise  is  lower  than  in  the  firs*  version. 

Systems  of  variable  structure  and  the  system  with  nonlinear 
filters  without  compensating  filters  will  have  a  lower  resistance  to 
interference  than  the  above  structure. 

(Control  Systems  with  Passive  Self-Adjustment. 

Auto-Oscillatory  Control  Systems) 

Auto-oscillstory  control  systems  are  also  parametrically  invari¬ 
ant  but  their  operational  regimes  differ  from  the  preceding  ones. 


Figure  6. 


Let  us  consider  Kokhenburger's  solution  (Figure  6). 

The  equation  of  the  oscillatory  loop  for  Kj  *  300,000  will  be 

Kt  (p)-= 


_ SPOOK, ft.  (pM-60-)p»-f  ICS  1 10/)*+ 15,28-  10*p+7C0-  K*‘ _ 

|p*+O74p!-HVI77'p',TlJ2-lO‘p*+‘J47-IOVl+8l0O-IO‘p-H»O5O-IO‘HP+Z00)  (16) 
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Under  those  conditions,  the  loop  becomes  unstable.  Due  to  the 
saturation  of  the  amplifier  (or  the  drive)  the  system  will  be  auto- 
oscillatory.  The  expression  in  the  brackets  can  be  factored  into  a 
product  of  two  polynomials 

(p3+667/>*+ 167  046/>+ 130,8 . 10%33 + 7,2/>3+61/7+145)  ( 17 ) 


The  first  polynomial  has  roots  with  positive  real  parts.  The 
instability  will  occur  at  the  frequency  o>  =  510.  If,  just  as  before, 
the  saturation  occurs  at  16  V,  then  the  amplitude  of  the  auto-oscil¬ 
lations  of  the  object  will  be  less  than  the  amplitude  of  the  noise 
xout.dr  *  0,001  v*  The  system  will  not  be  operational.  One  can 
increase  the  amplitude  of  the  auto-oscillations  at  the  output  from 
the  object  by  introducing  an  amplifier  after  the  nonlinear  element, 
and  a  divider  before  the  nonlinear  element.  However,  the  result  will 
be  that  the  amplitude  of  the  auto-oscillations  at  the  input  to  the 
object  will  be  very  high,  which  is  not  permissible  under  real  conditions. 

Ultra-Coarse  Control  Systems.  A  Structure  that  Remains 

Stable  for  K  •»  «> 


The  system  has  t  *  0.28  sec  for  Ko  *  0.1;  for  Ke  *  10'  the  system 
has  t  =09.02  sec.  The  over-regulation  for  K0  0  0.1  is  15%. 

r 


Figure  7. 


The  operational  ability 
in  the  presence  of  noise. 

Here  the  effect  of  the 

noise  x .  and  x^,  .  „  is 
in .  n  out .  n 

equivalent  and  similar  to 
what  was  used  earlier 


■v„p  200J  (p  +  IOO)5  . 
~  (w+^oiuo)1  ’ 


(18) 


Ja‘, .p  j = O.CWi  •  80  U00-  80  7 


(19) 
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At  16  V  the  drive  will  become  saturated  due  to  noise.  For  this 
type  of  noise,  vibrational  linearization  will  occur,  the  effective 
amplification  coefficient  of  the  drive  will  be  reduced  several  times, 
t  will  increase,  and  in  this  case  the  process  will  become  more 
oscillatory. 
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CERTAIN  PROBLEMS  OP  STATISTICAL  ESTIMATION  IN  CONTROL  THEORY 


A.  S.  Golubkov 

The  desire  to  satisfy  the  increasing  requirements  placed  on 
technical  systems  leads  in  a  majority  of  cases  to  the  use  of  optimal 
design  methods  in  the  construction  of  control  systems. 

When  one  solves  applied  problems  in  that  area,  there  always 
arises  the  problem  of  obtaining  reliable  information  that  would  be 
sufficient  for  the  construction  of  optimal  systems.  Attempts  to 
solve  this  problem  also  encounter  certain  specific  theoretical  diffi¬ 
culties,  particularly  when  the  problem  is  statistically  formulated. 

As  a  rule,  problems  of  that  type  cannot  be  tackled  by  a  single  theory. 
They  require  the  use  of  probability  theory,  mathematical  statistics, 
and  theory  of  optimal  processes. 

The  results  obtained  in  those  areas  indicate  that  to  construct 
high-quality  systems  it  1b  absolutely  imperative  that  measurements  be 
made.  If  this  is  impossible,  the  values  of  the  parameters  character¬ 
izing  the  state  of  the  processes  investigated  must  be  calculated  with 
high  accuracy  and  with  no  substantial  delay.  The  occurrence  of  various 
random  disturbances  and  measurement  errors,  and  changes  in  the  para¬ 
meters  and  characteristics  within  a  wide  range  in  a  way  impossible  to 
predict,  force  us  to  use  the  statistical  approach  in  the  evaluation 
of  processes.  As  always,  the  estimates  are  a  result  of  certain  opera¬ 
tions  on  the  measurement  date,  permitting  us  to  estimate  the  state 
of  the  process. 
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The  state  of  a  process  can,  as  we  know,  be  described  by  a  stochas¬ 
tic  model  to  any  degree  of  accuracy.  Moreover,  the  correspondence 
between  the  model  and  the  real  process  occurs  only  if  the  model  is 
statistically  sufficient.  Considering  only  one  class  of  processes 
—  namely,  the  class  of  controlled  Markov  processes  —  we  shall  assume 
below  that  their  space-time  behavior  is  described  by  the  trajectory 
x(t)  satisfying  a  model  written  in  the  form  of  a  stochastic  differential 
equation 


.v-/7(jc.0  +  ?(Jf. />»(/)  (1) 

where  to ( t )  is  a  random  process  of  the  type  of  white  noise  with  the 
parameters  E{w(t)}  and  E{<o(t)w(T)}  (E  is  a  symbol  denoting  the  aver¬ 
aging  operation). 

Information  about  the  state  of  the  coordinate  x(t)  can  be  obtained 
by  measuring  a  certain  y(t)  which  is  related  to  the  original  x  by 
equations  of  the  type 


y  (/)-//(*,/) 


(2) 


or 


y  (t)  =  H  (x,  t,  «>)  ( 3 ) 

and  is  given  as  a  solution  of  the  stochastic  differential  equation 

k(t)^A(z,y,t)  +  D(y,l)v(t)  (4) 

where  v(t)  is  a  random  process  of  the  type  of  white  noise  with  the 
parameters : 


£  {v  (t)}  0  and  £  { v  (()  v  (x)) 

This  type  of  description  allows  us  to  take  into  consideration 
the  peculiarity  of  the  measuring  process  which  includes  the  dynamic 
properties  and  the  additive  noise  of  measurements. 
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Assuming  for  simplicity  that  x,  y,  and  z  are  scalar  processes 
and  H(x,  t)  =  x,  the  equation  of  measurements  (4)  becomes 

z(t)**A(x.z.t)  +  lHx,t)v(t)  (5) 

According  to  the  terminology  accepted  in  theory  of  controlled 
random  processes  [1],  the  set  (x(t);  z(t)},  where  x(t)  is  the  unob¬ 
served  coordinate,  and  z(t)  is  the  observable,  came  to  be  called  par¬ 
tially  observable,  and  the  procedures  based  on  the  use  of  the  second 
coordinate  z(t)  alone  came  to  be  called  the  procedures  based  on 
incomplete  data. 

In  our  case,  the  set  (x(t);  z(t)>  satisfying  (1)  and  (5),  repre¬ 
sents  a  two-dimensional  Markov  process  in  which  the  first  coordinates 
x(t),  z(t),  being  fixed,  form  the  so-called  conditional  Markov  process. 

The  problem  as  applied  to  a  conditional  Markov  process  belongs 
to  the  class  of  problems  of  mathematical  statistics,  and  reduces  to 
the  question  of  how  one  can  best  estimate  the  component  x(t)  on  the 
basis  of  measurements  of  z(t).  The  technical  solution  of  the  problem 
involves  processing  large  amounts  of  data,  and  acceptable  solutions 
can  be  obtained  only  if  one  extracts  ’’sufficient  statistics"  based  on 
the  set  of  the  measurement  data  y(z),  that  are  the  central  question 
in  those  problems  [2]. 

One  of  the  principal  difficulties  in  this  approach  is  the  volume 
of  measurement  data,  which  increases  as  time  goes  on,  and  whose 
processing  in  real  cime  and  in  full  is  difficult  using  the  computers 
available  today. 

In  a  majority  of  problems  of  sequential  analysis,  where  the 
present  estimate  problem  also  belongs,  the  a  posteriori  probabilities 
turn  out  to  be  sufficient  statistics. 
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Referring  to  [3],  we-  can  state  that  the  use  of  the  Markov 
properties  of  the  set  {x(t);  z(t)}  enables  us  to  considerably  reduce 
the  volume  of  data  required  to  obtain  sufficient  statistics.  In  the 
case  of  a  finite  observation  interval,  the  a  posteriori  probability 
is  a  Markov  sufficient  statistic  if  the  following  conditions  are 
satisfied : 

1)  unobserved  process  x(t)  is  a  Markov  process; 

2)  observed  process  consists  of  a  sequence  of  conditionally 
independent  random  variables; 

3)  optimality  test  is  additive; 

4)  preceding  measurements  and  estimates  do  not  impose  any 
restrictions  on  the  following  estimates. 

If  the  conditions  (3)  [33  and  (4)  [4]  are  satisfied,  we  can  con¬ 
sider  the  synthesis  of  optimal  estimates  as  an  extremal  problem  of 
mathematical  statistics  in  relation  to  a  conditional  Markov  process 
subject  to  conditions  (1)  [1]  and  (5)  [53* 

The  problem  will  consist  of  making  decisions,  *(/)» t (?(•))  xOX. 
constructed  out  of  known  data  z(t),  zEZ,  0  <  t  <  T,  which  must  be 
optimal  in  the  sense  that  they  must  give  the  best  estimate  of  the 
parameter  x(t),  xEX,  t  v  T. 

When  the  Bayesian  approach  is  used,  such  an  estimate  (here  we 
have  in  mind  a  point  estimate)  must  minimize  the  mean  risk  [43 

R  --  )lo  .v, :  (:))  p  ( x ,  z)  ilxclz  ( 6 ) 

where  l o  [z,  x,  y  (z)]  is  a  loss  function. 

To  satisfy  (6),  it  is  sufficient  to  minimize  the  conditional  risk 

r  (x/zj  .-\/0  (z,  x,  x)  p  (xiz)  (lx  ( 7  > 


36 


at  each  point  individually. 


The  conditions  (1)  and  (5)  on  the  set  (x(t);  z(t)}  require  that 
one  consider  additional  conditions  in  the  extremal  problem,  in  parti¬ 
cular,  that  one  introduce  a  new  loss  function  L(x,  z,  x) 


n—\ 


L  ( X ,  z,  x)  -  l0  (x,  z,  *)  +  2  Vi  (•*.  X) 

/sat 


(8) 


and  minimize  the  mean  risk  R(x)  under  the  condition  that 

\ll(z,x,x)p(x,z)dxdz~cl  (9) 

In  view  of  (8)  and  (9)»  we  assume  that  the  estimate  has  the 
form  x  *  -y(z,  A),  where  the  parameter  A^  ■  A^(c^)  is  determined  by 

5  /|  (AT,  2,  T  (2,  A,))  p  ( X ,  2)  dxdz  — ct  (10) 

If  System  (10)  has  a  solution,  then  the  estimate  x  *  y(z,  A(c)) 
is  optimal  for  the  loss  function  (8)  and  is  the  desired  solution  of 
the  extremal  problem.  Next,  assuming  that  the  measurement  interval 
is  bounded  0«=<.t<.T<«,  and  the  process  in  which  we  obtain  and 
process  the  data  is  sequential  in  time,  the  introduction  of  time  in 
the  problem  leads  to  the  following  relations  for  the  loss  function  L 
and  the  estimate  x 

r 

L(x,z,  x,  7')  —  jj  L  {x  ( t ),  2  (f),  x  (0)  dt  +  Lr  (xr,  zt<  xt) 

o  (11) 


i(7’)  =  Ir(zW)  Ov) 

respectively.  The  optimal  solution,  satisfying  (11)  and  (12)  for- 
fixed  T,  can  be  found  from  the  minimum  condition 

ay (.v (T)/z  (-.))  =>  ^  L  (x (t),  2  (-),  x (t))  p (a*  (-))  dx 
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The  conditional  risk  (13)  depends  on  z(t)  0  <  t  <  T  through  the 
a  posteriori  distribution  p(x(t)/z(i))  0  .<  T  <  T,  which  contains  the 
entire  necessary  information  about  the  structure  of  the  optimal 
estimates. 

Since  (x(t);  z(t)}  is  a  conditional  Markov  process,  p(x(t)/z(i)) 
may  be  found  for  all  t  <  T  by  applying  the  theory  of  conditional  Markov 
processes  [5].  The  evolution  of  the  a  posteriori  probability  p(x/z) 
is  obtained  as  a  solution  of  the  stochastic  differential  equation, 
and  for  the  model  (1)  and  (5)  it  is  characterized  by  the  solution  of 
an  infinite-dimensional  system  of  stochastic  differential  equations 


i>x,t  it)  ~  Px/t  (t)  |i(0-  A  {»,  01  I* 

+{-,T  if*  w  Pf  n  +  t  5  w  <*.  v  p*  n] 

A  (z,  x ,  t)  -  J  A  (x,  z,  t)  px„  (t)  dx 


Suppose  that  P*/Z(t)  is  a  solution  of  (14).  Then  we  shall  assume 
that  the  estimate  of  x(t)  will  have  the  form  x(z,  t)  «  Y#(p*/Z(t). 

We  know  that  in  the  case  when  the  probability  density  Px/Z(t)  is 
uni modal  and  symmetric  in  its  expectation,  and  the  loss  function  is 
symmetric,  i.e.,  L(x,  x)  ■  L(x,  x),  then  the  conditional  expectation 
is  the  best  estimate.  Therefore,  the  optimal  estimate  will  always  be 
chosen  as 


x  *=  7  (?)  =  [  xp  (x/z)  dx  «= 
* 


J  xp  (x)  p  (z/x)  dx 
\p{x)p(z{x)dx 


(15) 


Such  an  estimate  being  asymptotically  unbiased  represents  an 
efficient  estimate  in  the  sense  that  the  deviation  is  minimized  [2,  4]. 

The  above  expression  for  the  a  posteriori  probability  clearly 
shows  that  the  process  of  finding  estimates  requires  a  solution  of 
parabolic  partial  differential  equations,  which  is  very  difficult. 
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Since  Equation  (14)  is  infinite-dimensional,  we  are  actually  not 


interested  in  it  as  much  as  in  the  fact  that  from  it  one  can  obtain 
equations  of  certain  characteristics  of  the  process  Px/Z(k)>  for 
example  moments,  characteristic  functions,  and  semi-invariants. 


However,  even  with  this  simplification,  the  equations  for  the 
first  two  moments  can  be  obtained  and  solved  only  for  an  exponential 
family  of  Markov  process  distributions,  and  then  only  in  the  case  when 
it  is  possible  co  establish  a  one-to-one  correspondence  among  the 
sets  X,  X  and  Z. 

As  an  example,  let  us  consider  a  Markov  process  described  by 
the  equation 


x=»  c  +  ma(/)  (IE) 

xo  *  x(0)  is  a  Gaussian  random  variable  with  the  parameters  E{x)  and 
E{x2}.  u>(t)  is  a  random  process  of  the  white  noise  type  with  the 

parameters  E{w}  «  0;  E{u2}. 

The  process  x(t)  is  observed  in  the  presence  of  white  noise  by 
a  system  of  measurements  described  by  the  equation 

* -*(/)  +  *>(/)  (2E) 

v(t)  is  a  random  process  of  the  white  noise  type  with  the  parameters 
E(vl  «  0  and  Etv2}. 

The  processes  w(t)  and  v(t)  are  statistically  independent.  On 
the  basis  of  the  measurements  z(t),  0  £  t  <.  T,  we  are  required  to 
obtain  optimal  estimates  of  the  coordinates  x(t),  which  are  optimal 
in  the  sense  that  they  minimize 


(2Ea) 


89 


Per  this  purpose,  we  have  to  determine  the  first  and  second 
moment  of  the  a  posteriori  density  Px/Z(x»  t). 


The  first  moment 


xto-\xp*$)dx 


(2Eb) 


is  an  optimal  estimate,  and  the  second 


(2Ec ) 


is  its  minimum  deviation. 


The  a  posteriori  probability  according  to  (14)  satisfies,  in  the 
case  of  Example  IE  and  2E,  the  stochastic  differential  equation 


dPx/,0) 

37 


(3E) 


with  the  initial  conditions  px/i(x,  0/0)“  p~rcxp  [—  • 

The  first  term  on  the  right-hand  side  of  (3E)  is  the  right-hand 
side  of  the  direct  Kolmogorov  equation,  and  the  second  is  a  result  of 
measurements. 


Equation  (3E)  for  px/,z(t)  yields  an  equation  for  x(t)  and  D(t) 

x(/)~c  +  D(t)(i(t)~x(i))  (i»  5) 

Z>(/)=.as-^  C»E) 
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The  solution  of  the  problem  is  thus  given  by  a  system  of  nonlinear 
differential  equations  whose  analytic  solution  in  closed  form  may  be 
obtained  only  in  certain  particular  cases.  Nevertheless,  analog  and 
digital  computers  can  be  used  without  too  great  difficulty  in  conjunc¬ 
tion  with  an  algorithm  for  the  processing  of  data.  Using  digital 
computers  the  operations  on  the  data  are  in  the  form  of  recursive 
computational  procedures,  using  at  each  step  only  the  current  data, 
which  means  that  the  required  estimates  can  be  obtained  with  a  delay 
equal  to  the  machine  time  of  computation  for  one  measurement. 
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AN  ANALYSIS  OP  A  TYPICAL  STRUCTURE  OP  AN  AUTOMATIC  CONTROL 
SYSTEM  AND  A  METHOD  OP  SELECTING  THE  TRANSFER  FUNCTION 
POR  A  STANDARD  MODEL  OP  A  SELF-ADJUSTING  SYSTEM 

N.  A.  Shokalo 

In  designing  a  self-adjusting  control  system,  one  has  to  analyze 
the  original  control  system  to  discover  the  possibilities  of  satis¬ 
fying  the  T.  T.  Z.*  with  respect  to  the  performance  indicators  of 
the  transfer  process,  which  may  be  in  the  form  of: 

1.  Maximum  over-regulation  of  the  bank  angle  (y )  no  greater 
than  +15*  of  Yglv* 

2.  The  length  of  the  transfer  process  no  greater  than  2-3  sec 
with  a  deviation  of  the  bank  angle  from  Its  steady  value  by  no  more 
than  ±5*. 

3.  Maximum  over-regulation  of  the  bank  rate  (wx)  no  greater 
than  +50*. 

Methods  of  the  theory  of  sensitivity  can  be  used  to  solve  this 
problem  for  a  system  of  any  order,  but  the  solution  is  hindered  by  the 
necessity  of  choosing  a  reference  transfer  function  for  the  control 
system  [1,  3]. 

*  Translator’s  Note:  Expansion  unknown. 


« 
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'Ue  purpose  of  this  paper  is  to  select  a  reference  transfer 
function  for  a  control  system.  For  this,  it  is  necessary  to  establish 
a  relation  between  the  performance  indicators  for  the  transfer  process 
with  a  distribution  of  the  poles  and  zeros,  and  the  coefficients  of 
the  transfer  function  for  a  typical  control  system  of  fifth  degree  in 
the  denominator  and  third  degree  in  the  numerator. 

„  AtP'  +  Atf'+Atp+A’.  (2) 

3c  r+/l)p‘+/»«p*-Mlp»+/l<p+/ll 


To  solve  the  problem,  we  have  to  vary  simultaneously  from  5  to 
8  coefficients  of  the  transfer  function.  The  second  problem  to  be 
solved  consists  of  selecting  the  optimal  coefficients  of  the  control 
system's  regulator  and  the  reference  model  such  that  certain  quality 
tests  will  be  minimized. 


Mutual  Relationship  amo 

Function  on  the  Boundary  of  the  Region  of  Transfer 
Processes  of  Given  Quality 


The  transfer  function  Is  studied  using  the  method  of  mathematical 
models.  The  problem  to  be  solved  can  be  simplified  by  fixing  coeffi¬ 
cients  in  the  denominator  A»,  A*,  A*,  which  contain  only  the  time  con¬ 
stants  of  the  object  and  the  director,  or  by  fixing  the  coefficients 
A*  and  As,  containing  the  typical  regulator  and  control  surface 
efficiency  coefficients.  The  coefficients  are  varied  successively 
from  the  average  value  in  both  directions,  and  the  variation  is  limited 
by  tolerances  in  the  given  quality  characteristics.  Experiments  serve 
to  establish  the  range  of  the  coefficients  of  a  normalized  and  unnor¬ 
malized  tranfer  function  of  the  system,  and  thus  we  could  construct 
the  grids  of  the  regions  of  transfer  processes  of  given  quality  in 
the  coordinates  of  any  two  denominator  coefficients.  The  effect  of 
each  numerator  coefficient  on  the  dimensions  of  the  regions  was 
determined,  and  the  restrictions  on  the  coefficients  were  selected. 


Figure  1.  Figure  2. 


The  relationships  thus  obtained  permit  us  to  claim  that: 

1.  The  transfer  processes  satisfy  a  given  performance  level 
for  a  wide  range  of  the  first  three  coefficients  Ai,  At,  At  (Figures 
1  and  2).  For  comparison,  we  plotted  the  volumes  "l"  and  "2,"  which 
define  the  maximum  possible  variations  of  the  coefficients  Ai,  At,  At 
depending  on  the  time  constant  of  a  real  director. 

2.  The  "lower"  boundary 
of  the  region  with  the  two 
last  coefficients  in  the  deno¬ 
minator  An  and  At  being  fixed 
is  determined  by  the  high  fre¬ 
quency  oscillation.  In  order 
to  assure  a  given  level  of 
performance,  the  coefficients 
At,  At,  At  must  be  selected 
such  that  the  boundary  of  the 
high  frequency  instability 
will  shift  beyond  the  area 
determined  by  the  performance 
indicators,  since  for  values  close  to  those  on  the  boundary  the  region 
of  performance  shifts  toward  the  origin  and  has  small  dimensions  as 
compared  with  the  stability  region. 


Figure  3. 
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3.  The  region  of 
transfer  processes  of  given 
performance  may  be  widened  in 
two  ways: 

a)  by  increasing  all 
three  coefficients  Aj,  A2, 

A*  simultaneously  (Figure  3); 

b)  by  fixing  the  first 
coefficient  Aj  and  simultan¬ 
eously  increasing  the  coef¬ 
ficients  A2,  Aj  (Figure  5). 


*1.  It  is  found  that  an 
increase  of  the  coefficient 
A i  in  the  coordinates  A*,  A* 
widens  the  region  of  transfer 
processes  of  given  performance 
only  in  the  direction  of  the 
coefficient  A$,  and  a  correc¬ 
tion  to  increase  the  coeffi¬ 
cient  Aj  is  necessary  (Figure 
6). 


Figure  5. 


Figure  6a,  b. 


5.  When  analyzing  the 
functional  relationship  along 
the  boundary  of  the  high- 
frequency  oscillation  between 
the  coefficients  A2  -  A*  and 
A i  -  As,  we  find  a  relationship 
which  is  close  to  linear.  Dy 
expressing  this  relationship 
as  a  formula,  one  can  reduce 
the  number  of  the  independently 
varied  coefficients  of  the 
equation.  The  modeling  scheme 
is  shown  in  Figure  7. 
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Function  Inside  the  Region  of  Given  Performance 


On  the  grids  of  the  selected  regions  of  transfer  processes  of 
given  performance,  we  have  plotted  the  regions  where  the  roots  are 
located  and  constructed  the  boundary  curves  that  separate  them  (Figures 
3  and  4).  In  constructing  the  regions  with  the  roots,  we  have  used  a 
procedure  described  in  [2].  The  classes  of  the  roots  were  determined 
as  proposed  by  A.  A.  Fel’dbaum:  2A  Is  a  class  with  one  real  root  and 
two  pairs  of  complex  roots,  the  nearest  complex  roots  with  respect 
to  the  imaginary  axis.  2B  Is  a  class  with  the  real  root  nearest  to 
the  Imaginary  axis  and  two  pairs  of  complex  roots.  0|  Is  a  class  with 
three  real  roots  and  one  pair  of  complex  roots,  as  well  as  the  real 
root  nearest  to  the  Imaginary  axis;  0  Is  a  class  in  which  all  five 
roots  are  real. 

The  method  is  essentially  as  follows.  If  the  characteristic 
equation  of  a  closed  system  has  the  form 


V  =  u,  where  4 «  1  (?) 

•) 

then  each  set  of  coefficients  Ak  (k  *  1 ,  2  . . .  n)  is  associated  with 
a  single  class  of  distribution  of  roots  of  Equation  (2).  We  know  from 
higher  algebra  that  the  necessary  and  sufficient  condition  for  Equation 
(2)  to  have  r  identical  roots  is  that  the  polynomial  in  (2)  and  its 
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r  h 

(r  -  1)  derivative  vanish  simultaneously  when  the  multiple  root  is 
substituted  in  the’".  To  simplify  the  problem,  we  set  r  =  2.  Thus, 
when  the  multiple  root  is  substituted,  both  the  polynomial  in  (2)  and 
its  first  derivative  should  vanish.  Thus,  we  have  two  equations 

«  n 

2  Ail1*  ***0;  2  A,(/l =» 0,  (3) 

M>  *.-0 


where  Ao  *  1. 


Letting  take  on  values  from  0  to  -»  we  obtain  various  sets  of 
the  values  of  coefficients  satisfying  System  (3)  when  no  is  substituted 
in  the  latter.  Each  of  the  sets  is,  in  the  hyperspace  of  parameters, 
a  figurative  point  lying  on  the  boundary  separating  the  subregions 
with  a  different  combination  of  roots.  The  boundary  of  the  subregion 
with  a  different  combination  of  roots  is  constructed  in  the  plane  of 
any  two  parameters,  where  the  remaining  parameters  are  fixed,  i.e., 
one  intersects  the  hyperspace  of  parameters  by  a  plane  of  two  para¬ 
meters  Aj  and  Am.  One  derives  an  expression  for  the  boundary  curve 
[4]  separating  the  subregions  with  a  different  combination  of  roots 
in  the  canonical  form,  and  its  behavior  is  considered  by  studying  the 
derivative  of  the  boundary  curve.  The  shading  of  the  boundary  curve 
from  uc  *  0  to  no  ■  -*  is  done  from  the  left  if  (m  -  l)  is  e'*en,  and 
from  the  right  if  (m  -  l)  is  odd  in  case  the  wm"  coefficient  lies  on 
the  axis  of  ordinates. 

The  distribution  of  the  zones  of  roots  in  the  coordinates  of  any 
two  normalized  coefficients  permits  us  to  establish  that  the  coeffi¬ 
cients  Aa,  Aj,  A*  have  the  greatest  effect  on  the  redistribution  of 
the  zones  of  the  curve  inside  the  selected  region.  It  is  found  that 
in  the  coordinates  A*,  Aj,  which  may  first  of  all  include  the  para¬ 
meters  of  the  regulator,  a  wide  performance  region  can  be  obtained 
in  two  ways: 

1.  If  the  director  has  a  large  time  constant  and  the  values  of 
Ai,  A 2 ,  A3  are  close  to  the  boundary  of  the  high-frequency  oscillation, 
then  fixing  the  coefficient  Aj  one  can  correct  the  coefficients  h2 
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end  A3  in  such  a  way  that  the  boundary  of  the  high-frequency  oscil¬ 
lation  will  move  beyond  the  region  of  transient  processes  of  given 
performance.  In  this  case,  the  values  of  the  coefficients  A2  and  A3 
should  be  equal  to  A2  *  2000,  As  ■  30,000.  The  region  of  the  transient 
processes  is  here  filled  with  roots  of  the  2B  type  (with  the  real 
root  nearest  to  the  imaginary  axis  and  two  pairs  of  complex  conjugate 
roots).  The  region  of  the  roots  02  is  small  and  is  close  to  the 
coordinate  origin  (Figure  6). 

2.  A  wide  region  of  transient  processes  of  given  performance- 
can  be  obtained  by  simultaneously  increasing  the  values  of  three 
coefficients  Aj,  A2,  A«,  i.e.,  by  bringing  their  values  closer  to  the 
values  corresponding  to  the  drive  constant  Tdr  “  0  .03  sec  and  the 
booster  constant  Tfe  *  0.015  sec.  In  this  case  the  entire  region  of 
given  performance,  with  the  exception  of  a  narrow  lower  strip,  is 
filled  with  roots  of  type  0!  (where  the  real  root  turns  out  to  be 
nearest  to  the  Imaginary  axis,  followed  by  a  pair  of  complex  conjugate 
roots,  and  by  two  real  roots).  It  is  important  to  know  that  by  varying 
the  values  of  Ai,  A2,  At  one  can  superpose  the  regions  of  the  distri¬ 
bution  of  03  roots  on  the  region  of  transient  processes  of  given 
performance  (Figure  5).  Thus,  having  the  grids  of  transient  processes 
of  a  given  performance,  and  given  the  aerodynamics  and  parameters  of 
the  real  director  of  an  aircraft,  one  can  determine  in  each  case  the 
values  of  the  nominal  coefficients  of  the  transfer  function,  the 
nominal  distribution  of  roots,  and  thus  the  nominal  transfer  function 
of  the  fundamental  system,  which  may  be  used  as  the  transfer  function 
of  the  reference  model  of  a  self-adjusting  system. 

Tests  of  the  Optimization  of  the  Coefficients  of  a  Control 
System  Regulator  and  a  Reference  Model 

The  optimization  of  the  coefficients  of  the  control  system  regu¬ 
lators  and  the  reference  model  of  self-adjusting  systems  is  accomplished 
by  calculating  the  partial  derivatives  of  the  phase  coordinates  in 
•„erm3  of  the  parameters  of  the  regulator  in  the  two  basic  regimes  of 
the  ccm-v rol  system:  damping  (pilot's  manual  control)  and  stabiliza¬ 
tion.  Partial  derivatives  can  be  best  calculated  by  using  the  method 
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of  che  sensitivity  function  models  which  may  be  linear  or  nonlinear. 

The  advantage  of  the  method  is  that  the  extremal  values  of  the  coef¬ 
ficients  are  obtained  fast,  and  it  is  not  necessary  to  vary  those  para¬ 
meters  by  trial  and  error.  The  regulator  coefficients  are  optimized 
in  the  presence  of  two  actions:  control  Ygiv»  and  wind  (fw).  As 
the  optimization  test,  we  propose  a  square  integral  test  determining 
the  minimum  loop  error  in  the  stabilized  regime 

T 

/  — J«? <*.  !l>  i) dt ,  Where  E, - Tjm _ 7  (k) 

y  is  the  regulator  coefficient  in  the  damping  loop,  i  is  the  regulator 
coefficient  in  the  stabilization  contour.  The  gradients  will  be 

T 

(5) 


r 

V/-2JE ,(*,  i)di  (6) 

The  extremal  values  of  the  coefficients  u*  and  i*  for  equal  values 
of  integral  estimates  (Figure  8)  correspond  to  the  values  of  y,  lying 
in  the  selected  region  of  transient  processes  of  given  quality,  and 
the  values  of  i*  lie  beyond  the  performance  region  —  for  example,  in 
the  middle  of  the  stability  region.  With  wind  disturbances  (of  the 
same  spectrum  as  the  control  ones)  the  optimal  values  of  y*,  i*  lie 
at  the  u^per  boundary  of  the  stability  region  (Figure  12).  Thus,  the 
test  determines  the  minimum  of  the  stabilization  error.  For 

a  flight “in  a  turbulent  atmosphere,  this  test  may  result  in  consider¬ 
able  overloads  acting  on  the  pilot  and  large  values  of  vane  deflections. 
For  aircraft,  the  damping  loop  may  usually  be  considered  to  be  auto¬ 
nomous  by  using  its  integral  performance  test  to  optimize  its  coef¬ 
ficient  y.  This  is  advisable  if  one  wants  to  adjust  the  coefficient 
y  independently  of  the  adjustment  of  the  stabilization  loop  coeffici¬ 
ent  i . 
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In  this  case,  we  can  take  the  following  Integral  tests: 

T 

/ - |  E  lilt  where  e^  -  X,„  -5  ( 7 ) 

T  T 

/  **  „dt  /  —  J  (u)J’d/ 

The  phase  coordinate  Ey  characterizes  the  loop  error  for  manual 
control,  the  second  phase  coordinate  8^  characterizes  the  deflection 
of  vanes,  and  is  of  practical  interest  since  it  can  be  measured  when 
forming  the  integral  estimate.  The  third  phase  coordinate  u>x  charac¬ 
terizes  the  minimum  overloads  acting  on  the  pilot.  In  the  case  of  a 
sufficiently  fast  director  (fdr  -  4  -  5  Hz),  these  estimates  turn  out 
to  be  close  (Figures  9»  11).  The  values  of  y*  and  i*  for  the  tests 
being  minimized  lie  within  the  region  of  transient  processes  of  given 
performance.  In  addition,  the  value  of  y*  is  inversely  proportional 
to  the  effectiveness  of  the  vanes  as  the  flight  regimes  change. 

T 

The  advantages  of  the  practical  test  dt  are  as  follows: 

1.  Independence  of  the  adjustment  of  coefficients  y  and  i. 

2.  Sufficient  stability  margins  relative  to  y  and  i,  corres¬ 
ponding  to  the  desires  of  the  pilot  'in  the  linear  regime), 
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Figure  11. 


3.  Constancy  of  the  stability  margins  as  the  regimes  of  flight 
change  for  both  the  control  (Ygiv»  Xf)  and  wind  disturbances  (of  the 
same  spectrum).  (Figures  9,  10,  and  12). 

4.  The  possibility  of  realizing  a  fairly  simple  self-adjusting 
scheme  (when  adjusting  the  coefficients  u  and  i  by  different  tests). 

5.  The  average  deflection  of  vanes  is  minimized,  which  is 
implied  by  the  test  itself. 

The  shortcomings  of  the  test  include  a  lower  accuracy  of  regu¬ 
lation  in  the  stabilized  regime,  and  a  higher  necessary  value  of  the 
coefficient  u  in  the  nonlinear  regime  (Figure  13). 
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Thus,  from  the  point  of  view  of  satisfying  the  requirements  on 
the  processes  of  given  performance,  of  minimizing  the  loads  acting  on 
the  pilot,  and  the  deflections  of  the  vanes,  we  see  that  the  efficient 
test  involves  optimizing  the  coefficients  of  the  automatic  pilot 

T,  T 

l  =  j  Eipdt  or  /-=  j  ( 8 ) 

r 

In  the  stabilized  regime  during  bombing,  the  test  I  I'idt  is  efficient 

0 

or  (in  the  case  of  wind  disturbances  or  nonlinear  control  systems)  its 
modified  version  is 

+  (9) 

To  achieve  independent  adjustment  of  the  coefficient  of  the  damping 
loop  regulator,  it  is  advisable  to  use  different  performance  tests  in 
the  stabilized  regime  for  the  adjustment  of  the  coefficients  p  and  i. 
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AN  APPROXIMATE  PROCEDURE  FOR  DESIGNING  A  SEARCHLESS  SELF-ADJUSTING 
SYSTEM  SENSITIVE  TO  THE  PROPER  FREQUENCIES  OF  THE  BASIC 

CONTROL  LOOP 

B.  V.  Viktorov,  G.  N.  Izmaylov,  B.  V.  Kirsanov,  V.  A.  Pokhvalenskiy 

The  analysis  and  synthesis  of  complicated  control  systems  is  more 
often  than  not  extremely  laborious  or  completely  intractable  analyt¬ 
ically.  For  this  reason,  it  is  of  great  practical  importance  to 
develop  approximate  procedures  based  on  permissible  simplifications. 

One  of  the  possible  ways  of  investigating  nonperiodic  processes 
in  complex  systems  is  the  method  of  separation  of  composite  motion. 

It  is  based  on  the  theory  of  differential  equations  with  a  small 
parameter  multiplying  the  derivatives.  Its  application  enables  us 
to  lower  the  order  of  the  original  system  of  equations  and  estimate 
the  accuracy  of  the  approximate  solution  [5]. 

The  system  under  investigation  is  shown  in  Figure  1.  This  system, 
known  in  the  literature  as  the  Marx  system  [1],  was  studied  in  quite 
a  few  papers  [  2,3,4  ]»  and  is  used  at  the  present  time  in  aircraft. 

In  this  paper,  we  do  not  investigate  the  dynamics  of  the  system, 
its  stability,  or  the  performance  of  the  adjustment  processes.  We 
shall  only  be  Interested  in  the  possible  ways  of  obtaining  the  analytic- 
relationship 

fir  ==  '*  0  v.i  >ui  '!  u.  7  n.  5m  7  6>  7  5,„  T)  (  2  ) 
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Figure  1. 


where  Madj  is  the  adjusted  value  of  the  amplification  coefficient, 

To,  Tdr,  Tup,  T1(JW,  T  are  the  time  constants  of  the  object,  drive, 
upper  filter,  lower  filter,  and  forming  filter,  respectively;  CQ,  5dr» 
£Up»  tlow  are  the  damping  coefficients  of  the  objecjt*  drive,  lower 
filter,  and  upper  filter,  respectively. 


To  solve  the  problem,  we  propose  to  make  a  number  of  simplifica¬ 
tions  that  reduce  to  the  following: 


1.  As  shown  by  simulation,  stand,  and  flight  tests  of  this 
type  of  system,  wind  disturbances  have  hardly  any  effect  on  the  adjust¬ 
ment  processes  and  on  the  steady-state  values  of  the  coefficient  m 
(with  the  possible  exception  of  the  takeoff  and  landing  regimes). 
Therefore,  the  pilot’s  signal  in  the  course  of  piloting  was  taken  as 
the  test  signal. 


?.  The  disturbance  used  as  the  input  into  the  control  system 
was  taken  in  the  form  of  a  single  function.  Stand  tests  showed  the 
correctness  of  this  simplification.  The  error  in  the  adjusted  value 
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Figure  2. 


of  the  coefficient  under  various  actions  (also  including  manual 
piloting)  did  not  exceed  10<. 


3.  The  dynamics  of  damping  loops  of  present-day  aircraft  is 
such  that  one  can  distinguish  two  qualitatively  different  types  of 
flight  regimes: 

a)  regimes  with  a  fairly  large  scatter  of  the  roots  of  the 
drive  and  the  aircraft.  The  root  hodograph,  illustrating  this  type 
of  regimes,  is  shown  in  Figure  2.  In  this  case,  the  control  law 
assumed  here  is  satisfied  with  sufficient  accuracy; 

b)  regimes  with  a  small  scatter  of  the  roots  of  the  drive  and 
the  aircraft.  In  this  case,  the  control  law  6  ■  yu  might  be  said  to 
change  to  6  *  not.  Here  6  is  the  angle  of  deflection  of  the  control 
surfaces  of  the  object,  w  is  the  angular  velocity,  a  is  the  angle  of 
attack.  This  type  of  regimes  is  illustrated  by  the  root  hodograph 
in  Figure  3. 


Regimes  of  this  type  do  not  satisfy  the  basic  requirement  on 
the  control  system  —  namely,  that  the  drive  should  act  faster  than 
the  control  object. 
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Figure  4.  Figure  5. 

Recently  there  has  been  a  tendency  to  increase  the  speed  of 
operation  of  the  drives.  In  view  of  this  fact,  we  shall  consider 
only  the  "slow"  regimes  or  "fast"  ones  in  which  filters  or  special 
control  laws  are  used  to  obtain  differences  In  motion. 

The  procedure  we  propose  is  based  on  extracting  from  the  general 
transfer  function  (of  fifth  degree)  two  component  second-degree 
functions  which  reflect  the  characteristic  types  of  motion:  aircraft 
("slow"  and  drive  ("fast"). 

Strictly  speaking,  three  types  of  motion  occur  In  the  loop  under 
consideration.  In  the  order  of  speed  of  action,  we  have  the  aircraft, 
drive,  and  booster  motions. 

The  fastest  motion  (booster  motion)  Is  not  taken  into  account 
when  considering  the  aircraft  motion.  When  considering  the  drive 
motion,  the  booster  motion  may  be  taken  into  account  in  the  form  of 
several  terms  in  the  expansion  of  the  transfer  function. 

The  zero-order  approximation  is  used  as  the  transfer  function 
generating  the  "slow"  component  of  composite  motion.  Structure  1  in 
Figure  4  illustrates  how  this  approximation  is  obtained.  In  principle, 
one  could  also  take  the  first-order  approximation  in  which  the  drive 
is  represented  by  the  expansion  (which  does  not  increase  the  order 
of  the  system) 

Wn(p)*l -2l,,Tnp 
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Figure  6. 

The  root  hodograph 
of  aircraft  roots  with 
real  and  inertlaless 
drives  is  shown  in 
Figure  5*  The  transfer 
function,  generating  the 


Figure  7. 


"slow"  motion,  is  written  as 


K(P) 


7y+2{.7Yp-H 
fjp*+ W*+ p+ i+ke(1 


(3) 


The  "fast”  component  of  the  composite  motion  is  conveniently  taken 
in  the  form  of  the  first-order  approximation.  Then  the  drive  (the 
unit  generating  this  type  of  motion)  is  closed  through  the  object 
taken  in  the  form  of  an  integrating  unit.  The  root  hodographs  of  the 
drive  roots  in  real  and  idealized  schemes,  confirming  this  approxima¬ 
tion,  are  shown  in  Figure  6. 

The  formation  of  the  "drive"  transfer  function  is  shown  in  Figure 
7.  Loop  *4  does  the  "generation,"  and  loop  5  extracts  the  desired 
motion. 
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Analytically  this  has  the  following  form.  We  assume  that  the 
transfer  function  of  loop  4 

W  (a)  -  M7y+a8"r»*+0 

A  (i  j/>5+^iiifiiP+ 

consists  of  a  sum  of  the  desired  and  the  transfer  function  of  loop 
5,  which  is  the  image  of  the  "slow"  motion 

U.  <5) 


The  denominator  of  the  unknown  transfer  function  is  determined  by 
the  expansion  of  the  function 


+  (JW.  -rW,+ 

+  +  2  a,/r< 


(6) 


in  a  Laurent  series.  Equation  (6)  can  be  obtained  by  dividing  the 
denominator  of  the  transfer  function  W*(p)  by  the  denominator  of  W5(p). 


Neglecting  the  remainder  in  (6),  we  obtain  an  approximate  expres¬ 
sion  for  the  denominator  of  the  unknown  transfer  function.  An  analysis 
of  the  transient  error  processes  for  the  system  in  Figure  1  shows  that, 
taking  into  account  the  transfer  function  for  the  "slow"  component, 
Wfe(p)  should  have  the  form 

_ _  k,p _ _____ 


where  kt  *  2£.T.„,  is  obtained  using  the  method  of  indeterminate 
dr  dr 

coefficients  and  the  condition  that  the  original  W*(p)  is  approximately 
equal  to  the  sum  of  the  approximate  transfer  functions  in  (5)  and  (7). 
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Thus,  the  reaction  of  the  original  system  has  been  decomposed 
into  two  superimposed  components,  given  by  the  transfer  functions  in 
(3)  and  (7). 

4.  Next,  we  assume  that  both  components  of  the  composite  motion 
are  sufficiently  different  in  speeds  of  action.  This  assumption 
enables  us  to  direct  each  ohe  of  them  to  the  ’’lower"  and  "upper" 
filters,  respectively  of  the  spectral  analyzer. 

The  forming  filters  are  considered  as  follows: 

The  value  of  the  root  corresponding  to  Loop  1  in  Figure  1  is 
chosen  so  that  its  modulus  lies  between  the  roots  of  the  drive  and 
the  aircraft,  and  the  value  of  the  root  corresponding  to  Loop  2  is 
much  greater  than  the  modulus  of  the  drive  root. 

In  view  of  the  above,  in  the  case  of  a  "slow"  signal,  the 
forming  filters  1,  2  in  Figure  1  may  be  accounted  for  by  the  differ¬ 
ential  loop  Tp,  2  in  Figure  4.  In  fact,  a  "slow"  signal  is  a  sum  of 
a  stepwise  component  and  a  component  that  increases  from  0  [see 
Equation  (3)].  The  stepwise  signal  is  separated  at  the  input  to  the 
low-frequency  filter  2  in  Figure  4  in  the  form  of  a  short  (as  compared 
with  the  time  of  the  transient  process  in  the  low-frequency  filter) 
pulse,  whose  effect  on  filter  2  will  be  practically  identical  with 
the  effect  of  a  delta  function  (since  the  areas  of  the  delta  function 
and  the  real  pulse  are  identical  and  equal  to  T).  However,  the  "flat" 
part  of  the  signal,  in  view  of  its  small  slope,  will  practically  be 
differentiated  by  Loop  1  (Figure  1)  and  passed  without  distortions 
by  Loop  2. 

The  "fast"  signals,  due  to  its  high  "speed  of  action,"  will  be 
passed  by  either  filter  with  hardly  any  distortion. 

The  above  is  reflected  in  the  structural  schemes  in  Figures  4 
and  7,  showing  two  separate  channels  along  which  the  input  signal 
proceeds. 
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In  addition,  to  account  for  the  effect  of  the  stepwise  component 
of  the  "slow"  signal  [see  Equation  (3)]  on  the  high-frequency  filter 
7  in  Figure  7,  we  introduced  an  additional  single  connection  6. 

5.  To  make  an  analytic  calculation  of  the  possible  scheme,  the 
moduli  were  replaced  by  quadratures  3  in  Figure  4,  and  8  in  Figure  7. 

a  simulation  of  the  complete  Figure  1  and  the  separated  Figure  4 
and  7  showed  a  good  agreement  between  the  results  with  respect  to 
the  adjusted  value  of  the  coefficient  uadj. 

An  error  in  the  value  of  uadj  was  no  greater  than  10*  in  various 
regimes. 

It  is  important  that,  when  the  convolution  integral  is  applied 
(in  a  complex  domain),  the  order  of  the  original  integrand  in  the 
exact  scheme  (Figure  1)  is  nine.  When  the  procedure  proposed  here  is 
used,  the  order  is  lowered  to  four. 

6.  The  subsequent  design  of  the  system  can  be  made  by  finding 
the  quadratic  estimate  at  the  output  from  each  of  the  channels.  This 
permits  us  to  obtain  an  analytic  expression  which  relates  the  basic 
parameters  of  the  system. 
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AN  INVESTIGATION  OP  ADAPTIVE  CONTROL  SYSTEMS 
FOR  RANDOM  ACTIONS 

M.  P.  Rosin  and  V.  I.  Ul’yanov 

At  the  present  time,  adaptive  control  systems  are  studied 
statistically  by  using  the  interpolation  method  [3],  the  method  of 
equivalent  perturbations  [4],  and  the  method  of  coefficients  [1]. 

Among  the  methods  in  which  the  output  signal  is  represented  by  poly- 
nomians,  the  interpolation  method  has  the  highest  accuracy.  In  this 
method,  the  interpolation  points  are  taken  to  coincide  with  the  roots 
of  the  orthogonal  polynomials,  with  the  weights  equal  to  the  coeffi¬ 
cient  distribution  densities  in  the  canonical  expansion  of  the  input 
signal  to  the  control  system. 

The  accuracy  of  the  probabilistic  characteristics  obtained  by 
the  interpolation  method  is  checked  by  comparing  the  results  for 
various  degrees  q  of  the  approximating  polynomial.  The  entire  compu¬ 
tational  procedure  is  repeated  again  for  a  new  value  of  q,  since  the 
optimal  interpolation  points  change.  In  contrast  with  the  interpola¬ 
tion  method,  the  method  of  coefficients  uses  the  results  of  previous 
simulations.  This  lowers  substantially  the  amount  of  work  involved 
in  the  calculation  of  the  probabilistic  characteristics  of  the  adaptive 
control  system.  For  example,  for  q  *  4  and  p  »  7  in  the  interpolation 
method  and  the  method  of  coefficients,  one  has  to  carry  out  2187  and 
113  simulations  of  the  control  system,  respectively. 


Ill 


However,  due  to  the  nonoptimal  selection  of  the  approximating 
points,  the  method  of  coefficients  is  inferior  in  accuracy  to  the 
interpolation  method.  One  can  increase  the  accuracy  of  the  coefficient 
method  by  taking  into  account  different  effects  of  the  coefficients  in 
the  canonical  expansion  of  the  input  signal,  and  by  combining  the 
coefficient  method  with  the  Interpolation  method. 

Determination  of  the  Probabilisitic  Characteristics  of  a  Control 
System  by  Taking  into  Account  Various  Effects  that  the 
Coefficients  in  the  Canonical  Expansion  of  the 
Input  Signal  Have  on  the  Output  Signal 

The  selection  of  the  degree  q  in  each  random  coefficient  in  the 
coefficient  method  is  achieved  by  improving  the  accuracy  of  the  proba¬ 
bilistic  characteristics  of  the  control  system,  and  does  not  require 
any  additional  simulations  of  the  control  system.  Let  us  select 
Chebyshev-type  points  in  such  a  way  that  their  number  will  be  suffi¬ 
cient  to  determine  the  coefficients  in  the  expansion  of  the  output 
signal  In  a  Maclaurin  series.  These  points  satisfy  the  requirements 
cf  the  coefficient  method:  they  lie  within  the  ranges  of  the  real 
values  of  the  random  coefficients  in  the  canonical  expansion  of  the 
input  signal,  and  the  output  signals  for  them  are  not  very  small. 

Let  us  consider  the  procedure  in  which  the  variance  of  the  output 
signal  from  an  adaptive  control  system  is  calculated  taking  into 
account  different  effects  of  the  random  expansion  coefficients. 

We  have 


O,(0-M  <T(/,  Xu  X, . X,,)) 


(1) 


Here 


T  (t,  Xt,  X, . A'rHK(/,  .V,  A1 . X„)  (01* 


(2) 
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X»,  X2,  ....  Xp  arc  the  random  coefficients  of  the  canonical  expansion 
of  the  input  signal  with  the  variances 


M  (*;M,  M  [All-o'i,  ...,M  tvjl-flg. 


(3) 


The  expansion  coefficients  will  be  assumed  to  be  distributed 
normally.  The  degrees  q2,  q2,  . ..,  qp  of  the  expansion  coefficients 
X,,  X2,  . ..,  Xp  in  the  polynomial  approximating  the  function 
4*(t,  Xi ,  X2,  ...»  Xp)  are  unknown.  To  a  first  approximation,  we 
assume 


According  to  the  coefficient  method,  we  find  the  first  approximation 
to  the  variance  of  the  output  signal  from  the  control  system 


Here 


D°y  (Q-¥(l,  0, 0 . 0)+2  l  |A*  (t,  **+A, 

*«"<  ‘  -V* 


V  it,  0,0,.. .,  0)—Y  (/,  xlt  Xi,...  xp)=*A0  (0 
for  . . .  *~Xp**0 


(5) 


(6) 


M  *•)-**</.  *,)-V(*.0,  0 . 0) 

(*•**)“¥(*.. *1.  x-i . xp) 


(7) 


for  .*,=>*2= 


’•**-» -.**+1* 


=xps0  and  xkr-0 


(8) 


Chebysher-type  points  are  chosen  as  the  concrete  values  of  the  random 
coefficients  [33. 
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0 
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2P 
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0 
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0 
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We  increase  the  degree  q  in  the  expansion  of  the  function  f(t,  xj, 

V  V  ln  a  Maclaurin  series  in  only  one  ("large")  random  coe:- 
ficient  xy,  so  that  the  highest  degree  of  the  approximating  polynomial 
of  the  function  ¥(t,  xit  x*,  . ..,  xy,  xp)  does  not  exceed  four. 

The  algorithm  in  which  the  random  coefficients  are  divided  into  "large" 
and  "smali"  is  given  in  Section  2. 


(A  . . .  •  •  -A>,r  A*i*'  fv)  \S:a,  (t, 

*»•  V 

.  ft  ft 

i'  /  xJ'Xl  •  \\ 


(9) 


I,  /Vi,  /-V 


'W*  A0-4>(*)M  ov(0‘*v+j  ^v(<)-4+ 

r,^v(0^+-1/i4(o4  (io) 


A*  (<’  x*)  ~Akl  (t)+ A\  (/).Xv+L  a4  (t)  4+1  a4  (<)  4 


(ID 


A,»,(*.  (f)jfv+I  yl 4  (t)xl 

*  Mi 


(12) 
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Substituting  Equations  (10),  (11),  (12)  in  Equation  (9),  we  have 


¥(/,*„  . . .  *v, . . .  *pHAiM+4v  W*+j  a4(().xI+ 

p 

+j,^(/)4+i^(<)4+  2  4h  (<)+<*(/)*+ 


*t1*V 


+y^(0^+i^»W4l**.+fS  2  Kh(<)+ 


*i=*i  *i  ->i 
*iv»vtn*v 


(13) 


+  K%  (*)  ^v+jA^  (t)  xJl  XMl 


Applying  the  expectation  operator  to  Equation  (13 ),  we  obtain 


D\  (0-0$  (/)+ AD*v  (t)amAo  (/)+ lyl4  (t).M  {4)+Jj  a4«)'M  Uu- 

P 

+l2A4(i)M{x*)+l  2  (K„(/)+j^v  (0Ai{4>IA/{4,) 

•  kx**l  *•*» 


(1*0 


The  coefficients  k  (t)  are  determined  from  the  known  formulas  by 

Z 

using  the  results  of  the  previous  simulations  of  the  control  system. 
The  coefficients  Aw  (*)  are  calculated  using  the  formulas 


*v*  **i)+Av,*,(A— Xv— '**.)+■ 

-Xv-Xj-(-AVl*t(/,  Xv-X4l)I-j[Av(h  A’v)  I-  (15) 

+K(t,  -~.Vv)]---  K(/.  *J+A,,  (/,  -.vj| 


The  number  of  simulations  of  the  control  system  for  one  "large"  random 
coefficient  as  compared  with  the  first  approximation  to  the  variance 
D®(t)  is  increased  by 


A  4  (p — 1). 


(16) 


Since  the  Chebyshev-type  optimal  points  were  chosen  as  the 
approximation  points,  some  of  the  terms  in  the  variance  (1*1)  can  be 
calculated  using  the  Interpolation  method 


115 


(17) 


M  <'F  (t,  0, 0 . .tv . 0)-4>  W+7  a}  (i)  M  W  i* 

+ J-,  A'Z(t)  At (4)-'F(/, 0, . . ..0)ft,+V(/.0 . 1,732 . 0jp,+ 

+Y(/, 0, . . . , —*  1 ,732 flvi 

Here  p0,  Pi,  Pa  are  the  Christoff el  numbers  used  in  the  normal  distri¬ 
bution  of  the  random  coefficient  with  the  following  values  [3] 


p0  **0,666;  p,-0,166;  0,166. 


(18) 


A  partial  application  of  the  interpolation  method  permits  us  to  reduce 
the  number  of  simulations  of  the  control  system  by  2»K^;  is  the 
number  of  "large”  coefficients  in  the  expansion.  If  the  increase  in 
the  variance  of  the  control  system  AD*v(t),  caused  by  an  increase  of 
the  degree  q^  relative  to  the  "large"  random  coefficient  xy,  does  not 
exceed  K*D®(t),  then  in  the  case  of  the  random  coefficient  xy,  we 
limit  ourselves  to  the  obtained  value  of  qy. 


The  coefficient  K  characterizes  the  accuracy  to  which  the  variance 
is  determined.  It  can  be,  for  example,  taken  as  equal  to  K  »  0.05. 

If  ADy\>(t)  >  KD®(t),  then  one  should  consider  —  in  the  expansion  of 
the  output  signal  —  terms  of  a  higher  degree  q^  in  the  random  coef¬ 
ficient  xy.  Upon  finding  AD*v(t),  we  proceed  to  consider  two  "large" 
random  variables  xy  and  x^.  For  this  case,  in  the  formula  for  the 
variance  of  the  output  signal 

O;  (t)~Dl  (/)+**£•  W+ aoJh  (/)  (19) 

we  add  the  terms 


Ti  (0  M  {4),  {■  Afa  (/)  M  {xl)  M  {.vi,) 


whose  coefficients  are  determined  using  the  formulas  similar  tc 
Equations  (7)  and  (8). 
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The  number  of  additional  simulations  of  the  control  system  in 

order  to  calculate  the  variance  D*(t)  as  compared  with  the  variance 

Dl(t)  is  increased  by 

y 

A'L,  ~4(p-2)  (20) 

The  increase  in  the  variance  AD*y(t)  due  to  the  random  coefficient 
X  is  compared  with  the  value  of  K*D®(t).  If  AD*p(t)  <  KD°(t),  we 

r  v  v  v 

limit  ourselves  to  the  value  s  otherwise  the  degree  y^  in  the 
random  coefficient  is  increased  to  six. 

The  increases  AD  (t)  and  degrees  q  needed  to  take  into  account 
the  successive  random  coefficients  are  calculated  in  the  same  order. 
For  p  *  5,  p  =  10  and  q  **  4,  we  give  below  the  necessary  number  of 
simulations  of  the  control  system. 

Number  of  simulations  using  the 
coefficient  method  considering 
different  effects  of  the  expan¬ 
sion  coefficients 


ia 

I 

2 

3 

4 

5 

G 

7 

8  ' 

9 

- - 

10 

p-5 

27 

39 

-17 

51 

51 

! 

G1 

p— 10 

57 

89 

117 

141 

1C1 

177 

189 

197 

201 

201 

1 

221 

Thus,  the  number  of  simulations  of  the  control  system,  taking  into 
account  different  effects  of  the  random  coefficients  in  the  canonical 
expansion  of  the  input  signal,  is  much  lower  than  the  number  of  simu¬ 
lations  as  compared  with  the  usual  coefficient  method;  or,  for  the 
same  number  of  simulations,  the  accuracy  of  the  coefficient  method 
taking  into  account  different  effects  of  the  coefficients  is  higher 
than  the  accuracy  of  the  ordinary  coefficient  method. 

2.  The  algorithm  for  dividing  the  random  coefficients  into 
"large"  and  "small." 


Coefficient  method 
without  considering 
different  effects 
of  the  expansion 
coefficients 
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Let  us  consider  the  algorithm  in  which  two  coefficients  in  the 
canonical  expansion  Xj  and  X2,  normally  distributed,  are  divided  into 
"large"  and  "small." 


Using  the  coefficient  method  for  Chebyshev-type  points 


0,000  0,000 

+  1,732  0,  +  1,732  0, 

-  1,732  o,  -  1,732  a, 


we  find  the  components  of  the  first  approximation  to  the  variance  of 
the  output  signal. 


M  <M'o (*,  0)j  -  V(t,  0.  0)  +  ^  Aj(t)-M  {xj} 


(21) 


0.  x,)}-^,  0,  Ol+jAjW.Ai;,-). 


(22) 


"Using  the  same  results  o^  a  simulation  of  the  control  system, 
we  use  the  interpolation  method  to  find  the  expectation  of  the 
functions  i|<l(t,  xj,  0),  <i*(t,  o,  xa)  for  q*  «  qa  *  4 

7W  {'!» (/,  x, .  0)}  0,  0)+i  AK  M  {xi}  +  i*, (/)M  {*!}  -  ( 2 

-  ¥(/,  0,  0)?u  +  V(f-  1,732 o„  0) p,  +  V(i,  -l,732o„  0)p, 


{V  {/,  0,  .v,)}  ~  Y(/.  0,  0)  +  i  +  ^  /i;,(/)/W{A.ij-: 

-V(7.  0,0)p0  +  Y(/,0,  l,732o,)p,  +  V(t,  0,  - 1,732  o.)Pi 


We  calculate  the  differences 


M  (/.  x„  0))  -  M  < '!«>(/,  x„  0)) « 

At  l’1"  <f'  0.  -  M  {'F(y,  0,  x,)}  ~ax, 


(24) 


(25) 

(26) 
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The  coefficients  aXj,  characterize  the  effect  of  the  fourth-deeree 
terms 

A*  {■*!}• 


on  the  value  of  the  variance  of  the  output  signal  from  the  control 
system. 

The  random  coefficient  with  the  largest  value  of  ax  is  taken  as 
the  "large"  random  coefficient  when  going  from  the  first  approximation 
D^(t)  to  Dy(t).  The  transition  from  "small"  to  "large"  random 
variables  is  done  in  the  descending  order  of  a 

xv 

The  transition  from  qv  ■  4  and  qv  *  6  is  done  similarly. 
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EQUATIONS  OP  MOTION  OP  A  SYSTEM  OP  BODIES 
OP  VARIABLE  MASS  AS  A  CONTROL  OBJECT 

Zh.  S.  Agayev,  B.  V.  Viktorov,  and  I.  S.  Ukolov 

A  majority  of  objects  in  aircraft,  missile,  and  Bpace  technology 
are  mechanical  systems  that  can  be  represented.  If  one  neglects  the 
elasticity  of  their  construction,  as  ensembles  of  a  number  of  Ideally 
rigid  bodies  of  variable  mass  whose  relative  motion  and  the  variation 
of  mass  are  assumed  to  be  known.*)  Such  objects  Include,  for  example, 
spacecraft  with  cosmonauts  or  other  movable  mechanical  masses  moving 
in  them,  aircraft  of  variable  geometry,  and  a  number  of  other  objects. 

Due  to  the  shifting  of  the  centers  of  mass  of  the  component  bodies 
the  gyroscopic  linkage  between  rotating  bodies,  the  shifting  of  the 
axes  of  Inertia,  the  change  of  the  moments  of  Inertia  of  the  ensemble 
In  its  entirety,  etc.,  the  dynamics  of  the  above  mechanical  system 
differ  substantially  from  the  dynamics  of  an  Isolated  rigid  body  of 
variable  mass.  Numerically  this  difference  is  greater  if  the  mass  of 
the  moving  bodies  is  greater  compared  to  the  mass  of  the  main  body, 
and  If  the  movements  of  the  masses  are  faster. 


*The  term  has  the  following  meaning:  the  relative  kinematics  of 
bodies  and  the  character  of  mass  variation  of  each  one  of  them  are 
known  functions  of  time  or,  in  a  more  complicated  case,  are  the 
phase  coordinates  of  another  more  complicated  dynamic  system. 
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It  is  important  to  note  that,  for  the  particular  case  of  bodies 
of  constant  mass,  the  equations  of  this  type  (only  for  rotational 
motion)  were  obtained  in  [1]. 

The  mechanical  system  considered  here  has  of  course  six  degrees 
of  freedom.  Therefore,  it  is  necessary  and  sufficient  to  set  up  six 
independent  dynamic  equations  of  first  order,  with  the  variables  in 
the  form  of  the  components  of  the  linear  velocity  vector  of  an  arbi¬ 
trary  instantaneous  center  of  mass  and. the  components  of  the  angular 
velocity  relative  to  that  center. 

In  the  present  paper,  the  center  of  reduction  is  the  instantaneous 
center  of  inertia  of  the  system.  The  scalar  equations  for  transla¬ 
tional  and  rotational  motion  are  written  relative  to  the  instantaneous 
central  axes. 

The  hypothesis  of  near-action  [2]  is  used  as  the  initial  assump¬ 
tion,  as  is  done  in  a  majority  of  papers  on  the  mechanics  of  variable 
masses. 


The  Coordinate  System  and  Kinematic  Relations 

Let  us  consider  the  case  when  the  system  consists  of  the  main 
body  "0"  and  a  number  of  movable  bodies  "i"  (Figure  1).  We  shall  use 
the  following  coordinate  systems: 

CnS  —  inertial} 

XoY«Z#  —  attached  to  the  main  body  ”0''; 

XiYiZi  —  attached  to  body  "i"} 

XYZ  —  movable  system  whose  origin  coincides  with  the 
instantaneous  center  of  Inertia. 

Since  the  relative  motions  of  bodies  and  the  change  in  their  mass 
characteristics  are  considered  to  be  given,  the  following  quantities 
are  also  assumed  to  be  given:  Pi  which  is  the  radius  vector  giving 
the  position  of  the  origin  of  the  XYZ  system  relative  to  XoYoZ#;  the 
Euler  angles  giving  the  relative  orientation  of  the  XYZ  and  XoYoZo 
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axes;  is  the  angular  velocity 
vector  of  the  system  X^Y^^  rela¬ 
tive  to  XaY«Za;  p#i  is  the  radius 
vector  giving  the  origin  of  the 
system  X^Y^  (i  *  1,  2,  ...» 
n)  relative  to  XoY«Zi;  the  Euler 
angles  giving  the  relative 
orientation  of  the  XjY^Z^  ana 
X|Y|Zt  axes;  ft*  is  the  angular 
velocity  vector  of  the  system 
XYZ  relative  to  X#Y#Zc*;  n1  is 
the  angular  velocity  vector  of 
the  system  XYZ  relative  to 

*iYiZi*  moreover 

Oo  —  iuj  +  U,  (1.1) 

Figure  1. 

are  the  Euler  angles  coordinating  the  orientation  of  XjY^Zj  and  XoY0Z». 
Of  course,  a  knowledge  of  the  orientation  of  XYZ  and  X»Y0Za,  X^Z, 
and  XYZ  uniquely  determines  the  orientation  of  X«Y*Z»  and  XjYjZj. 

p;  is  the  radius  vector  giving  the  location  of  the  origin  of  XYZ 
relative  to  X^YjZ^  (i  ■  1,  2  ...  n)j  is  the  radius  vector  giving 
the  position  of  the  instantaneous  center  of  inertia  "0M  relative  to 
XYZ;  Is  the  radius  vector  giving  the  position  of  the  instantaneous 
center  of  inertia  of  a  body  "i"  relative  to  XYZ;  p#c  is  the  radius 
vector  giving  the  position  of  the  instantaneous  center  of  inertia  of 
the  body  ”0"  relative  to  XaY|Z«;  pic  is  the  radius  vector  giving  the 
position  of  the  instantaneous  center  of  inertia  of  a  body  "i"  relative 
to  X,Y,Z,;  in  and  m«  are  the  instantaneous  masses  of  the  body  "0"  and 
”i,"  where  by  assumption 


lm  I 


0 


(1.2) 


•Translator’s  note:  Incorrectly  given  as  ,!Y"  in  foreign  text. 
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and 


n 


m  =  m0  -|-  2  in,; 

tr-l 


( 


*» 

X  • 


3) 


mj0  and  are  instantaneous  masses  of  "current  points  J"  forming 
the  bodies  "0"  and  "i,"  respectively,  where 


:  m,.;  2  nt 
I 


■n 


a 

m,;  22  nt,,  =*  in. 

1*0  I 


(1.4) 


The  above  points  are  referred  to  XYZ  using  the  radius  vectors 
and  respectively,  where 

Ay  i  *»  Ay  (I  +  y j,d  -i-  tj,n  (/•»!...«)  (1.5) 

IX\'i*  ^yvi»  I2Zi*  IXYi*  IYZi*  ^ZXi  are  the  in8fcantane0UB  moments  of 
inertia  of  a  body  "i"  in  the  central  axes  XYZ.  Ixx,  IYY>  I„z,  I  , 

Iv„,  I.,v  are  the  instantaneous  moments  of  inertia  of  the  ensemble 
relative  to  the  central  axes  XYZ. 

q0  is  the  local  derivative  of  p»  in  the  X0YoZc  system,  and  it 
characterizes  the  velocity  of  the  origin  of  XYZ  relative  to  the  boun¬ 
daries  of  the  body  "0,",  qi  is  the  local  derivative  of  pi  in  the 

X.Y.Z.  system,  and  it  characterizes  the  velocity  of  the  origin  of  XYZ 
111  § 

relative  to  the  boundaries  of  a  body  "i";  qp  is  the  local  derivative 

* 

of  q#  in  the  X«YoZp  system;  qA  is  the  local  derivative  of  qi  in  the 
XlYiZl  system»  9*1  is  ioca1  derivative  of  p o in  the  XpYoZo  system; 
q0i  is  the  local  derivative  of  q#^  in  the  XYZ  system;  RQ  is  the  radius 
vector  giving  the  origin  of  the  XYZ  system  relative  to  cn£;  V  =  d/dtRQ 
is  the  derivative  of  RQ  in  the  cnC  system;  d/dt  V  is  the  derivative 
of  V  in  the  cn£  system;  Z,  d,  n  are  the  unit  vectors  along  the  XYZ 
axes.  From  the  vector  triangle  in  Figure  1,  we  can  write 


fo  “  Pol  +  Pi 


(1.6) 
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Then  the  local  derivative  of  pg  in  the  XgYgZg  system  is 


tfo  “  9oi  +  +  rt»Xpi 


(1.7; 


and 


q0  *=  Mo/  + “/  X  Pi  +  “i  X  (“i  X  Pi)  +2'>,  '<  Qi  -r  Hi> 


(l.B) 


« 

where  is  the  local  derivative  of  in  the  XgYgZ#  system  equa  ’  '  '*■ 
the  local  derivative  of-  toi  in  the  X^Y^Z^  system.  Clearly, 


didtV  =  V  +  K  +  U )  x  V 


(1.9) 


and 


V  —  V xi  +  Vyd  -1-  Vztt, 


(1.10) 


where  V  is  the  local  derivative  of  V  in  the  XYZ  system. 

The  total  angular  velocity  of  a  body  wi"  and  its  angular  accel¬ 
eration  are  equal 

“in  «■  (1.11) 


dldt w(„  ^  dldt  (o0  +  dldt u>(, 


(1.12) 


respectively,  where  uo  is  the  angular  velocity  vector  of  the  X0YoZo 
system; 

* 

dldt  “  u»i  +  %  X  H>c 

(1.13) 


Furthermore, 


dldt  *■*  •«,,  +  (w„  Ql)  ,  <  iii^, 


(1.1*0 
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where  w0  is  the  local  derivative  of  uq  in  the  XYZ  system.  Finally 


(licit  %  =  «»„  -f  Q0  X  «v 


(1.15) 


It  will  be  noted  that 


U)„  »■*  ‘#q.\7  •!*  +  ‘‘’(i/W* 


(1.16) 


where  wqX,  oj0Y,  w0Z  are  the  projections  of  u>q  on  the  XYZ  axes.  By 
analogy,  we  can  also  write 


(l /(It  <0,*  =*  (ily  +  X  %  ‘i-  «»<  T  K  -I"  -()  X  »/. 


(1.17) 


where  is  the  local  derivative  of  in  the  XYZ  system. 


The  Force  Equation 


In  accordance  with  Equation  (22)  from  [3]  and  the  notation  used 
here,  the  vector  equation  of  the  translational  motion  of  a  body  of 
variable  mass,  written  for  the  body  "0,"  has  the  form 


muti:dtV  «  F„  -i-  R,,  h  m0  ( -  <1  'ill  >■>„  x  %  ■ 

■r  +  </«  ’I*  H  K«J. 


(2.1) 


where  Fq  is  the  resultant  of  all  the  extex’nal  forces  acting  on  the 
body  "0."  This  may  include  forces  which  are  internal  relative  to  ’.he 
ensemble.  Ro  is  the  resultant  of  all  the  external  (with  respect  tc 
the  body  "0")  reaction  forces.  It  may  also  include  forces  that  are 
internal  relative  to  the  ensemble. 


The  structure  of  the  expression  in  (2.1)  remains  the  same  for  a 
body  " i . "  One  must  only  change  the  indices  and  write  out  the  to’.al 
angular  acceleration  and  velocity  vectors  in  more  detail. 


In  view  of  (1.13),  (1.14),  and  (1.15)  after  a  number  of  trans¬ 
formations,  we  obtain 


m,d/dt  V  —  F,  +  R,  +  mt  \-dldt^X^ci  —  K  A/H  +  u>oA,»i  + 
+  m,  |— «|XAf<  —  («tlArt)ui|  +w<At,J  +  2mi(uVo,)Aw  - 
-  2m,(mQ-bcl)<a,  +  m,(qt  +2u>0xq,  -I-  2«>,Xq<). 


Adding  Equations  (2.1)  and  (2.2),  we  obtain  in  view  of  (1.2)  and 

(1.3) 


where 


md/dtV  «*  F  +  R  +  2  mt  f— «/XAri  —  («><Af()  ml  +  uijAf(|  + 


/«l 


+  2  2  ("V'-/)  —  2  2  mt  (“V  A,/)  + 


/-i 


<~i 


+  |’nu4i  +  2  j  +  2u»o  X  |«c<lo  +  2  ntiQi'j  + 


+  2  2  mi  (wj  X  di)i 

/r.| 


(2.3) 


F-Ft  +  %F,  and  /?  -  +  2 /?,. 

/-i  (-1 

It  should  be  noted  that  F  and  R  do  not  now  oontain  any  components 
internal  to  the  ensemble. 

Equation  (2.3)  is  a  vector  equation  determining  the  motion  of 
the  instantaneous  center  of  inertia  of  the  system.  If  another  point 
is  taken  as  the  center  of  reduction,  the  structure  of  the  right-hand 
side  of  (2.3)  will  turn  out  to  be  mu'.h  more  complicated. 
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Let  us  project  the  vector  Equation  (2.3)  onto  the  axes  of  the 
moving  system  XYZ.  Here  in  order  to  simplify  the  structure  of  the 
scalar  equations,  we  shall  represent  only  the  principle  terms  from 
(2.3)  in  the  projection  onto  the  XYZ  axes.  They  characterize  the 
dynamics  of  the  translational  motion  of  an  instantaneously  "fixed" 
system.  The  remaining  terms  in  an  expanded  form  will  be  written 
separately,  outside  of  the  equations. 

Thus,  denoting  the  additional  terms  in  (2.3)  by  a  vector  Qi  and 
taking  Equations  (1.9)  and  (1.10)  into  account,  we  find 

m  \Vx  +  Ky  **“ “<jy) V z  -  (w02  +  2o i)  Vv\  *=‘  Fx  +  %x 
>n !^y+Hz+2«) Vx  -  Kc+SW/,1  -  Fy+R.,-\  Q.y,  < 2  •  “ ) 

*  l^+Kx+^x)^  -(v+ V^xi  “ 

where  w0X»  <*>oY»  «#Z»  8#X»  fi«Y»  A#Z  are  the  components  of  wt  and  in 
the  XYZ  system;  Fx,  Fy,  Fz,  flx,  R^»  Rz,  Qxx>  Qiy,  Qiz  are  the  compon¬ 
ents  of  F,  R,  and  Qi,  respectively,  in  the  XYZ  system. 

Let  us  write  the  expanded  scalar  expressions  for  the  vector  Qi 
in  the  form  of  separate  groups  of  terms. 

The  first  term  is 

« 

/  r\  r-  \ 


In  view  of 


(2.6) 

Equation  (2.6),  considering  (2.1),  becomes 


^  —  <0,  +  8,x>«< 


(2.7) 


Substituting  (2.7)  and  (2.5)*  and  considering  that 
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bel-x(ll+yeld+2trn, 


we  get 


2  »h ((-«,y *ci+ 

+ (“>«  -  •«  *aw+ V  fi/x  -  "/A-  •  ‘  V.  '  yw + (“?y+“?«)  *wi  /+ 

It 

+  2  mj ((—  i/Jt  —  •,r*a«+-«*fln--  V-«)^  + 

/“i 

+("«  ~‘“/y'a»2+u,«a/y  ~  Vu,«)^+(u’«+"iA)yrtlrf+ 

/I 

+  2  *i  I(-«Vr -»,* ■a/y+u>/y'a/z-“:z-V  Vd+ 

i  *1 

+(lV  —  u>/2,y«+u'ix  ,aiz  -  */z  •*/*)  J{«i+(®«f +**y)  *wl  «• 


The  second  term  is 

"  a 

2  2  mi  (uVu,i)  ^ci  *~  2  2  m,  (“b* At/)o»4 «« 

(-1  Ul 

H 

"  2  2  IKy ,<#<y +%4 - \Y"*txy<i - %t ••«*dl  *+ 

M 

+  3  2  *»i  K*  ••«+<**••«)  y£(  ~  ~ u>oJr'“/y,JCdl  d+ 

II 

+ 2  2  «i  lKy**<y+“«x  'vix)zci  —  tt,ox"#/z  ~ 

i-i 

— ®oy ‘®  n* 

The  third  term  is 

»wo+2  «'?<• 

W 


(2.8) 


(2.9) 


In  view  of  (1.1)  and  the  scalar  expressions  for  the  local  derivatives 
of  qo  and  qA  in  XYZ 

?0  ■“/0Xl+1tY<l+ioZn> 

<hm-i,xl+qiv(l+qixn, 


we  obtain  the  following  scalar  expression  for  the  third  term 


{"'o  (Qox  +  ^02  '  2oy  (1cy'~oz)  + 


ft 

+  2  "it  (V«  +  (-oy  -  <VHz  -  (Qoz-(,);z)'?/y]}l  + 

tel 

(/tig  (^Oy'^’^OX  '  ^0  Z  Qoz  ‘  ~0.y)*^ 

» 

+  2  foy  +  (S0Z  -  w/z)  ?/*  -  (20r  -  “»)  ^*l>d+ 

let 


+  {/»o  (?oz  ?ov  ‘  2ox  ~  ?ox  ‘  2oy) 

n 

+  2  «<  (V,z  +  (Qox  “  “/x)  <liY  ~  (s0 y  -  “ly)  <?/xl>  «• 

tel 


(2.10) 


The  fourth  term  is 


(It  /I 

"Wo  +  2  Mtfj)  "■  2u)0  X  2  “• 

/-I  1-0 

n 

1 2  S  «|(«0 Y-lit-Oot-liv)1  +  2 2  »i(“oz'9«  -  *0 x'1tz)d  + 

jsO 

A 

+  2  S  »i  Kx  •  ^ly  ~  “oy  ‘  ?<x)  «• 


(2.11) 


/=-o 


The  fifth  term  is 


A  A 

2  2  «<(<»*  X  id  -  2 s  *i(»ir  *fa - ^z-9/y) 1  + 

/•^t  /Ml 

A 

+  2  S  /«,  (mu  -qlx  -  y -?,2)  d  + 

A 

+  2  2  m, (\x-qiY  -  «/K •  qlx)  n. 


(2.12) 


Using  Equations  (2.8),  (2.9)»  (2.10),  (2.11)  and  (2.12),  one  can 
make  a  numerical  estimate  of  each  additional  term  in  Equation  (2.3). 
In  addition,  these  expressions  can  be  combined  with  the  corresponding 
equations  in  (2.4)  in  order  to  obtain  a  more  exact  solution  for  V^, 
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The  Moment  Equation 


According  to  Equation  (42 )  in  [33  and  the  notation  used  here, 
the  vector  equation  of  the  rotational  motion  of  a  body  of  variable 
mass,  set  up  for  the  body  "0W  relative  to  the  origin  of  XYZ,  has  the 
form 

2  faifao  X  [d/dt<x ip  X  Ay0  -f  fafaj  u>0  —  w£Ay0)  —  Alf,  +  M#,  + 

+  no  [dldtV  X  ^  +  Afo  x  q0+  2A(0  X  fa  X  qJJ,  ( 3  * 1  * 


where  Mp#  is  the  moment  of  all  external  (relative  to  M0")  forces; 

M_  is  the  moment  of  all  external  (relative  to  "0")  reaction  forces, 
n# 

Mp#  and  Mr#  may  include  components  which  are  internal  relative  to  the 
ensemble. 


I 


i 


The  structure  of  (3.1)  remains  the  same  for  a  body  "i."  One 
must  merely  change  the  indices  and  write  the  total  angular  accelera¬ 
tion  and  velocity  vectors  in  more  detail  [see  (1.11),  (1.12),  and 
(1.13)].  After  a  number  of  transformations,  the  torque  equation  for 
a  body  "i"  will  become 

2  fafa/  | d!dt%  X  &ji  -f-  —  wjjAyjJ  »  Wfa  M«(  + 

+  \d:diV  X  Ki  -I-  A,/  X  q,  •!•  2Ae.  X  K  X  qi)J  •»-  2 X 

X  fa  X  q()  -  2  2  fa/  W^j,)  (Ay,  X  «»,)  -  ( 3  •  2 ) 

J 

—  ^  fa, Ay,  X  fa  X  Ay,  +  fa,  Ay,)  U>,  —  W’Ay,). 


Adding  Equations  (3.1)  and  (3*2),  we  obtain  in  view  of  (1.2) 

2  fafai  X  {d/dtUQ  X  Ay;  -j-  faAy,)  -  U/jjAy,)  a.  /Vfa  +  MR  + 

*  t  «  , 

/To  !  tl  ^  ^  **  X  (‘"0  Xq,)  + 2  2  faAf,  X  fa  X  q,)  — 


*  » 

-  2  2  2  fa/  faoAy,)  (Ay,  X  «>/)  “22  fa/ Ay,  X  fa,  X'  Ay;  -(- 

M  J  /»!  ! 

+  faAy,)W;-lUjAy,|t 


(3-3) 
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The  vectors  Mp  and  MR  do  not 


where  2  =  Mr~j\  M.  and  M*  =  j)  M»  • 

<=•■  0  /  HO  ‘  HO  ' 

contain  components  internal  to  the  ensemble. 

Equation  (3*3)  is  a  vector  relation  determining  the  rotation  of 
a  mechanical  system  about  the  center  of  reduction,  in  this  case  the 
<  instantaneous  center  of  inertia.  Just  as  in  the  case  of  Equation  (2.3)* 

.  the  structure  of  the  right-hand  side  of  Equation  (3-3)  will  be  more 

complicated  if  another  point  is  taken  as  the  center  of  reduction.  In 
j  particular,  it  is  precisely  for  this  reason  that  Equation  (3*3)  does 

j  not  contain  parameters  of  the  translational  motion. 

i 

j  Thus,  if  Mp  and  MR  are  functions  of  time  and  parameters  of  rota- 

|  tional  motion  alone.  Equation  (3.3)  nay  be  integrated  independently 

;  of  (2.3).  Thus,  the  total  system  of  equations,  composed  of  Equations 

'  (2.3)  and  (3-3)  and  determining  the  total  motion  of  the  mechanical 

system,  in  view  of  the  remarks  with  respect  to  the  functional  depen- 

i 

dence  of  F,  R,  Mp  and  MR,  is  connected  in  only  one  direction.  It 
;  should  be  noted  that  Equations  (2.3)  and  (3*3)  are  not  total  time 

i  derivatives  of  the  momentum  and  angular  momentum  vectors,  respectively. 

This  is  the  peculiarity  of  the  mechanics  of  bodies  of  variable  mass. 

Let  us  project  the  vector  Equation  (3*3)  onto  the  axes  of  the 
moving  system  XYZ.  In  so  doing,  we  shall  give  the  components  In  the 
XYZ  system  of  only  the  basic  terras  in  (3.3).  They  characterize  the 
dynamics  of  the  angular  motion  of  an  instantaneously  "fixed"  system. 

i xx  (“ox  +  ~oyo>oz  —  S —  (/yy  —  Izz)  —  hr  (l'W  — 

“Vll/C  +  Su^uio.x  —  ii[>Xwvz)  —  l xx  ((U0<J  +  '"UX-'W  ~ 

-  &>rw«x)  +  hz  K*  -  u>lr)  =  Mfx  +  A/«a.  +  Q.X, 

1  hr  (“o  y  +  2«“ux  —  2o*woz)  —■  Vzz  —  hx )  “W^x  — 

—  I YZ  (wo z  —  +  2oX“0 y  —  Goj'iOjx)  —  lyX  (ui0X + W,>)'Woz  +  (  3. 1|  ) 

+  Qqyu>oz  —  2o zmy)  +  Izx  (“;;*  —  w-^)  ™  MFy  -j-  Mry  -f.  Q„r, 

/ zz  (“o z  4*  2ox“oj'  —  Sci'Wox)  -  (hx  ~~  hr)  w«iuor  — 

—  hx  (luox  —  +  Qitf'Wox  —  Qoz^y)  —  h.r  (>»o y  ~r  *“oz“ox  4* 

+  2oz«X),v  —  2ox“o z)  +  hr  (%V  ~  “J.v)  “  4^''z  +  Mp.z  -ir  Q:z, 
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I  wawwnww 


Tr»r»Mw.n  >■■>.■.»■  *- 


where  Mp  ,  Mp  ,  Mp  ,  Mp  ,  Mp  ,  Mp  ,  Q2x,  Q2y,  Q2g  are  the  components 
X  Y  Z  X  Y  Z 

of  the  vectors  Mp,  MR  and  Q2  in  the  XYZ  system. 

The  vector  Q2,  similarly  as  Qi  in(  2.4),  includes  the  additional 
components,  and  expresses  the  specific  features  of  the  mechanical 
system  under  consideration. 

Let  us  write  the  scalar  expressions  for  the  vector  Q2  in  the  form 
of  individual  groups  of  terms.  The  first  term  is 


* 

2  X  9,  -  2  «A<  X  (g,  +  a,  X  q,)  - 

M)  jsO 

n 

“2  m,  { (q-,z  +  Qixqir  —  Qirqix)  Vet  -  fay  +  %zqix  —  2ixqiz)zJ  /+ 

iiivO 

n 

+  2  mi  [(9/x  +  Qarqtt  —  Qizqiy)2ei  —  (Riz  +  Qixqir—%rqix)xJ  d+ 
/« 0  1 

n 

+  2  mi  [(911'  +  2«9/x  —  %xRiz)  xct  -  (Rix  +  QiYqiz—8izqir)ycl ]  n. 


(3-5) 


<:*o 


In  (3.5)*  it  is  possible  to  make  a  change  of  variables  with  the 
help  of  Equations  (1.1)  and  (1.7)  -  (1.8).  The  second  term  is 


2  2  mAi  X  H  X  g,)  - 

/mo 

n 

2  2  m,  K^xgiy  —  woyqix)  )>ei  +  («o xR'iz  —  <"«9/.Y)*rj!  /  + 

/  s.U 
ft 

+  2  2  Ml  \{"*YthZ  —  Wo Z<Uy)ZcI  +  (<»0Y<llX  —  rf  + 

iwO 

It 

+  2  2  WJ  —  a'0.v9^)  *r  (wa*9/K  —•  *j yRu)  y,,| «. 


(3.6) 


If  it  is  necessary,  one  can  make  a  change  of  variables  also  in  (3.6). 
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The  third  term  is 


2  2  "tAi  X  fa  X  q,). 


By  analogy  with  (3.6),  we  obtain 

n 

2  2  »I  -^yg,x)yct  +  fa*?*  -  «V*7/X)2W1  /  + 

rt 

T  2  2  rnt  ((<* >/y?/2  -  ««?*)  zei  +  fay?,*  _  u^y)  *f/]  + 

It 

+  2  2  m,  \(wl2q,x  -  mlxqlz)xct  +  fa*9/y  _  <alYq,2)ytl\  n. 


The  fourth  term  is 

»  „ 

-  2  3  S  mjt  (»Ay)  (AyiX“()  --22  l/*w«o.w;r  +  /yz,  fazw/z  - 

-  «**«>, y)  -  /<x|HMHfy  +  0,5  (/«,  +  /**,  _  /m)  - 

rt 

-  0,5  ( l xxi  +  Iyyi  -  hz,)  uio^y)  1-22  \lrzmy»ix  +  hxi  X 
X  (Wx  -  «*>, *»*)  -  /xrMr%>  +  0,5(/**,  +  /m  _  /»)  Wjc  _ 

/I 

-  0,5  (/  m  +  /«,  -  / xxi)  «o.wz)  d  -  2  2  ( Woz«>/y  +  /*„  x 

X  (“ov“/y  -  witm*)  -  /yz,u>oz««  +  0,5 (/YVI  +  /„,  /AAl)Wv  _ 
—  0,5  (Izzi  4-  / XXj  _  / yy ()  uJoyu)/AJ  jj. 


The  fifth  term  is 

It 

~~  S  2  mjAjl  X  faXAy,  -f  (v>Ajl)l»,  —  UJjAyj]  — 

It 

“  ~2  Uxxi  fa*  -I-  ~IY'»/Z  —  2/2UI,'}.)  —  (/ yyi  —  IZZI)  miyw,x 

—  1  x\'i  ('“/y  —  w«w*  +  U/zuiyx  —  «,*<«*)  +  /yz/faf*  —  «»;v.)  — 

/» 

-  /*X/faZ  +  «a*(ir  -»•  2,AU>,y  -  G,yW/A)|  l  _  2  |/y,.,(«,y  Q,^/A— 

!*>  1 

—  ~ix^,z)  —  (/ yzi  —  I  xxi )  «•/*«/*  —  Iyzi  (viz  —  w/j-w/.y  +  Qi.v«ir— 


(3.7) 


(3.8) 


(3.9) 
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—Qiyu<x)  +  1  zxi  (*JX  —  ">«)—  /xv/(m+«/y“i2+2/jw/z~2,z.i/r)id  — 

ft 

.  “  2  [/«/  (a'/z  +  -.A'w.y  —  -iy<*ix)  —  (I xxi  —  /  m)«aw/v  -*/  z.vi  X 

•  i 

X  (“/.v  —  w/z^y  *1*  -/v“/z  —  Q/z^r)  +  /xv/  ('“Jy  —  w;*)  — 

—  fyzi  (u J/y  r  *»/z»/.v  t  2/2ui,x  —  S2/*«>,*)J  n. 

By  summing  the  terms  (1)  -  (5)  [see  Equations  (3.5)  -  (3.9)]  using 
the  same  unit  vectors  1,  d,  n,  and  by  including  them  in  (3.4),  we 
obtain  a  system  of  three  differential  scalar  equations  determining 
the  rotational  motion  of  the  system  under  consideration. 

Using  Equations  (3.5)  -  (3-9)»  we  can  make  a  numerical  estimate 
of  each  additional  term  in  Equation  (3.4).  In  order  to  obtain  a  more 
accurate  solution  for  w»y  and  u«g»  the  above  expressions  should 
be  complemented  by  the  expressions  in  (3.4). 

If  the  scalar  equations  are  set  up  with  respect  to  the  principal 
central  axes  of  inertia,  then  of  course  the  structure  of  the  equations 
will  become  simplified.  The  transition  from  the  components  of  the 
vector  b)0  in  one  coordinate  system  to  another  can  be  made  using  the 
matrices  of  direction  cosines. 

Supplement :  The  orientation  of  the  XYZ  system  relative  to  the 

stationary  Cnt  system  is  given  by  the  vector  R_  which  determines 

8 

the  origin  of  XYZ  relative  to  and  by  the  Euler  angles  </,  d,  and 
Y  * 


The  relationship  between  variables  Vx,  Vy,  Vz  and  the  components 
of  V  in  the  system  is  the  usual  one:  it  is  formed  by  means  of 
the  direction  cosines  expressed  in  terms  of  the  angles  ij>,  u,  y,  and 
yields  three  equations.  Since  the  projections  of  the  total  angular 
velocity  vector  of  the  frame  XYZ  are  (ut  +  U«)Y»  («o  +  fto)v,  (“o  +  iio)v» 
respectively,  the  variables  )|/,  0  and  y  must  be  expressed  in  terms  of 
these  projections  using  well-known  relations.  This  shows  the  differ¬ 
ence  between  the  structure  of  the  given  kinematic  angular  Euler  equa¬ 
tions  and  the  analogous  Euler  equations  used  in  the  mechanics  of 
bodies  of  constant  mass:  in  the  latter  the  projects  are  equal  to 

W«X»  u*y»  U#Z* 
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It  will  be  noted  that  by  recalculating  the  projections  a>ox,  woy* 
u>oz  in  the  aoYoZo  system  (attached  to  the  body  "0")»  one  can  use  the 
ordinary  kinematic  angular  equations.  In  this  case,  one  of  course 
assumes  that  the  angles  0,  ip  and  y  specify  the  orientation  of  the  X 
X0Y0Z0  axes  (and  not  XYZ).  Any  of  these  approaches  will  determine 
the  last  three  deficient  equations. 

In  obtaining  Equations  (2.3),  ( 3 • 3 ) »  and  their  scalar  analogs, 
the  vector  quantities  po,  p^,  A#c,  A1(J,  qQ,  qit  qQ,  q^  etc.  were 
used  as  parameters.  However,  p#c,  p#**  w^,  Plc»  and  their  derivatives 
are  of  course  the  starting  vector  quantities  which  are  determined  by 
the  specific  relative  kinematics  of  the  bodies,  distribution  of  mass 
in  them,  rate  of  change  of  the  mass,  and  other  design  factors.  Similar 
remarks  can  be  made  also  with  respect  to  the  moments  of  inertia  of  the 
bodies.  It  is  more  natural  to  consider  the  latter  as  given  functions 
of  time  relative  to  the  X#Y#Zo  system  (for  "0")  and  X^Y^Z^  system  (for 
bodies  "i"),  but  not  relative  to  XYZ.  These  quantities  can  be  easily 
converted  using  well-known  relations. 

Conclusions 


As  a  result  of  the  present  study  carried  out  in  a  fairly  complete 
form  (considering  the  dynamics  of  the  variability  of  mass  characteris¬ 
tics  and  configurations),  we  obtained  a  mathematical  model,  in  parti¬ 
cular,  of  a  flight  vehicle  as  a  control  system  which  can  be  used  in 
the  synthesis  and  development  of  control  systems. 

We  have  shown  the  fundamental  difference  between  the  dynamics  of 
systems  of  variable  geometry  and  mass  and  the  dynamics  of  bodies  of 
constant  mass,  and  this  must  be  taken  into  account  when  designing  a 
number  of  objects. 
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ON  THE  PROBLEM  OF  SYNTHESIZING  SEARCHLESS  SELF-ADJUSTING 
CONTROL  SYSTEMS  FOR  AIRCRAFT 

V.  D.  Yeliseyev 

It  Is  advisable  to  design  self-adjusting  control  systems  by 
minimizing  a  certain  statistical  performance  index.  The  performance 
Index  for  aircraft  should  usually  Include  the  error  determining  the 
accuracy  of  the  system,  and  the  coordinates  characterizing  the  energy 
losses  and  the  accelerations  acting  on  the  craft. 

The  control  system  may  contain  the  following  adaptive  blocks: 
blocks  of  Identification,  blocks  of  self-adjustment,  and  blocks  of 
learning.  After  the  system  Is  used  many  times,  the  blocks  of  learning 
make  it  possible  to  establish  and  use  the  average  functional  relation¬ 
ships  between  the  adjusted  coefficients  and  the  readings  of  the  flight 
regime  sensors. 

In  practical  applications,  the  control  system  may  be  considerably 
simplified.  Thus,  the  learning  block  may  be  absent,  and  the  functional 
relationship  between  the  adjustable  coefficients  and  the  regime  sensor 
readings  may  be  fixed  In  the  design  stage.  To  reduce  the  level  of 
complexity,  the  Identification  blocks  may  also  be  removed.  However, 
the  self-adjusting  blocks  can  be  realized  quite  easily  if  one  uses 
the  gradient  method  in  conjunction  with  simplified  models  of  sensitivity 
[13. 
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Figure  1.  A  scheme  for  optimizing  the 
error  produced  by  the  angular  velo¬ 
city  circuit  for  the  aircraft 
fuselage. 

Wd  —  drive  operator  with  a  rigid 

feedback;  W&  —  aircraft  operator; 

Wang  —  operator  of  the  angular 

velocity  sensor;  L  —  operator  of 
transforming  filter;  Wr  —  operator 

of  the  correcting  filter;  f  —  wind 
disturbance  resulting  in  a  deflection 
of  a  control  surface;  dotted  lines  — 
links  necessary  in  designing  a  non¬ 
linear  sensitivity  model. 


Thus,  in  the  case  of  a 
single-loop  system,  in  a  num¬ 
ber  of  cases  the  sensitivity 
model  may  be  taken  to  have  the 
form  of  a  stationary  filter 
maintaining  an  energy  balance 
between  the  high  and  low  error 
signal  frequencies,  which  is 
proper  for  an  optimal  system. 
The  possibility  of  using  a 
stationary  filter  is  related 
to  the  stability  of  the 
system’s  pass  band  or  adjusted 
values  of  the  regulator 
coefficients. 

A  study  of  a  single-loop 
self-adjusting  system,  using 
a  stationary  filter  as  a 
sensitivity  model,  shows  that 
the  adjusted  value  of  the 
amplification  coefficient 


depends  less  on  the  spectral 

composition  of  the  input  signals  than  in  the  case  of  an  optimal  system. 

We  shall  give  a  simplified  procedure  for  designing  a  single-loop 

self-adjusting  system  with  one  adjustable  coefficient  assuming  that 

the  basic  loop  operates  quasilinearly .  First,  we  find  the  structure 

of  the  optimization  scheme  (Figure  1)  to  calculate  the  current  value 
31 

of  the  gradient  of  the  performance  index  I,  taken,  for  example,  in 
the  form  of  the  absolute  value  of  the  transformed  error  signal 

0/  08t 
du~^:,'Ene" 


3  i 


J 


where  I  ■  e»;  p  is  the  adjusted  parameter;  ei  *  L(p)e;  L(p)  is  the 
transfer  function  (operator)  of  the  transforming  filter;  e  *  x  -  y  is 
the  error  signal,  where  x  is  the  input,  and  y  is  the  output  from  the 
system;  4^-  is  the  sensitivity  function  of  the  signal  e»  to  the  vari¬ 
ation  of  the  parameter  y,  obtained  using  the  sensitivity  model. 

For  a  given  type  of  random  or  determined  signals,  the  optimal 
values  of  the  adjustable  parameter  y  are  calculated  in  all  flight 
regimes  using  the  adjustment  law  ^  *  -A-|^,  where  A  is  a  sufficiently 
small  parameter. 

By  properly  choosing  the  transforming  filter,  we  can  eliminate 
the  low-frequency  components  (if  the  system  is  non-astatic),  change 
the  optimal  value  of  the  coefficient  y  (stability  margin),  and  reduce 
the  effect  of  noise. 

A  part  of  the  optimization  scheme  consisting  of  the  transforming 
filter,  sensitivity  model,  and  the  multiplication  unit  may  be  called 
the  spectral  analyzer,  since  the  signal  produced  by  it  depends  on  the 
ratio  of  the  power  of  high  and  low  frequencies. 

When  constructing  a  real  self-adjusting  system,  the  nonstationary 
model  of  sensitivity  may  be  in  many  cases  "frozen"  at  an  adjusted 
value  of  y  in  a  certain  average  regime  of  flight,  i.e.,  it  may  be 
taken  in  the  form  of  a  stationary  linear  filter.  If  the  accuracy  of 
adjustment  is  insufficient  in  certain  regimes  that  differ  greatly 
from  the  chosen  selected  regime,  it  is  possible  to  use  nonstationary 
filters  by  introducing  a  readjustment  of  the  filter  parameters  in  the 
function  p  (for  example,  in  the  sensitivity  model,  one  can  change  the 
analog  of  the  coefficient  y  and  simultaneously  the  analog 
trol  surface  efficiency  coefficient  which  is  proportional 
Thus,  the  structure  of  a  self-adjusting  system  is  similar 
structure  of  the  optimization  scheme  (Figure  2). 


of  the  con- 

to  it).  * 

to  the 
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Figure  2.  A  scheme  of  a  self-adjusting  system  controlling  the 
angular  velocity  of  an  aircraft. 

1  —  transforming  filter;  2  —  a  model  filter  of  an  optimal  basic 
loop;  3  —  multiplication  unit;  ^  —  filter-amplifier;  5  — 
integrating  unit. 


Figure  3.  An  example  of  a  multi-loop  self-adjusting  system 
involving  an  adjustment  of  the  coefficient  of  each  loop. 

Wi  —  operator  of  the  open  drive  depending,  for  example,  on 
the  hinge  moment  of  the  control  surface;  W2  —  angular  velocity 
operator  of  the  aircraft;  W$  —  operator  which  transforms  the 
angular  velocity  into  either  an  angular  deflection  or  a  linear 
acceleration;  W  _  —  operator  of  the  angle  or  linear  acceleration 
sensor;  1,  2,  3  — »  self-adjustment  blocks. 


After  selecting  the  parameters  of  the  analyzer  of  the  self- 
adjusting  system,  one  must  assure  the  desired  dynamics  of  the  self- 
adjusting  circuit,  in  particular,  one  must  select  a  coefficient  X 
which  can  be  adjusted  depending  on  y  and  on  the  average  value  of  the 
signal  modulus  e». 


It  will  be  noted  that,  to  simplify  the  system,  one  may  sometimes 
use  the  control  surface  deflection  signal,  for  example,  instead  of 
the  error  signal. 
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Figure  4.  An  example  of  a  two-loop  self-adjusting  system  with 
the  adjustment  of  the  damping  loop  coefficient  and  a  functional 
readjustment  of  the  outer  aircraft  control  loop  coefficient. 

1  —  nonlinear  converter;  2  —  inner  loop  self-adjustment  block. 

This  procedure  may  be  extended  to  the  case  of  multi-loop  control 
systems  for  which  it  is  advisable  to  use  a  separate  self-adjusting 
circuit  for  each  loop  (Figure  3)*  Here  the  self-adjusting  circuit  of 
each  loop  uses  the  error  of  its  loop  and  contains  a  simplified  model 
of  the  sensitivity  of  this  loop.  Calculations  of  a  two-loop  system 
controlling  the  angular  motions  of  an  aircraft  show  that  the  adjusted 
value  of  the  internal  circuit  coefficient  is  in  many  cases  practically 
independent  of  the  value  of  the  outer  circuit  coefficient.  In  addi¬ 
tion,  often  the  outer  circuit  coefficient  may  be  readjusted  depending 
on  the  inner  circuit  coefficient,  which  requires  only  one  self- 
adjusting  circuit  (Figure  4). 

Thus,  the  design  of  a  searchless  self-adjusting  system  may  be 
based  on  the  property  of  optimal  systems  which  maintain  the  energy 
balance  of  the  frequency  components  in  the  error  signal  and  can  be 
carried  out  using  the  gradient  method  with  simplified  sensitivity 
models. 
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A  CRITERION  FOR  COMPARING  ADAPTIVE  CONTROL  SYSTEMS 
V.  I.  Kozlov  and  N.  I.  Savkln 

Recently  there  have  appeared  many  different  schemes  of  adaptive 
control  systems,  which  sooner  or  later  have  to  be  compared. 

As  we  know,  a  comparison  of  systems  reveals  the  advantages  of 
some  and  shortcomings  of  others.  The  work  done  to  eliminate  those 
shortcomings  may  bring  the  systems  being  compared  to  the  same  level. 
For  example,  historically  the  systems  with  constant  adjustments  of 
the  regulator  ceased  to  satisfy  the  technical  requirements,  and  were 
replaced  by  parametrically  invariant  adaptive  systems.  Due  to  the 
method  for  designing  high-accuracy  structures,  developed  by  Professor 
Meyerov,  M.  V.,  and  the  rational  physically  realizable  structures, 
developed  by  Professor  Sokolov,  N.  I.,  systems  with  constant  adjust¬ 
ments  have  acquired  "coarse"  properties  with  respect  to  the  change  in 
the  parameters  of  an  object.  In  this  connection,  it  would  be  inter¬ 
esting  to  compare  the  above  adaptive  systems  with  constant  adjustments 
with  the  self-adjusting  systems. 

Unfortunately,  the  experience  accumulated  thus  far  in  the  area 
of  the  appraisal  of  adaptive  control  systems  is  grossly  inadequate, 
due  to  the  complexity  of  the  systems  themselves  and  the  existing 
criteria  for  their  comparison. 
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SKB-3  MAY*  made  a  comparison  of  such  systems  when  designing  an 
adaptive  stabilization  network  for  controlling  the  aircraft  bank 
channel: 

Here  K  =  const  (2)  and  the  time  constant  is  T  *  0.1  -1.5  sec, 

depending  on  the  regime.  The  desirable  low-frequency  transfer  function 

of  the  system  corresponded  in  all  regimes  to  the  aperiodic  block  f^+~\ 

with  the  time  constant  T  «  0.15  sec. 

c 

The  system  was  designed  in  the  constant-adjustment  and  self¬ 
adjustment  versions  using  the  frequency  characteristics  of  the  elements. 
Simultaneously,  an  attempt  was  made  to  find  a  criterion  for  comparing 
adaptive  systems  that  would  be  simultaneously  fairly  simple  and 
objective.  During  the  design,  a  comparison  was  made  of:  a  system 
with  constant  adjustment  and  a  relay  auto-oscillatory  self-adjusting 
system;  a  relay  system  operating  in  the  forced  oscillatory  regime, 
and  a  system  on  the  boundary  of  stability. 

The  designed  adaptive  systems  in  various  versions  satisfied  the 
requirements  placed  on  them.  However,  it  turned  out  that  the  infor¬ 
mation  about  the  frequency  characteristics  of  an  object  (aircraft, 
drive,  and  sensors),  which  is  necessary  in  designing  the  systems, 
varies  from  system  to  system,  and  is  as  follows; 

0  -  100  Hz  —  for  systems  with  constant  adjustment; 

0  -  30  Hz  —  for  a  relay  self-adjusting  system  in  a  forced 
oscillatory  regime; 

0  -  25  Hz  —  for  a  relay  auto-oscillatory  self-adjusting 
system; 

0  -  15  Hz  —  for  a  self-adjusting  system  on  the  boundary  of 
stability. 

"Translator 's  note:  This  designates  the  Special  Design 
office-3  of  the  Moscow  Aviation  Institute. 


The  slightly  larger  amount  of  information  needed  by  a  self- 
adjusting  system  in  the  forced  regime,  as  compared  with  a  relay  auto- 
oscillatory  system,  is  due  to  the  fact  that  the  frequency  of  forced 
oscillations  u>b  of  the  first  system  is  chosen  to  be  larger  than  the 
frequency  of  the  auto-oscillations  to  of  the  second  system.  This  is 

9. 

how  one  achieves  the  greater  stability  under  interference  in  the 
relay  system. 

These  results  lead  us  to  conclude  that  the  amount  of  information 
about  an  object  that  is  needed  for  its  design  may  be  utilized  as  a 
criterion  in  comparing  adaptive  control  systems.  This  criterion  wjll 
enable  us  to  find  systems  whose  design  (design  of  correcting  blocks) 
requires  information  in  the  frequency  range  in  which  the  specific 
features  of  an  object  begin  to  exert  influence.  We  are  talking  about 
features  that  were  not  taken  into  account  by  the  differential  equa¬ 
tions  describing  the  elements  of  the  system. 

For  example,  at  high  frequencies  the  bending  oscillations  of  an 
aircraft  become  important,  and  the  transfer  functions  for  the  object 
as  an  "ideally  rigid  body"  cease  to  be  valid. 

At  high  frequencies  also,  the  transfer  functions  of  the  system 
drive  may  turn  out  to  be  invalid  due  to  additional  lags  introduced  by 
parameters  of  the  drive  that  were  not  considered  before. 

Finally,  there  are  frequency  bands  in  which  the  frequency  charac¬ 
teristics  of  the  drive  —  obtained  for  example  by  means  of  experiments 

—  may  no  longer  be  utilized  in  the  design  of  the  correcting  devices 
using  the  ordinary  methods,  since  the  real  drive  may  exhibit  nonlinear 
behavior.  Of  course,  a  system  —  for  which  the  design  of  the  correc¬ 
ting  devices  was  made  in  the  frequency  range  for  which  one  has  no 
reliable  information  about  the  frequency  characteristics  of  the  drive 

—  may  in  practice  turn  out  to  be  very  sensitive  to  the  neglected 
small  parameters  of  the  drive. 
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Thus,  we  propose  a  criterion  for  comparing  adaptive  control 
systems  in  the  form  of  the  frequency  range  in  which  the  information 
about  an  object,  necessary  for  the  design  of  correcting  devices,  is 
concentrated. 

This  criterion  permits  us  in  the  design  stage  to  exclude  those 
systems  which  require  information  about  the  control  object  in  the 
frequency  region  in  which  that  information  is  practically  inaccessible. 


The  criterion  enables  us  to  see  if  it  is  necessary  to  consider 
certain  small  parameters  of  an  object,  and  gives  limits  within  which 
the  effect  of  those  parameters  must  be  conside. 
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DESIGN  OP  THE  STRUCTURE  OP  A  DIGITAL  SELF-ADJUSTING  SYSTEM 

WITH  A  MODEL  FOR  AN  AIRCRAFT  AND  ITS  IMPLEMENTATION  IN  AN 
ON-BOARD  DIGITAL  COMPUTER 

N.  D.  Litvinov 

Design  of  the  Structure  of  the  Control  Program 
for  the  On-Board  Computer 

Let  us  consider  a  control  system  consisting  of  a  continuous  part 
(composite  object)  and  an  on-board  digital  computer. 

The  continuous  part  will  include  an  interpolator,  usually  of 
zero  order,  an  amplifier  with  the  actuator  (control  surface  machine), 
and  the  control  object  which  has  varying  parameters. 

It  will  be  assumed  that  the  control  object  is  linear  and  has  a 
given  structure. 

Digital  sensors  will  be  used  as  th'  sensors  of  the  data  cn  the 
coordinates  of  the  object  and  the  control  signals. 

The  control  signals  involve  specific  and  random  inputs. 

The  continuous  part  or  the  composite  object  may  be  described  by 
a  linear  difference  equation  with  time-dependent  coefficients. 
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Figure  1.  A  block  diagram  of 
the  control  program. 

1  —  control  program;  2  — 
composite  control  object; 

3  —  optimization  program; 

4  —  program  representing 
the  model  of  the  desired 
closed  system. 


We  shall  consider  quasi-stationary 
processes  for  the  cases  when  the  para¬ 
meters  (coefficients  of  the  equation) 
of  the  composite  object  change  con¬ 
tinuously  and  slowly  or  stepwise, 
and  then  remain  constant. 

It  is  required  to  design  a 
searchless  adaptive  control  program 
for  the  on-board  digital  computer  as 
applied  to  the  composite  linear  non¬ 
stationary  control  object,  that  would 
result  in  the  desired  dynamic  charac¬ 
teristics  of  the  closed  system. 


We  shall  represent  the  control  program  in  the  form  of  a  control 
program  —  generating  the  control  signal  in  the  on-board  computer, 
and  the  optimization  program  —  which  optimizes  the  closed  system  by 
minimizing  the  index  used  to  evaluate  the  performance  of  the  entire 
system  by  adjusting  the  coefficients  of  the  control  program. 

The  control  program  may  be  represented  in  the  form  of  a  block 
diagram  (Figure  1). 

Let  us  determine  the  structure  of  the  control  program. 

The  equation  of  the  control  object  may  be  written  as 


*-i  *~i 

!/(«+*)+  23«/ <«)!/('!+/)-  y,bi{n)  II  (n+ J),  (1) 

i>  u  i  a 


where  yn  is  the  output  from  the  control  object,  uR  is  the  control 
signal  generated  by  the  on-board  computer. 

For  the  basic  regime  of  operation  of  the  control  object,  we  can 
assume  that  the  control  coefficients  a^n)  and  bj(n)  are  constant 
and  given. 
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Let  us  determine  the  structure  of  the  control  program  using  a 
variational  method  [1],  For  this  purpose,  we  obtain  a  difference 
equation  of  the  control  program  from  the  conditions: 

a)  stability  of  the  closed  system; 

b)  astatism  of  first  order  relative  to  external  inputs,  i.e., 
exact  generation  of  a  step  input  in  a  steady  state  motion; 

c )  optimality  of  the  system  in  the  sense  that  the  quadratic 
estimate  of  the  transient  error  process  is  minimized. 

For  this  purpose,  we  write  Equation  (1)  for  the  error  in  the 
closed  system,  and  take  the  first  difference  of  both  sides  of  the 
equation.  Setting  x(n)  *  xol(n),  we  obtain  the  following  equation: 

*1  ■  *-! 

v  £.<■’(*+*)-- Ytyn  («+/), 

t  II  /.  o  It/ 


where 


/?*■=<**-» — «>.  (<!»•>» 1 1  ajs+t^a-i^O)  (s=>  1,2,  1) 


(3) 


A«  (ll)rrtt,  (/,), 


(4) 


in  which  the  initial  conditions  of  Equation  (1)  are  replaced  with 
certain  new  conditions  corresponding  to  the  action  of  the  perturbation 

* 

a,xtoin+l+i), 

/*•  o  1 0 ) 


where 


o  (<i+Y)«AI  («+Y). 


(6) 
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We  represent  Equation  (2)  in  the  form  of  an  equivalent  system  of 
first-order  difference  equations  in  generalized  coordinates  in  terms 
of  the  error  with  certain  new  initial  conditions 

(|t=O.I . ft  +  1).  (7) 


Now  we  can  formulate  the  variational  problem  stating  that  it  is 
necessary  to  determine  the  structure  of  the  control  program  for  which 
the  closed  system,  described  by  the  system  of  first-order  difference 
equations,  moves  from  point  7  to  the  origin 

*„(e°)«0  (8) 

while  minimizing  the  functional 

oo  v*+l  \ 

f  («.)=*  2  [  Sv»i(“)+P«i(#)).  («SiP*=  const),  /q\ 

rt«*0  /  \jJ 


which  represents  a  generalized  quadratic  estimate  of  the  error 
transient  process. 

The  problem  is  variational  and  Involves  a  conditional  extremum. 
However,  it  can  be  easily  reduced  to  finding  an  unconditional  extremum, 
which  —  when  solved  employing  the  usual  methods  —  will  lead  us  to 
an  optimal  control  law  in  the  form 


*+i 

«i  (")=*  2)  Vn  <rt>> 

|l«l 


(10) 


where  are  fully  determined  coefficients  expressed  in  terms  of  the 
coefficients  of  the  equation  for  the  composite  control  object  in  the 
coefficients  and  0,  which  in  turn  are  determined  using  the  method 
of  the  standard  difference  equations  [1],  starting  with  a  given  per¬ 
formance  index  of  the  transient  error  process  in  the  closed  system. 
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Substituting  in  (10)  the  expressions  for  ey(n)  and  Ui(n)  in 
terms  of  e(n)  and  u(n),  we  obtain  a  difference  equation  which  will 
determine  the  structure  of  the  control  program.  The  structure  is 
stable,  since  it  was  determined  from  the  condition  in  (8). 

The  coefficients  in  the  difference  equation  of  the  control 
program  will  be  adjusted  by  the  gradient  method  with  the  aid  of  the 
optimization  program. 

In  order  to  estimate  the  approximation  of  the  desired  process, 
given  by  the  model,  as  compared  with  the  real  process,  we  introduce 
a  performance  index  of  the  adjustment  in  the  form  of  the  functional 

/=»£(£).( p(e),  (11) 

where  P(E)  is  a  certain  difference  operator,  E  is  the  forward  shift 
operator,  $ ( e )  is  the  differentiated  function  of  the  quadratic  form 
in  the  deviation  e(n). 

Fo.  the  gradient  method  of  adjustment,  the  variation  of  the 
coefficients  of  the  control  program  is  giver,  by 

(«)+*/ 0-2) 

where  i  specifies  the  adjusted  parameter,  are  coefficients  which 
generally  depend  on  the  deviation  of  the  position  of  extremum  and 
gradient  components  [2], 

It  will  be  assumed  that  the  coefficients  have  been  selected 
in  such  a  way  that  a  steady  adjustment  process  is  obtained. 

We  obtain  the  following  expression  for  the  gradient  for  the  case 
of  the  above  block  diagram  (Figure  2)* 


if/  (c) 
</<•; 


-£<£) 


<f<l  (e) 
c/e 


</'{£.  *,  c)  (I  —  0  (£,a,  c)J. 


i _ d\v„yi.  c) 

u',a£.  c)‘  dci 


•*(/!), 


"Translator 's  note:  Foreign  text  has  no  Figure  2. 


(13) 
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where  4>(E,  a,  c)  is  the  operator  of  transformation  of  the  closed 
system,  dependent  on  the  adjusted  coefficients  of  the  control  pro¬ 
gram,  c^,  and  the  time-dependent  coefficients  of  the  control  object 
a,  (E,  c)  is  the  operator  of  transformation  for  the  control  program, 
Wob(&,  a)  is  the  operator  of  transformation  for  the  composite  object. 

Since  for  the  deviations  e(n)  equal  to  zero  we  have  *(E,  a,  c)  * 
♦(E) ,  where  *M(E)  is  the  operator  of  transformation  for  the  model, 
the  algorithm  for  adjusting  the  coefficients  c^  of  the  control  program 
may  be  written  as 


a  («+ ‘)-cj  (n)-hP  (£)*5r  **  i  «-♦*  (£>1  WZTeTc)* 


dWn(H,  c) 
dct 


*(«)• 


(14) 


The  algorithm  in  (14)  represents  an  optimization  program,  and 
is  described  by  a  nonlinear  difference  equation. 

To  form  the  algorithm  in  (14),  we  can  use  the  intermediate 
coordinates  of  the  closed  system.  For  this,  we  use  the  relation 

U— 0(C,«.e)l  «(<!)««  <«>.  (15) 


Thus,  the  control  program,  the  optimization  program,  which 
adjusts  coefficients  of  the  control  program,  and  the  program  with  a 
model  of  the  desired  closed  system  will  altogether  constitute  a 
searchless  adaptive  control  program  for  the  on-board  digital  computer. 

A  Method  of  Implementing  the  Control  Program 

In  those  cases  where  a  composite  control  object  is  described 
by  a  high-order  difference  equation,  the  structure  of  the  control 
program  is  also  of  high  order.  From  the  equation  of  the  optimization 
program,  we  can  see  that  its  order  will  be  determined  by  the  orders 
of  the  equations  of  the  model  and  the  control  program.  Of  course, 


15 
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for  a  given  computer  speed  the  elementary  time  interval  of  the  control 
program  will  be  very  high,  which  will  affect  the  accuracy  and  quality 
of  control.  On  the  other  hand,  for  a  given  elementary  time  interval 
of  the  control  program,  we  need  an  on-board  digital  computer  of  very 
high  speed  in  order  to  realize  the  program,  and  this  is  usually 
impossible  to  achieve  in  practice.  Therefore,  other  ways  of  realizing 
the  control  program  are  necessary. 

As  one  of  the  possible  ways  of  realizing  the  control  program, 
one  can  use  the  following  method. 

Due  to  the  fact  that  one  considers  quasi-stationary  processes  in 
which  the  coefficients  of  the  equation  for  a  composite  control  object 
change  continuously  and  slowly  or  stepwise,  and  then  remain  constant, 
by  analogy  with  continuous  self-adjusting  systems  one  can  assume  that 
the  rate  of  the  process  in  a  real  closed  system  is  higher  than  the 
rate  of  the  process  of  adjustment  of  the  control  program  coefficients, 
which  will  be  determined  by  the  length  of  the  range  in  which  the  para¬ 
meters  of  the  composite  control  object  stay  constant.  Therefore,  one 
can  organize  the  operation  of  the  closed  loop  which  includes  the  con¬ 
trol  program  and  the  composite  control  object,  using  the  length  of 
the  elementary  time  interval  selected  from  the  conditions  determined 
by  the  requirements  on  the  accuracy  of  control,  by  the  dynamic  proper¬ 
ties  of  the  composite  control  object,  and  the  frequency  properties  of 
external  inputs.  The  operation  of  the  optimization  program  and  the 
program  of  the  model  of  the  desirable  closed  system  may  be  organized 
for  the  length  of  the  elementary  time  interval  determined  by  the 
length  of  the  interval  in  which  the  parameters  of  the  composite  con¬ 
trol  object  stay  constant.  Of  course,  we  have  to  impose  the  condition 
that  time  intervals  be  multiples  of  the  elementary  time  interval  thus 
obtained,  due  to  the  discreteness  in  the  deviation  e(n). 

When  the  operation  of  the  control  program  is  as  described  above, 
one  will  meet  the  requirements  for  accuracy  of  control,  and  will 
obtain  lower  requirements  on  the  speed  of  an  on-board  digital  computer. 


Conclusion 


In  designing  a  searchless  adaptive  control  program  for  an  on-board 
digital  computer  as  applied  to  a  linear  nonstationary  object,  the 
design  must  be  such  that  the  requirements  on  accuracy  and  quality  of 
control  will  be  met.  The  above  method  of  realizing  the  control  pro¬ 
gram  for  multiples  of  the  elementary  time  interval  enables  us  to 
meet  the  requirements  for  accuracy  and  quality  of  control  in  a  wide 
range  of  the  control  object  parameters. 
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A  DESIGN  OP  A  SELF-ADJUSTING  AUTOMATIC  PILOT 

A.  A.  Kotenko  and  V.  A.  Serdyukov 

The  behavior  of  the  rigid  fuselage  of  an  object  may  be  approxi¬ 
mately  described  by  the  equation 


The  blocks  characterizing  the  elastic  oscillations  of  the  fuselage 
are  represented  by  equations  of  the  type 

vj+aiiV<+a/,V|=ifl/s5;  /*>!, 2.  /~\ 


where 


t 


VlTT 


n 


1*0 


(3) 


We  assume  that  the  dynamics  of  the  control  organ  is  described  by  the 
equation 
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In  Equations  (1),  (2),  (3),  (4)  we  used  the  following  notation: 
v  is  the  pitch  angle  of  the  elastic  object;  v0  is  a  component  of  the 
angle  v  characterizing  the  pitch  deflection  of  the  rigid  fuselage  of 
the  object;  vA  (i  ■  l,  2,  . ..,  n)  are  components  of  the  angle  v  charac¬ 
terizing  the  elastic  oscillations  of  the  object;  a^j  (i  *  1,  1,  ...» 
n,  j  »  0,  1,  2)  are  time-dependent  coefficients;  6  is  the  output  of 
the  control  organ;  6  is  the  input  of  the  control  organ;  b  is  a  con- 
stant  coefficient. 

It  is  assumed  that:  1)  during  flight  the  dynamic  characteristics 
of  the  elastic  object  undergo  changes  (coefficients  a^  change  arbi¬ 
trarily  with  time,  but  slowly  enough  so  that  within  a  given  time  inter¬ 
val  one  can  take  a^  »  const,  which  is  the  condition  of  quasi-statlon- 
arity);  2)  measuring  apparatus  has  zero  inertia;  3)  the  coordinates 
v^,  1  =  0,  1,  . ..,  n  may  be  measured. 

It  is  required  to  design  a  self-adjusting  autopilot  that  will 
compensate  for  the  effect  of  the  changing  parameters  of  a  rigid  object 
and  the  elastic  oscillations  of  the  fuselage. 

The  control  law  is  as  follows 

6 e~(D+b)  ^(*/iD+*/,)  v/j  (  5  ) 


Here  g  is  the  control  signal,  kj#,  kj j ,  ks  are  adjustable 

coefficients. 

We  substitute  (5)  in  (1)  end  (2),  taking  (3)  and  (4)  into  account, 
and  obtain 


n  n 

V/ + o„V| +OJ.V)  •»  M/I«— 23  'V— 5  “i  »*/«v/  (f^-o.l . n.)  (g) 

Jtm.0  jnO 

We  approximate  the  desired  dynamic  characteristics  of  the  basic  loop 
of  the  system  ''object-autopilot”  by  a  standard  filter  (model)  of  the 
form 
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Here  aj^  (i  *  0,  1,  2)  are  constant  coefficients  of  the  model. 

We  set 

e“v~v"=i,v'~v*  (8) 

We  shall  sum  Equations  (6)  and  subtract  (7)  from  this  sum.  In  view 
of  (8),  we  obtain 


e  H-«,  t+a,t 


-3  S  ( ail—ai+na„ku)\ 

/t*»0  /i-0 


(9) 


We  assume  that  the  coefficients  of  the  control  law  (5)  are  adjusted 
by  the  algorithm 


{m0,|,.,.,n;  /«=- 0,1. 


We  represented  (9)  and  (10)  in  matrix  form 


e  — li, 

k 


(10) 


(11) 


where 


A~ 


m  _  \ 

\t—Q 


0 

I.  " 

k-  — «*!+«<*/**//)  v)y) 

1  /«U/~0 


/ 


(a 


♦«. 

*1# 

fi. 

*•1 

;  ♦= 

?.v 

*8 

(12) 


The  adaptation  algorithm  (10)  will  be  formed  using  the  condition  of 
the  existence  of  the  Lyapunov  function  for  System  (11).  The  Lyapunov 
function  is  assumed  to  have  the  form 

V  xe' Pe + ).  V)  £  (ail- “) + ««/t *//)*  +  X,  Y  (auk, -a, )» 
l-o  i-o  So 


(13) 


*  1'VWJWW 


Here  tc,  X2  and  X  are  certain  coefficients,  and  P  is  a  square  positive 
symmetric  2x2  matrix;  e'  is  a  transposed  matrix. 

In  view  of  Equations  (11)  and  the  condition  that  the  coefficients 
ai3  »  ■  ».  1,  ...»  n;  j  ■  0,  1,  2)  be  quasi-st at ionary,  the  derivative 
of  (13)  will  be 

V  =ixe'Me+2x  (P}ie+PfSe)^2)o/t*,— ojj  g— 

jS (0//~ «/+"a(«fci/),ui/>j  +2X»  (jS^1***-*)  (jS,0'*)  ?5+  (14 ) 

/i*0  io-0  \/*0  / 

Here  M  is  a  2  x  2  negative  matrix.  The  elements  of  the  matrix  M  which 
are  specified  arbitrarily  are  related  to  the  elements  of  the  matrix  P 
in  (13)  by  the  equation 

A'P+PA**M.  (15) 

If  the  real  parts  of  the  roots  of  the  characteristic  equation 


M— Xfi!»  0  (lb. 

are  negative,  then  according  to  [2]  to  the  definite  negative  form 
e'Mc  in  (14)  there  always  corresponds  a  definite  positive  form  e'Pe 
in  (13) .  Consequently,  in  order  for  the  adaptation  process  in  the 
"object-autopilot"  system  to  be  stable,  it  is  sufficient  that  the 
following  conditions  be  satisfied 

* (P:.«  +  P«J*) | ^ -< 0«  j  K-  V  naltk,i)  j  + 

+2).»S(«/l*,-5)fs«.-,)  V,+2ky,y  (a, /-ai+na^k.Dx  (1?: 

\rb  /  "0“ 

VU<0 

Since  changes  in  the  deviations  of  the  coefficients  are  independent 


/—a/— ««/-*//),  (>]<»/,*,— o,  ),  (<--0.  I,  ...  n\  j™ 0,1) 

\  /r-0  / 
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we  have  the  following  sufficient  conditions  for  the  inequality  in 
(17)  to  be  satisfied 


w« 0  / 

x  iPue+Ptti)vj/l'=i  ■*)**. 

(**0,  1,  ....n;  /«0,  1. 


(13) 


We  take 


*-“TT — • 

S«a« 


Then  the  algorithms  for  the  adjustment  of  the  coefficients  of  the 
control  law  (5)  become 

k.«=-x  (Pji® + Put)  g, 

K|/-x (J>„B+P,te) oj«,  /=0,  1 . «;  /— 0.  I.  (2°) 

In  (20)  the  factor  k  may  be  arbitrary  (including  equal  to  infinity). 

In  this  case  with  the  coordinate  or  K2  restricted  in  absolute  value, 
we  obtain  a  relay  switching  of  the  coefficients  of  the  control  law 


Ki=  ~B  sgn  (Ptit+Pni)  sc n  g 
A,‘//«=>fisgn(i,,le+P«ie)sgn o}*’  i=»0,  l . .  /«•(),  I. 


(21) 


According  to  the  design  of  the  algorithms  (20)  and  (21),  the  self- 
adjusting  loops  for  the  system  "object-autopilot,”  designed  according 
to  (20)  or  (21),  assure  the  stability  of  motion  of  the  output  relative 
to  the  output  of  the  standard  model,  and  consequently,  assure  the 
invariance  to  within  e  of  the  dynamic  characteristics  of  the  system 
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"object-autopilot"  with  respect  to  the  changing  parameters  of  the 
object  and  the  elastic  oscillations  of  the  fuselage. 
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ADAPTATION  OP  SYSTEMS  WITH  CONSTANT  ADJUSTMENTS 


M.  I.  Savkin 

We  know  of  adaptive  systems  in  which  a  change  in  the  parameters 
of  the  control  object  is  compensated  for  by  a  corresponding  change 
in  the  adjustment  of  the  regulator.  If  a  system  contains  a  nonsta¬ 
tionary  object  and  it  includes  special  correcting  devices  with  con¬ 
stant  parameters,  one  can  reduce  the  range  of  its  performance  indices. 
As  a  result,  such  systems  with  constant  adjustment  acquire  properties 
that  are  analogous  to  adaptation. 

This  will  be  shown  using  an  example  of  a  simple  system.  Here, 
we  assume  that  conditions  are  satisfied  so  that  the  system  may  be 
regarded  as  quasi-stationary . 


«Mp) 


Suppose  we  are  given  the  transfer  function  of  the  object 

Ku'  ,  the  amplification  coefficient  *o 11 ,  and  the  time  constant 
r^'p+i 


that  satisfy  the  following  inequalities  in  the  design  regimes  (1,  II) 


«5,M»i;  *o'»Ko;  *0 


Suppose  that  we  must  select  the  structure  and  parameters  of  a 

sequential  correcting  device  WK(p),  such  that  the  cutoff  frequency 
I  II  I  II  A 

w  *  of  the  system  *  *  (p)  closed  by  a  single  negative  feedback, 

will  be  greater  than  or  equal  to  a  certain  desired  value  u>d,  where 


kJ  kJ1 

7T  <«k<*57T* 

,U  *0 


If  WK(p)  »  K  is  taken  such  that  »  w,,  then  we  obtain 

7\i  G 


(1) 


where 


xxi' 


7?ir 


In  this  case 


(2) 


Here  the  performance  indices  of  system  (1)  will  differ  greatly  depen¬ 
ding  on  the  design  regime. 


However,  if 


then 


K(Ttp+ 1). 
7>+l  • 


K-K'0 

tI 

Jo 


1 7",  =»  o.' 


a  ' 
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and 


Pu»  Tj  1 

«'  “  r>> 


By  properly  selecting  the  value  of  a  (for  example,  based  on  the  per¬ 
missible  stability  margin),  one  can  bring  the  cutoff  frequencies  w1 
and  w11  closer  by  a  factor  of  o,  and  thus  achieve  the  same  effect  in 
the  range  of  the  performance  indices  of  the  system.  This  effect  may 
be  amplified  if  one  takes  a  "closing”  correcting  device 


(ry+^,r,p+i)"' 


The  example  considered  above  is  illustrated  in  Figure  1.  In  conclu¬ 
sion,  we  note  that,  if  at  the  input  to  the  system  one  places  a  low- 

frequency  filter  with  constant  para- 


a -d3. 4> (deg.  ) 


meters 


"  \\ 
v5 VC 

10  ,T> 


WHl’ 


i  \  \ 

'  \  V. 


V  *'■* 

,  '>A~ 


- f-twf  8> v Ml.h/sec fhere  T  s  <2  "  5)“d  «  then  the 

>.  'v  system  thus  obtained  will  have 

s  ju  \  v  \  \ 

\  50  s\  *\  approximately  constant  dynamic 

\s  \\  characteristics.  It  is  important  to 

»»  _  ' - note  (see  the  preceding  article) 

u0  V*v  \  /  that  such  a  system  with  "closing" 

m  correcting  devices  requires  an 

iso  identification  of  the  object  W0(p) 

m  in  a  smaller  frequency  range  than  a 

system  with  fixed  adjustments  without 

Figure  1.  Logarithmic  frequency 
response  characteristic  of  the  such  devices, 
system. 

Thus  we  have  considered  a  way  of  designing  a  single-loop  system 
with  constant  adjustments  that  possesses  a  property  similar  to 
adaptation.  This  method  can  be  applied  to  stabilize  the  dynamic 
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properties  of  multi-loop  systems  with  constant  adjustments.  In  this 
case,  the  "closing"  correcting  devices  must  be  put  in  each  loop  of 
the  system  that  includes  given  elements  with  variable  parameters. 
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